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ABSTRACT 


This report describes a user oriented computer program SISCAP (SIS 
Computer Analysis Program) for analyzing SIS mixers. The program allows 
arbitrary impedance terminations to be specified at all LO harmonics and 
sideband frequencies. It is therefore able to treat a much more general 
class of SIS mixers than the widely used three-frequency analysis, for 
which the harmonics are assumed to be short-circuited. An additional 
program, GETCHI, provides the necessary input data to program SISCAP. 

The program SISCAP performs a nonlinear analysis to determine the SIS 
junction voltage waveform produced by the local oscillator. The quantum 
theory of mixing is used in its most general form, treating the large 
signal properties of the mixer in the time domain. A small signal linear 
analysis is then used to find the conversion loss and port impedances. 

The noise analysis, includes thermal noise from the termination 
resistances and shot noise from the periodic LO current. Quantum noise is 
not considered. 

Many aspects of the program have been adequately verified and found 
accurate. However, this project was terminated before the program was 
extensively employed, and it is possible that errors may occur in various 
untested parameter ranges. 

An example is given, which forms part of some initial investigations 
into the effect of decreasing the SIS junction capacitance on the 
performance of a simple SIS mixer circuit. Early indications are that, 
for this circuit , a high WpR^Cj product is desirable. 
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LIST OF SYMBOLS USED IN THIS REPORT 


B bandwidth 

0 thermal energy parameter (1/kT) 

Cj SIS parasitic junction capacitance 

A energy gap parameter 

6 Kronecker delta function 

61 small signal current flowing in SIS junction 

61' small signal current in augmented small signal circuit 

6V small signal voltage across SIS junction 

6i thermal noise current 

T 

6i shot noise current 

S 

6t' time increment 

6V voltage increment 

e electronic charge 

f cycle frequency 

h Planck's constant divided by 2tt 

Irk Kraraers-Kronig function 

^DC DC SIS tunnel junction current flowing in absence of LO 

I^L nonlinear SIS large signal current 

iLIN linear large signal current 

iLIN nth Fourier harmonic of I^IN 

n 

J SIS small signal response function 

k Boltzmann's constant 

Lij conversion loss to i^-th sideband from the j-th sideband 
m index parameter 
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m frequency index 

M numerical accuracy parameter 

P convergence parameter for voltage update method 

Pq noise power out 

<|> barrier phase 

Q numerical accuracy parameter 

Rjq identity element resistance for voltage update method 
Rfl normal state SIS resistance 

t SIS normalized time response parameter 

t time 

T SIS physical temperature 

Tf4 mixer input noise temperature 

Tdsb mixer input double sideband noise temperature 

Tggg mixer input single sideband noise temperature 

U large signal SIS time domain response variable 

V SIS tunnel junction voltage 

updated SIS tunnel junction voltage 
V n nth Fourier harmonic of V 

noise output voltage 
Vro LO source voltage 

Vqc DC source voltage 

Vgap SIS gap voltage 

^TOTAL total voltage range required of discretized I-V curve 
Wrq small signal analysis variable 

o)p angular LO frequency 

a) Fourier transform variable 
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m-th sideband frequency 
X SIS junction impulse response 

Y small signal intrinsic admittance matrix 

Y' small signal augmented admittance matrix 

Z' small signal augmented impedance matrix 

Z e embedding impedance 

Z c impedance of the SIS junction capacitance Cj 

Zq transmission line impedance for multiple reflection technique 

Z m input impedance of the mixer at the m-th sideband 

Z' center row of small signal impedance matrix 
0 
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2. Introduction 


As a result of work by Tucker [1] on the quantum theory of mixing, it 
is now possible to predict the performance of SIS (superconductor- 
insulator-superconductor) quasi-particle mixers with reasonable accuracy. 
This theory has been quantitatively verified in its three frequency 
approximate form [2] . The purpose of this report is to make available a 
user oriented computer program SISCAP (SIS Computer Analysis Program) for 
determining the performance of an SIS mixer with arbitrary terminations at 
the LO (local oscillator) harmonics and harmonic sidebands. The program 
uses the Tucker theory in its most general form. It is much more 
versatile than earlier analyses which employed the three-frequency 
approximation, for which all LO harmonics above the fundamental and 
sidebands above the upper sideband of the LO fundamental are assumed 
short-circuited. 

The general outline of this report follows that of a similar work by 
Siegel and Kerr [3] , which described a computer program to analyze 
Schottky diode mixers using classical mixer theory. Reference [3] has 
proved very valuable in the analysis of Schottky mixers, and it is hoped 
that this report will serve the same purpose for SIS mixers. 

The mixer analysis program SISCAP is subdivided into three parts. 

The most difficult part of the program is the nonlinear large signal 
analysis which is used to determine the LO voltage waveform at the 
junction. This analysis is done iteratively to find a self-consistent 
solution and requires a significant amount of CPU time. Following the 
nonlinear analysis, a linear small signal analysis computes the conversion 
loss between the input (signal) and output (IF) ports as well as the input 
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and output impedance match. Finally, the program determines the 
equivalent input noise temperature of the SIS mixer. As an additional 
feature, this computer program performs the three frequency model analysis 
if desired. The run time for the example listed in Section A. 3 is 0.5 
mins, on an AMDAHL V6. 

Although this report is written explicitly for a mixer using single 
SIS junctions, the program may be used to analyze any single-particle 
tunneling device in the quantum regime whose I-V curve is resistive at 
large voltages. In particular, the program reduces a series array of SIS 
junctions to its single junction equivalent by assuming the junctions are 
identical. In any situation however, the DC I-V curve must always be a 
monotonic function and the user should ensure that experimental 
measurement errors do not violate this requirement. 

This report also describes the program GETCHI, which must be run 
prior to SISCAP. GETCHI processes an experimental SIS DC I-V curve to 
provide input data for SISCAP in proper form, and determines the value of 
various control parameters for SISCAP. Further information on this 
program is given in Section 3. An empirical study of the digitization 
constraints required to maintain numerical errors to an acceptably low 
value is also described. 

1.1 SIS Junctions : Definition of Terms 

This work is not meant to be a general reference on the theory and 
operation of SIS mixers. It is assumed that the reader has a solid 
working knowledge of this topic, and explicit detail is given only where 
necessary to describe the computer program and the choices made in its 
design. For general reference, the reader is referred to a recent review 
paper [4] . 
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Fig. 1; Typical SIS DC I-V Curve 
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This section outlines the terms and symbols to be used in this 
report. The units used here are SI units. A typical quasiparticle DC I-V 
curve of an SIS junction is shown in Fig. 1. Note that, for this program, 
superconducting pair tunneling currents are ignored. The rapid rise in 
current seen in Fig. 1 occurs at the gap voltage Vg ap , which is related to 
the energy gap parameter A by the formula: 


'gap 


2A 

e 


( 1 ) 


At large DC bias voltages, the I-V curve approaches the resistance Rm, 
which is also the resistance of the junction in the normal state. The SIS 
junction has the geometrical form of a parallel plate capacitor, whose 
capacitance Cj acts to shunt the nonlinear quasiparticle current. 

In order to describe mathematically the behaviour of the SIS 
junction, the time domain equations given by Tucker [1] will be used: 


V(t) t 

xNL(t) - + Im (U*(t) / x(t-t') U( t ' ) dt'} 

R N 


( 2 ) 


2 “ . ,fK0. 

X(t) = - / [Idc(-) “ — 

7r 0 e eR N 


] sinoit da) 


(3) 


U(t) = exp{i<f>(t)} 


(A) 


e c 

<f>(t) = - - / V(t ' )dt ' 

Tl 


(5) 


Eq. (1) gives the instantaneous quasiparticle current iNL(t) as a function 
of the junction voltage V at all time t. The impulse response function 
x(t) of the SIS junction is related to the nonlinear DC I-V curve through 
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the sine Fourier transform as shown in Eq. (3). U(t) represents the 
driving force applied to the -junction and is related via the barrier phase 
<f>(t) to the instantaneous voltage V(t) in Eqs. (4) and (5). 

For the small signal analysis, the SIS junction is treated in the 
frequency domain. The general frequency domain I-V relation is the 
Fourier transform of Eq. (2), which is complex. The appropriate small 
signal response function is then: 

J(v) = I KK (v) + i Ipc( v ) (6) 


where 


V 

[idc(v) - -] 
1 “ % 
t kk( v ) = p - / 

^ —oo V - V 


dv 


(7) 


Finally, the problem of analyzing series arrays must be considered. A 

2 

property of arrays is that their saturation power level increases as N , 
where N is the number of junctions in the array. The conversion 
gain and noise temperature are independent of N, both in theory [4] and 
apparently in experiment [5] . Because of these useful properties, 
the program described in this report has been written to analyze such 
arrays. The array is reduced to its single junction equivalent on the 
assumption of identical junctions, by dividing both the current and the 
voltage of the array DC I-V curve by N [5] . 
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2. Outline of the Necessary Theory 

Given the bias voltage, LO amplitude and the DC I-V curve, an SIS 
mixer analysis should output the following quantities: conversion loss, 

input impedance, output impedance and mixer noise temperature. These 
quantities will be calculated by a small signal analysis once the large 
signal voltage and current waveforms at the SIS junction are known. The 
first step in the mixer analysis algorithm therefore is to determine these 
waveforms . 

2. 1 Large Signal Analysis 

The most difficult part of this project is the large signal nonlinear 
analysis. All of the Fourier coefficients of the LO waveform must be 
determined in a self-consistent manner. Several papers have appeared in 
the literature on this problem in the context of the Schottky diode mixer. 
Early workers such as Torrey and Whitmer [6] assumed the voltage across 
the diode junction to be sinusoidal, with the implication that the 
harmonics of the local oscillator were short circuited. Fieri and Cohen 
[7] removed this assumption but their work was restricted to simple lumped 
element embedding networks. More realistic embedding networks were 
analyzed by Egami [8] and by Gwarek [9], using harmonic balance 
techniques. However, these schemes frequently had difficulty converging 
for some circuits of practical interest. In particular, the method of 
Egami [8] tended to diverge when more than three harmonic terms were 
considered, and the convergence of the method of Gwarek [9] was a strong 
function of the initial solution estimates. 

In a bid to overcome these convergence deficiencies, Kerr [10] 
developed the multiple reflection technique, in which a 
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hypothetical lossless transmission line of arbitrary characteristic 
impedance is introduced between the nonlinear diode and the linear circuit 
elements. This method converges in all cases tested and the convergence 
rate [11] depends upon how close the harmonic embedding impedances are to 
the characteristic impedance of the transmission line. There is no 
dependence upon one's estimate of the initial conditions. 

Another technique, developed by Hicks and Khan [11,12], consists of 
two dual methods, the voltage update and the current update methods. This 
has also converged for all cases tested, and the rate of convergence 
depends upon the proximity of the embedding impedances at each LO harmonic 
to either short circuits or open circuits. In both this technique and the 
multiple reflection method, the algorithm operates in the time domain when 
considering the nonlinearity and in the frequency domain when dealing with 
the linear embedding network. 

Most practical SIS mixers have rather large -junction capacitances 
(WpRNCj 1 products are usually in the range 1 to 10), and so the embedding 
impedance seen by the nonlinear junction conductance approaches a short 
circuit with increasing frequency. This condition therefore lends itself 
to the use of the voltage update method of Hicks and Khan [11] and this 
method has been adopted for use in the program SISCAP. Such a choice will 
lead to a rapid convergence rate. Moreover, the use of the voltage update 
method permits the nonlinear SIS equations to be solved in their simplest 
form. This is a voltage-input-current-output mode, and requires only a 
direct calculation from the SIS equations. The multiple reflection 
technique requires the SIS nonlinear equations to be solved on a 

= angular pump frequency, R^j = normal-state resistance of the 
SIS device and Cj = junction capacitance. 
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current-input-voltage-output basis, which adds an extra iteration loop to 
the process. To improve convergence further, the voltage update program 
can incorporate a parallel identity resistance technique as outlined in 
[11]. Typically, identity resistances of the order of 0.5 Rjj have been 
found to give satisfactory performance. 

Finally, for additional flexibility, the multiple reflection method 
[10] has also been added to this program as an additional option in the 
large signal analysis section and should be used where low WpRjqCj products 
(less than 0.5) are encountered. In such cases, considerable experience 
is necessary for correct use of the voltage update method and the 
relatively "hands free" multiple reflection technique is preferable. 

2*1.1 Solution of the Large Signal Problem 

The equivalent circuit of the mixer is shown in Fig. 2. For analysis 
of this mixer, the circuit is bisected at the linear-nonlinear interface 
with each half treated separately. The -junction capacitance, Cj, is 
treated as being part of the linear circuit. 

2 

The expressions governing the behaviour of the intrinsic SIS 
junction are given by Tucker [1] in the time domain as follows: 

V(t) t 

iNL(t) = + im (u*(t) / x(t-t') U(t') dt'} (2) 

R N 

2 00 hw hu) 

X(t) = - / [Idc(~) “ ] sinmt du (3) 

* 0 e eR N 

U(t) = exp{i<f>(t)} (4) 

2 

The expression intrinsic junction is used in this paper to refer to the 
nonlinear tunneling admittance of the junction, considered separate from 
the junction capacitance. 
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U- LINEAR NETWORK -H JUNCTION 


Fig, 2: SIS Mixer Equivalent Circuit 
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e r t 

♦ (t) - - - / V(t ’ )dt ' (5) 

ft 


where Iqq (V) Is the measured DC I-V characteristic of the SIS 

junction with no LO power applied; 

Rfl is the normal resistance of the junction; 

V(t) is the instantaneous voltage across the junction; and 

iNL(t) is the instantaneous quasiparticle tunneling 

junction. 

The conditions imposed on the steady state response of the mixer 
circuit by the embedding network can be more conveniently expressed in the 
frequency domain. Referring to Fig. 2: 

V 

n 

” = [z e (naip) // Z c (nu>p) ] n — 2,3,.... (8) 

LIN 

I 

n 

V - V 
1 LO 

= [Z e (u p ) // Z c (w p >] (9) 

LIN 

I 

1 


V - V 
0 DC 

= [Z e (0) ] (10) 

LIN 

I 

0 

LIN 

where V n and I n are the amplitudes of the Fourier components of V(t) and 
iLIN(t) a t frequency nojp, V^q and Vqq are the amplitudes of the Thevenin 
equivalent LO and DC voltage sources, and Z e (no) p ) is the impedance of the 
equivalent external embedding circuit at frequency ntOp, and Z c (nidp) = 
-i/nUpCj. 
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2.1.2 Voltage Update Algorithm 

The iterative voltage update algorithm [11] works as follows: 

(i) An initial estimate is made of the large signal steady state 
voltage waveform across the nonlinear tunnel junction, V(t); 

(ii) Eqs. (2) - (5) are used to determine the current iNL(t) which 
is produced by the voltage V(t) across the intrinsic junction; 

(iii) Using Kirchhoff's current law, I LIN (t) = -I NL (t); 

(iv) A fast Fourier transform gives I L ^(m) from I^IN(t); 

(v) This current, which flows into the embedding network, 

must be maintained by a voltage V^(m) whose various Fourier 
components are given by Eqs. (8) - (10); 

(vi) Using an inverse fast Fourier transform, V^(t) is obtained from 

Vt(m); 

(vii) V^(t) and V(t) are compared; if "equal”, the iteration is 
complete. If not, a new V(t), equal to pV^(t) + (1-p) V(t), 
is defined, and the cycle repeated from step (ii). The 
convergence parameter, p, is normally a value in the range 

0 < p < 1 [11] . 

This procedure is called the "voltage update method", because the 
mechanism for approaching the self-consistent solution in a controllable 
fashion is the updating of the junction voltage waveform with each 
iteration in step (vii). It may be shown [11] that this method converges 
when, at all the harmonics of interest, the magnitude of the impedance of 
the linear embedding network is smaller than the magnitude of the 
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impedance of the nonlinear circuit. This is the case in an SIS mixer by 
virtue of the junction capacitance, considered here as part of the 
embedding network. The value of the parameter, p, is chosen on the basis 
of previous experience [11]. Essentially, its value represents a tradeoff 
between speed of convergence and the risk of divergence. In this work, p 
was fixed at unity since the use of an identity element (to be discussed 
below) proves sufficient to guarantee convergence. 

The computer program allows parallel resistive identity elements [11] 
to be inserted at the linear-nonlinear interface. A resistive identity 
element consists of the parallel combination of a resistance Rjd and a 
resistance -Rid* The net effect of the two parallel resistances on the 
circuit performance is zero since one cancels the effect of the other. 
However, in the large signal analysis, convergence is improved since the 
resistance Riq is lumped in parallel with the linear embedding network and 
in a similar manner, the resistance -Rjd I s added in parallel with the 
nonlinear circuit [11]. The effect of such an addition is to enhance the 
voltage update properties of the circuit by increasing the effective 
impedance of the tunnel junction and simultaneously decreasing the input 
impedance of the embedding network. From experience, setting Ri D equal to 
0.5Rfl has been found effective. 

2.1.3 Multiple Reflection Technique 

The multiple reflection algorithm [10] provides an alternative to the 
voltage update technique discussed in the previous section. In this 
method, a lossless transmission line of characteristic impedance, Zq, is 
inserted at the linear-nonlinear interface. Its length is an integral 
number of LO wavelengths. The steady state operation of the modified 



13 


network is identical to that of the original circuit, but the approach to 
steady state differs from that of the original circuit. Also, because the 
solution approaches steady state in a physically realizable way, the 
multiple reflection method should not be susceptible to numerical 
instability. 

The iterative procedure proceeds as follows: 

(i) A right propagating wave from a source impedance Zq is emitted 
from the embedding network and impinges on the SIS device. From the 
nonlinear circuit equations, the periodic voltage and current produced by 
the incident wave at the SIS device terminals may be calculated. 

(ii) The nonlinear SIS junction generates a left propagating wave 
containing, in general, all harmonics of the LO. After a delay, the left 
propagating wave impinges on the embedding network. 

(iii) From the reflection coefficient of the embedding network, a new 
right propagating wave may be deduced. 

This iterative cycle of steps (i) - (iii) continues until convergence 
is reached. At convergence, the SIS junction voltage waveform remains 
unchanged with increasing iteration number. This method should converge 
for any value of Zq, the transmission line characteristic impedance, 
although the rate of convergence is affected. Typically, Zq equal to % 
has been found to give good results. 

2.2 Small Signal Analysis 

The small signal analysis follows precisely the theory as formulated 
by Held and Kerr [13] for the analysis of Schottky diode mixers, with the 
conversion admittance matrix elements given by Tucker [1], The small 
signal and noise theory given here is based largely on reference [14]. 
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Using the results of the large signal analysis together with a knowledge 
of the embedding impedances at the various sideband frequencies, a linear 
small signal conversion admittance matrix may be formulated for the mixer. 
This matrix gives the relationships between the small signal sideband 
currents and voltages at the nonlinear element, and hence gives all the 
small signal properties of the mixer. 

2.2.1 Frequency and Subscript Notation 

If a mixer is pumped at frequency top and has an intermediate 
frequency oiq, then, under the linear assumption, the only small signals 
which can produce an IF response are at the sideband frequencies 
tog + mojp, m = 0, ±1, ±2,... . Following Saleh [15] it is useful to define 
the sideband frequencies by: 

(Ojj = cog + m(0p m = 0, ±1, ±2, ±3,... (11) 

It may be seen that lower sideband frequencies are represented by negative 
terms (m < 0) . A brief comment on the meaning of these negative frequency 

3 

terms is given in the footnote. 

Saleh’s frequency notation leads to a considerable simplification of 
the mathematics of mixer theory. Using this notation all upper sideband 
frequencies (<00 + ]m|( 0 p) are considered positive, while all lower sideband 
frequencies ( coq - |m|< 0 p) are negative. The sideband frequency index m is 

3 

Electrical quantities are frequently described by a single complex 
quantity associated with some frequency, assumed positive. For example, 
a voltage of frequency to may be described simply by its complex half- 
amplitude V, implying an instantaneous voltage v(t) = V eJ a)t + V* e~i wt . 
It is just as meaningful to work with a negative frequency (-to) and the 
conjugate of the complex half-amplitude (V*), provided the convention is 
clearly understood. Impedances and admittances are then simply the 
conjugates of their conventional positive frequency values, i.e. 

Z(-co) = V*/I* = Z*(co). 




15 


used as a subscript to the various electrical quantities and hence the 
upper sideband, intermediate, and lower sideband frequencies are: 

0 ) + i = cog + Wp, cog, and = oig - a) p ; and V+j , Vg and V_i represent 
voltages at these frequencies. 

2.2.2 Small Signal Conversion Matrix 

Using the sideband notation described in the previous section, let 61 
and SV denote the vectors of the small signal sideband currents (6l n ) and 
voltages (6v n ) at the terminals of the intrinsic SIS junction. Then 

fil = [..., 5lj, 6I 0 , 6I_! , ...] T (12) 

and 

6V = [..., fiVj, 6Vg, 6V_i, ...] T . (13) 

/n. /x 

Torrey and Whitmer [6] have shown that 61 and 6V are related via a 
conversion admittance matrix Y defined by 

61 = Y 6V . (14) 

If the row and column numbering of the square matrix Y correspond with the 
sideband numbering, Y can be written out as: 



with element values given by [1]: 


- ( - i) 


2ho) m i n,n'=-“ 


I ^LoCntOp) WLo( n 't)p) 5 m _ m t >n »- n 


{ [ J(no) p + eVo/h) - J(mu p - (%i + eVo/h)] 

- [j*(n'o) p + (%. + eVg/h) - J*(n'w p + eV 0 /h)]} 


(16) 


where W^g is defined by: 


e t °° -inuipt 

exp{-i - / [V(t') - Vg] dt'} = l W L0 (na) p )e ‘ 

“ — oo n== — co 


(17) 


V(t) is the instantaneous large signal voltage across the tunnel 


junction; 
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6 is the Kronecker delta function; 

J is defined in Section 2.2.5; 

Vq is the DC component of the junction waveform; and 
* denotes the complex conjugate operation. 

The matrix Y can be regarded as the admittance matrix of a 
multifrequency multiport network, as shown in Fig. 3, in which there is 
one port for every sideband frequency . If the parallel combination of 
the embedding impedances and the parasitic junction capacitive impedance, 

Z e / / Z c , corresponding to the sideband frequencies, a^, are now 
m m 

connected in parallel with the intrinsic junction, an augmented 
network is formed as shown by the broken line in Fig. 3. The ports of 
the augmented network correspond to the terminals of the intrinsic SIS 
junction at the various sideband frequencies and do not represent 
physically accessible ports in the real mixer. The augmented network can 
be described by the admittance matrix Y' , defined by: 


61' = Y' 6V (18) 


where 


61' = [..., 6IJ , 6r , 51' , ...] T 



18 


and 

/v 

6 V = [..., 6 V ls 6 Vq , 6 V_!, ...] T , (19) 

6 V m and 61' are the small signal voltage and current components, 
m 

at sideband <^= 0 ) 0 + majp (port m) of the augmented network. The 
elements of the augmented admittance matrix Y' are given by: 


Y' = Y m * n 

mn mn 

and 

Y' = Y + [Z e // Z c ] -1 , m = n 
mm mm m m 


Inverting (18) gives 


6 V = Z' 61' , 

where 


( 20 a) 


(20b) 


( 21 ) 


Z’ = (Y ') _1 . (22) 

The impedance matrix Z' enables us to calculate the conversion loss and 
the input and output impedances of the mixer and is also needed in 
computing the noise temperature. 
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2»2»3 Mixer Port Impedances 

The impedance ^m of any port of the intrinsic SIS junction (see Fig. 
3) can be found by open circuiting the corresponding embedding impedance 

i'V 

Z e and then forming the Z f matrix defined by (21). The desired port 
impedance is given by the mm-th element of the newly formed Z’ matrix: 

Z = Z’ co (23) 

m mm’ ’ 


where the subscript 00 indicates that Z' has been formed with Z e open 

m 

circuited. In such an operation, the corresponding mixer input impedance 

seen by the embedding circuit therefore is: 

^in = Z m = Z ' oo . (24) 

m mm 

In particular, the IF output impedance is given by: 


Zip = Z± n = Zq = Z' c 
out 0 00’ 


(25) 


2.2.4 Conversion Loss 

The conversion loss from sideband j to sideband i in a mixer is 

A power available from the signal source at sideband w-: 

L . = ( j . 

converted power from the signal source, delivered to 

the load at sideband 

Consider for the moment only the intrinsic SIS junction shown in Fig. 3. 
The power available from impedance (z e // Z c ) at the j-th sideband is 

j j 


^available = 1/4 I SI’ | 2 Re[z e // Z c ] . (26) 

j 1 j 

The power delivered to a load impedance (z e ^ // at sideband i is, 

using (21) : 
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^delivered 


! z ij | 2 l«Ij| 2 Re[z e> 1 1 Z C- ] 
J J 11 

I z e . // z c.l 2 
1 1 


(27) 


Dividing (26) by (27) 
tunnel ■junction: 


gives the conversion loss L' of the intrinsic 

ij 


^available 


^delivered 


|z e // Z c I 2 | Z e // z c 1 2 

1 1 3 J (28) 

L' = 

ii ' 2 

4| Z ±i | Re [z e // Z c ] Re[z e // Z c ] 
j j i i 


The conversion loss of the actual mixer, L , is equal to L' since 

ij ij 

the linear parasitic junction capacitance dissipates no power. 

2*2*5 Calculation of J 

The value of the function J is required in the calculation of the 
small signal conversion parameters. J is a complex function of a real 
argument, defined as: 

J(v) = I^K ( v ) + i Idc( v ) (6) 

where i = / -1 , Idc( v ) i s the current in the SIS junction at DC 
voltage v, and 1^ is the Kramers-Kronig transform of Ido defined 
mathematically as follows: 
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V 

[i D c(v) - - ] 

1 “ R N dV (7) 

I^Cv) = P - / 

If —00 V ~ V 

The following properties of Eq. (7) are worth noting: 

(i) The Kramers-Kronig transform is identical to the Hilbert 
transform [16]. 

(ii) Note that Eq. (7) contains a singularity at V = v. To enable 
a finite integration to be calculated, the Cauchy principal value is used 
wherein the singularity is approached equally from both directions on the 
V axis. 

(iii) The ohmic portion of the I-V characteristic is subtracted out 
of the numerator of the integrand in the above integral such that the 
integral is finite. Fig. 4 shows the area under the I-V curve is infinite 
unless the linear or ohmic portion is removed. The linear or ohmic 
portion of the DC I-V curve contributes only a constant term to the value 
of Irk* As only subtractions of Ij^ terms are relevant to the calculation 
of the value of the mixer properties, the linear portion of the DC I-V 
curve may be subtracted out without error [17]. 

(iv) The calculation of the Kramers-Kronig transform does not proceed 
directly according to the definition given above in Eq. (7), but via the 
Fourier transform as follows [16] : 



I 


5 


2 3 4 

VOLTAGE (mV) — ► 



h} SIS DC * -v Curve; (b) Same I-V Curve With Ohmic Portion 
Removed. 
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(i) Form the Fourier transform of the function [ipc - — ] ; 

R N 

(ii) For values in the positive Fourier domain, multiply by -i; for 
values in the negative Fourier domain, multiply by +i; 

(iii) Invert the resulting transform using a fast Fourier transform. 

This method represents a significant saving in CPU time, particularly 
in view of the need for fast Fourier transform routines in other parts of 
the computer programs GETCHI and SISCAP. 

2*3 Mixer Noise Theory 

The noise observed in a SIS tunnel junction comes from three sources: 
(i) shot noise due to the statistical nature of the current flow across 
the junction, (ii) thermal noise due to the random motion of the charge 
carriers in the embedding network, and (iii) quantum noise, which occurs 
by virtue of the zero point energy. 

A theory capable of predicting the contribution of the quantum noise 
in a circuit analysis computer program has yet to be derived. Indeed, the 
problem of quantum noise in tunnel junction mixers, discussed in Refs. [4] 
and [18], has not been clearly resolved, and it will be ignored here. 
General considerations [19] require that a “high gain linear amplifier", 
such as the SIS mixer, add at least a half photon of fluctuation energy, 
referred to its input, to any incoming signal. This results in a minimum 
noise temperature = hu)/2k, which is small at the frequencies of 
interest (2.8K at 115GHz). 
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The equivalent circuit of the SIS tunnel -junction, including noise 


2 2 

sources, is shown in Fig. 5. 6i and 6i are the mean-square values of 

T S 

the thermal and shot noise currents in the frequency range f to f + Af. 


These current sources can be regarded as generating a multitude of quasi- 


sinusoidal frequency components, each with its own amplitude and phase. 


In the multifrequency multiport equivalent circuit of the mixer (Fig. 3) , 


the noise sources can be included by connecting a noise current source at 


each sideband frequency to the appropriate port of the augmented network. 


2*3*1 Thermal Noise 


Thermal noise generated in any embedding resistance has components 


which are uncorrelated at the various sideband frequencies. Let Si' 

T 

m 

represent the quasi-sinusoidal component at sideband frequency o> m of the 


thermal noise current source in Fig. 5 and let 6V be the sideband noise 

T 

m 

voltage produced by 51' . The noise voltage produced at the IF port of 

T 

m 

the augmented network (Fig. 3) by the thermal noise at all the sidebands 


can be found using Eq. (21): 


5V 


T 


0 


Z' 51' 
0 T 


( 29 ) 


where 6l' = [..., Si' 

T T l 



Si’ ,...] is the vector of input 
T-l 


thermal noise currents at the sideband ports of Fig. 3 and Z' is the 

0 


center row of the augmented impedance matrix Z'. From (29): 
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Fig. 5 : SIS Noise Equivalent Circuit 
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6V • 6V* = Z' 61' • fz' 6l')t = 2' 61' 6 1 ' + Z't . (30) 

T q T q 0T OT OTTO 


where t Implies the complex conjugate transpose of the vector. Taking 
the ensemble average gives: 




~ ~ t ~ t 


<61* 61' > Z 
T T 


0 


~ ~ t 


(31) 


The square matrix <6I' 61' > is the thermal noise current correlation 
matrix. Since the thermal noise components at the various sideband 
frequencies are uncorrelated, the matrix is diagonal. Including the 
Planck quantum correction to the Rayleigh-Jeans blackbody formula, this 
matrix is: 


<6l' 

T 


61 > = 0 


m 


m' 


m * m' (32a) 


4fuo m B Re(l/Z e ) 

j- m for m = m' (32b) 

<6l' 6 1 ' > = 

T m T m [exp(hm m /kT) - l] 

2*3*2 Shot Noise 

The shot noise in a mixer arises from the current produced in the 
tunnel junction conductance by the local oscillator and DC bias. The 
formal derivation for the output shot noise follow the thermal noise 
derivation. The mean square shot noise voltage produced at the IF port 
of the mixer is given by an equation analogous to Eq. (31) but with 
subscript T replaced by subscript S. 
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The shot noise can be considered as white (Gaussian) noise, amplitude 
modulated by the LO waveform. Dragone [20] and Uhlir [21] have 
investigated the properties of this modulated noise and have shown that 
there is a partial correlation between the quasi-sinusoidal components at 
the various sideband frequencies. The correlated components at these 
sidebands are down converted in the diode to the intermediate frequency 
where they add vectorially. Therefore, the shot noise current correlation 
matrix is not diagonal. It is given by Tucker [1] for quantum mixer 
theory as follows: 

t 00 * 

< 61 ' 6I ' > = e l W L0 (nm) W L0 (n'u>) S m _ m . >n ._ n 

m V n , n ' =-°° 

{coth[B(eVQ + n'hrn + hai m i)/2] Ipc^O + n 'fWe + hio m »/e) 

+ coth[ 8(eVg + nhco - ha^O/2] IdcK v 0 + n ^“/ e “ fi^m'/e)} (33) 

where 8 = 1/kT. 

2»3*3 Total Mixer Noise 

The total output noise of the mixer is obtained by combining the 
thermal and shot noise components: 


/v /v /%. 




<|V N | 2 > = Z' [ <5 1 * 6l't> + <6x ' 6l' t >] Z ,! 
0 OSS T T 0 


(34) 


The noise power delivered to the IF load from the mixer is: 


p 0 = <|Vn 0 I 2 > *e[z eo ]/|Ze o |2 


(35) 
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The equivalent input noise temperature T^ of the mixer is the 
temperature to which the signal source conductance must be heated to give 
the same output noise from a noiseless but otherwise identical mixer as 
the actual mixer would produce when its signal source conductance was 
maintained at absolute zero temperature. Thus, 

T m A P 0 L 01 /kB . (36) 

Note that Eq. (36) does not use the quantum corrected form of Johnson 
noise. The quantum corrected form used as an alternate definition would 
give a slightly different numerical result. In either case the noise 
power is unambiguous. 

Tm, as defined above, is the single sideband noise temperature, Tggg. 
When describing the performance of a mixer whose physical input port 
is coupled to both the signal and image frequencies, it is more convenient 
to talk in terms of a double sideband noise temperature Tpgg. Tggg is the 
temperature to which the signal and image source conductances must be 
heated to give the same output noise from a noiseless but otherwise 
identical mixer as the actual mixer would produce when its signal and 
image source conductances were maintained at absolute zero temperature. 

For mixers in which the conversion loss from the signal and image 
frequencies, Lqi and Lq-j, are equal: 

t DSB = T SSB/ 2 (37) 

and in general when Lqi * Lq-i 

t SSB 

t dsb = 

Lqi (38) 

1 + 

L 0-l 
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2*4 Comment on the Effect of Various Program Control Parameters 

In transforming these procedures into a workable computer program, 
there is a practical limit on the number of harmonics of the local 
oscillator which can be computed in the nonlinear large signal analysis. 
Also, the small signal admittance matrix Y will be truncated above some 
finite sideband number. 

For the voltage update large signal analysis, this truncation is 
equivalent to terminating the intrinsic junction in a short circuit at all 
frequencies higher than the truncation frequency. For a large enough 
number of harmonics, this approximation is reasonable since the cjpRj^Cj 
product of the junction is usually greater than unity. These same remarks 
apply to the truncation of the small signal admittance conversion 
matrix Y. 

For the multiple reflection large signal analysis, the 
justification for the truncation is not so clear. In this method, the 
truncation is equivalent to terminating the intrinsic SIS junction in the 
characteristic impedance of the transmission line at higher frequencies. 
Although the parasitic capacitance clearly leads to short circuit 
terminations at these higher harmonics, it has been found in practice that 
the errors so caused by the truncation are not significant. 

It should be noted that, in spite of the high a) p R N Cj products, a 
large number of points (e.g. 64 harmonics of 113.9GHz LO) still needs to 
be considered by the program since the intrinsic response of the junction 
is still exceptionally fast (e.g. 0.5 psec rise time for 3mV gap voltage 
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device). However, embedding impedance measurements do not have to be made 
right up to the 64th harmonic, since the shorting effect of the junction 
capacitance say above the 10th harmonic makes the mount impedance external 
to the capacitance irrelevant; the junction merely sees a short circuit. 
The number of harmonics treated in the program is automatically equal to 
one half of the number of time intervals into which each LO cycle is 
divided (Sec. 3.4. v.i). 

2*5 Nonlinear SIS Device Equation Calculations 

Three problems arise in putting the nonlinear quasiparticle current- 
voltage equation, 

NL v ( c ) . ,t 

I (t) = + Im (U*(t) / xCt-t'Mt^dt'} (39) 

R N 

in a form suitable for numerical solution. These are: 

(i) How to approximate the lower limit of -» in the integration; 

(ii) How small must the increment St' be (6t' is the numerical 
approximation to dt'); and 

(iii) As the two terms on the right hand side of the above equation 
are approximately equal in magnitude and opposite in sign for all DC bias 
voltages below the energy gap voltage, relative errors in I NL are likely 
to be large in this region. 

The problem of the lower limit of - 00 is identical to the problem of 
deciding how much of x is required before its asymptotically decreasing 
contributions at higher arguments are insignificant. The solution to this 
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M 2M 

problem is obtained by numerically evaluating / x( £ ) dt and / x(t) dt. 

0 0 

M is changed until the two integrals agree with each other to within 0.1%. 
This procedure was done for both an ideal SIS junction and also a typical 
real SIS junction. The results of this investigation are as follows: 


M ideal > 95 r 

(40) 

^real ^ 16 r 

(41) 


h 

where t = — is a normalizing time constant, with h = Planck’s 
2 A 

constant and A = the energy gap. It is emphasized that the value of M rea i 

will depend on the sharpness of the DC I-V curve of the device to be 

analyzed. The curve so used to obtain the value of 16 t was relatively 

sharp as may be seen in Fig. 6. Should a user have a significantly 

sharper DC I-V curve, then by trial and error, the value of the program 

control parameter NOCHI should be varied to obtain the appropriate value. 

It should be noted at this stage that the integral to be evaluated is 
t 

/ x(t~t') U(t') dt' and not / x(t') dt' as used above. By experiment, 

0 

it has been found that selecting a lower limit based on an 0.1% error in 
the latter integral leads to no more than 0.1% error in the former. 

The second problem, the size of St, is solved in the following 
manner. The following two numerical integrals, expressed as sums, are 


computed: 
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VOLTAGE (mV)— »• 


Fig. 6: I-V Curve Used to Determine Numerical Constraint Parameters 

in Section 2.5 
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N 2N kfit St 

l x(k<$t)6t as l x( ) — 

k=l k=l 2 2 

fit is adjusted until the two sums are within 0.1% of each other. The 
value obtained for fit is the same for both real and ideal SIS DC I-V 
curves, since the rise time of the x function in each case is the same. 

The value of fit obtained is: 

fit 

— < 0.064 (42) 

T 

for both real and ideal SIS DC I-v curves. 

The third problem with the nonlinear SIS circuit equation lies in the 
addition of the two nearly equal and opposite terms. These terms, for 
voltages below the energy gap, are approximately equal in magnitude but 
opposite in sign. This is so since the sum must equal the leakage current 
in the case of zero LO power applied. Even with LO power applied, the DC 
pumped current is often of the same order of magnitude as the leakage 
current. 

The following table indicates the accuracy required in calculating 
t 

the term Q = Im{u*(t)J x( t-t ') U(t')dt'} to ensure a 1% error bound in 

—CO 

I^(t) for the SIS junction of Fig. 6: 
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VOLTAGE 

% MAX ERROR IN Q 

lmV 

0.0062 

2mV 

0.024 

3mV 

0.66 

4mV 

0.86 

5mV 

0.92 

6mV 

0.94 


TABLE 1 


As expected, the error requirements on Q are most severe in the areas 
below the energy gap if a constant percentage error is to be maintained in 
I NL (t). However, it would seem intuitively correct that we could relax 
the percentage error in I^Ct) for very low values of I NL (t) without 
appreciably affecting the calculated conversion loss or noise temperature 
of the mixer. 

A summary of the digitization constraints on x based on this 
discussion is given in Table 2, for a real and the ideal SIS I-V curves. 
For comparison, the values required to perform an analogous time-domain 
analysis of the ideal Josephson pair tunneling currents, from Gayley 
[22] , are also presented. 



IDEAL 

SIS 

REAL 

SIS 

IDEAL 

J0SEPHS0N[ 22] 

M 

95x 

16t 

80t 

6t 

0.064t 

0.064t 

0.020t 

No. of Points 
of x 

1484 

250 

4000 

• 


TABLE 2 
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It should be noted that the x function derived from the SIS real I-V curve 
decays much more rapidly than that x function corresponding to the SIS 
ideal I-V curve (16t compared to 95x). Fig. 7 illustrates this. The 
decay rates for the Josephson pair current x function and the ideal 
quasiparticle (SIS) x function may be shown analytically to be the same 
and the two independent results obtained for M in Table 2 verify this. 

The much smaller step 6t required for the Josephson pair current case 
compared to that for the SIS calculations is a result of the shorter 
response time for the Josephson currents. Typically, the Josephson re- 
sponse time is approximately three times smaller than that for the quasi- 
particle currents. As may be inferred from the "number of points" values 
in the table, a substantial saving in computer memory is obtained by 
analyzing mixers containing real SIS junctions rather than ideal SIS 
junctions. 

Having established the digitization constraints on x> those for the 
SIS I-V curve automatically follow. x(t) is a sine Fourier transform of 
the I-V curve (minus the ohmic part) , from which it follows that the 
digitization constraints for the I-V curve are: 


v TOTAL > 8 v gap 


(A3) 


<5V < 



(44) 


where: the DC I-V curve is digitized from 0 to VTOTAL volts using an 
interval of SV; 


V 


gap 


A 


2A 

= — ; and 

e 

= energy gap value in Joules. 
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TIME (psec)— *• 


Fig. 7 : Comparison Between CHI Function Derived From: (a) An Ideal 

SIS DC I-V Curve; (b) A Real SIS DC I-V Curve 
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It should be emphasized that these digitization constraints were 
obtained from one particular DC I-V curve, albeit a sharp I-V curve. 
Clearly, different I-V curves will lead to differing digitization 
constraints. However, on the basis of our experience, the above 
digitization constraints have ensured sufficient numerical accuracy with 
all the measured I-V curves we have tested. 

The above digitization constraints are embodied into the program 
GETCHI which is described in Sections 3.1 and 3.2. 
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3. Description of the Mixer Analysis Program SISGAP and Program GETCHI 

Using the theory given in the previous sections, a user oriented 
computer program was written to analyze the performance of SIS mixers. 

The programming language used was FORTRAN. The program consists of two 
parts, SISCAP and GETCHI. 

The program GETCHI must be run first, since it provides information 
on the storage requirements necessary for an accurate analysis in the 
program SISCAP. Also, GETCHI outputs the discretized I-V curve and the 
discretized CHI curve in a form suitable for input to the program SISCAP. 
GETCHI requires the following as inputs: (i) the I-V curve in some 

arbitrary yet accurate digitized form; (ii) the number of junctions 
represented in the input DC I-V curve; (iii) the LO frequency; and (iv) 
the gap voltage of the single equivalent device. Further details on the 
input/output data structures can be found in Section 3.2. 

The program SISCAP requires as inputs the following: (i) the 

embedding impedances seen by the SIS tunnel junction (excluding the 
parasitic junction capacitance contribution) at each harmonic of the local 
oscillator and at the sidebands; (ii) the SIS DC I-V curve (minus the 
ohmic part, as provided by the program GETCHI); (iii) the SIS x 
curve (which is derived from the I-V curve by the program GETCHI, to 
be described later); (iv) the various operating conditions for the tunnel 
junction mixer, e.g. DC bias voltage, LO amplitude, the frequencies of the 
LO, RF and IF; (v) the normal resistance % and junction capacitance Cj of 
the SIS tunnel junction, and finally (vi) the parameters used to control 
the various numerical calculations within the program. 

The output includes (i) the conversion loss, (ii) the RF Input 
impedance; (iii) the IF output impedance, and (iv) the single sideband 
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mixer noise temperature. 

The remainder of this section explains both programs in detail and 
illustrates the steps used for running them. A complete documented 
listing of the two programs appears in Appendix A and a general flowchart 
is given in Fig. 8. 

3. 1 GETCHI Program Implementation 

The main routine of program GETCHI calls various subroutines 
described below in the sequence indicated. 

Subroutine EVEN provides an evenly-spaced (along the voltage axis) 
internal version of the input I-V data which is used as the reference data 
for subsequent processing by the program. Linear interpolation is used to 
calculate these points. If the input data is for an array, subroutine 
EVEN reduces this to its single junction equivalent. Subroutine EVEN 
requests as input the number of tunnel junctions arrayed. If not an 
array, unity should be entered. The given I-V data must be input on 
device 7. 

Subroutine REFINE then takes the I-V internal working data produced 
by subroutine EVEN and extracts the ohmic or linear portion of the I-V 
curve . 

Subroutine INEXP then interpolates the internal working I-V data 
stored within the program such that the new spacing and number of points 
correspond to the correct spacing required for further processing prior to 
the conversion to CHI data. The LO frequency (FREQ) and the single 
junction gap voltage (VGAP) are required as input so the program can 
provide lower limits for the number of LO cycles (NOCYC) and number of 
points per cycle (NOPNTS). These limits ensure that an accurate 
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Fig . 8(a): Flowchart of Program GETCHI 
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Fig. 8(b) ; Flowchart of Program SISCAP 


OUTPUT SMALL SIGNAL LARGE SIGNAL INPUT 

DATA AND NOISE ANALYSIS ANALYSIS DATA 
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calculation Is performed in the program SISCAP. 

Following this, subroutine REFINE is again called to ensure that the 
last point of the I-V data is zero current. 

Subroutine IVTOCH uses a fast Fourier transform and converts the 
working I-V data to the equivalent CHI data file for output on device unit 
number 8. 

Finally, subroutine REEXP outputs suitably digitized I-V data for use 
in evaluating the current Ipc and the complex valued parameter J in the 
SISCAP program. To output this I-V data, the program requests both the 
number of points required (NOIV) and the spacing between the points 
(DELTAV). Suggested values for these parameters are NOIV = 1024 and 
DELTAV = 0.02 mV for a 3mV gap SIS device. For other gap voltage values, 
DELTAV should be linearly scaled to suit. The DC I-V data is output on 
device unit number 9. 

Output to the console from program GETCHI is a summary of the values 
of key parameters to be input to program SISCAP. The parameters are: 
NOCHI, NOIV, NOPNTS and RN. The first three parameters are also used to 
determine the memory allocation requirement in SISCAP as per Appendix 

A. 5. 

3.2 Running the GETCHI Program 

A listing of the program GETCHI appears in Appendix A along with the 
output of a run. The comments in the listing provide a description of the 
FORTRAN coding. 

Prior to running program SISCAP, GETCHI must be run to provide 
information on memory storage allocation for SISCAP, values of key 
parameters in program SISCAP as well as to provide the CHI and* the DC I-V 
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curve, each with the appropriate discretization for input to program 
SISCAP. GETCHI requires as input to it, the DC I-V curve of the SIS 
junction. Multiple junction I-V curves may be input to GETCHI but the two 
output files from GETCHI containing CHI and the DC I-V data will reflect 
those values appropriate to the equivalent single junction device. The 
program will request the number of junctions (NOJUN) in the series array; 
a value of unity should be input if a single junction is used. 

The number and distribution of the points given in the input DC I-V 
data file to program GETCHI are at the discretion of the user but must 
conform to the following guidelines: 

(i) the first I-V data point must be the origin (0.0, 0.0); 

(ii) the input I-V points need not be evenly spaced; 

(iii) the span of the input points must be such as to extend well 
into the linear portion of the I-V curve, i.e. the points must adequately 
cover the entire nonlinear portion (typically to a voltage four times the 
gap voltage); and 

(iv) the I-V points must represent a monotonic function. 

In order to obtain the CHI function, it is necessary to input the 
frequency of the local oscillator (FREQ), the gap voltage of the single 
junction equivalent SIS device (VGAP), the number of points per LO cycle 
(NOPNTS) to be used in the calculation and the number of cycles of points 
required to accurately describe CHI (NOCYC) in accordance with the 
conclusions of section 2.5. The program offers assistance in deciding 
what the minimum acceptable values for NOPNTS and NOCYC are. Both these 
numbers must be powers of 2 . The parameter, NOCHI, whose value must be 
known for input to the program SISCAP, is equal to N0PNTS*N0CYC. As such, 
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the value of NOCHI is always greater than or equal to NOPNTS. The program 
also outputs the value of RN to be input to the program SISCAP. The 
program requests both the number of points for output (NOIV) to the file 
for program SISCAP and the voltage spacing between these points (DELTAV). 
The number of points (NOIV) must be a power of 2 and 1024 has been found 
to give satisfactory results. DELTAV, the voltage spacing between the 
points, is typically 0.02 mV for a 3mV gap device. The DC voltage span 
covered by these typical parameters is 20.48 mV (= 1024 x 0.02mV) which 
clearly covers the nonlinear portion of the I-V curve. 

3.3 SISCAP Program Implementation 

The program commences in the MAIN routine and inputs all the 
parameters necessary to control the operation of the program. Next, the 
values of CHI (x) and the SIS DC I-V characteristic are read in. At this 
point, control is transferred to the subroutine GETJ in order to obtain 
the Kramers-Kronig transform of the DC I-V curve and hence the complex 
function J. Control is then transferred back to the MAIN routine to input 
the embedding impedances at the local oscillator harmonics (subroutine 
GETZ) and sideband frequencies (subroutine GETSSB) . 

All information has now been input and the program is now ready to 
perform the large signal analysis. There are three methods which may be 
used: 

(i) the voltage update method for lOpR^Cj products in excess of 0.5; 

(ii) the multiple reflection technique for other oipR^Cj products; and 

(iii) the three frequency approximate analysis. 

The voltage update method proceeds as follows. An estimate of the 
junction voltage is made. A call is made to subroutine SISNL which 
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calculates the current flowing into the SIS device, given the input 
voltage. Next, subroutine LINEAR calculates the voltage required across 
the embedding network terminals consistent with the current calculated in 
subroutine SISNL. Finally, COMPAR performs a comparison between the new 
voltage waveform and the initial estimate of the voltage waveform to 
determine if a solution has been obtained. 

Alternatively, the multiple reflection technique may be used. This 
proceeds by a call to subroutine MRT which performs the overall iterative 
algorithm. It in turn calls subroutine SISCHA which sets up various 
vectors such that the succeeding call to subroutine SIS will solve the 
nonlinear SIS tunnel junction equations. 

If required, a three frequency model calculation can be performed 
instead of the more general nonlinear analysis routines described above. 

Following the completion of the large signal analysis calculations by 
whatever large signal analysis method, a call to TOPS initiates 
calculations to obtain the small signal and noise results. TOPS is the 
principal subprogram controlling the signal and noise analysis programs. 

It initializes the working matrices to zero and calls SIGNL and SIGLIN. 
SIGNL calculates the matrix conversion elements which are computed from 
the SIS device parameters together with the large signal analysis 
waveforms. The embedding impedances at all the relevant sidebands are 
calculated by subroutine SIGLIN and, by virtue of the linearity, the 
matrix so calculated is diagonal. The combination of the subroutines 
SIGLIN and SIGNL produce the augmented admittance matrix. The combination 
of the two subroutines INVRT and LOSCAL invert the augmented admittance 
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matrix and perform the arithmetic to calculate the conversion loss and the 
output impedance. The noise analysis follows next using calls to SISNOI 
and RESNOI. SISNOI computes the contribution from the shot noise of the 
SIS device while the thermal or Johnson noise is computed from subroutine 
RESNOI. Subroutine PROCES then takes the complete noise matrix and, using 
the previously calculated impedance conversion matrix, computes the 
equivalent input noise temperature of the mixer. 

If a three frequency analysis is invoked, the small signal and noise 
analyses are calculated as follows. Subroutine THEORY provides the 3x3 
matrix elements describing the small signal conversion and noise 
correlation process. Subsequent calls to subroutines SIGLIN, INVRT and 
LOSCAL calculate the conversion loss of the mixer. The single-sideband 
mixer noise temperature is calculated by calls to subroutines RESNOI and 
PROCES. 

For a more detailed description of the operation of the program, the 
reader is referred to the comments in the program listing and the general 
flow chart of Fig. 8. 

3*4 Running the Mixer Analysis Program SISCAP 

A listing of the SIS mixer analysis program SISCAP appears in 
Appendix A along with the output of a run. Using the IBM Fortran 77 
compiler, the execution time for this particular listing is 0.5 mins on an 
AMDAHL V6 computer. The comments in the listing provide a description of 
the FORTRAN coding. To run the program, the following information must be 
supplied as files accessed via the following FORTRAN unit numbers: 

(i) UNIT NO. 15: The embedding network impedances at the LO 

« 

frequency and the higher harmonics. Preceding the impedance values, there 
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must be a number indicating the number of impedances in the file. It is 
read in under 15 format. All impedance values are to be in ohms and are 
read in as complex variables using format 2E15.8. The program expects 
there to be NOPNTS/2 impedances in this file. If there are less than 
this, the program assumes that those impedances given are zero. It is 
important to note that these impedance values should not Include the 
effect of the parasitic junction capacitance . 

(ii) UNIT NO. 14: The sideband impedances in complex form, again in 
ohms. The ordering of information In this file is as follows: number of 
impedances in the file, IF impedance, fundamental USB impedance (z e (uji)), 
fundamental LSB impedance (z e (oj_i)), second LO harmonic USB impedance 
(z e (a> 2 )) and the second LO harmonic LSB impedance (z e (u)_ 2 )], etc. Lower 
sideband impedances should not be input in their complex conjugate form , 
the program does this operation. As for the previou? case, a number 
indicating the number of impedances in the file must be submitted as the 
first entry in the data file, using the format 15. The complex impedances 
are read in under 2E15.8 format. The program expects there to be (NOPNTS 
- 1) impedances in the files. If there are less than this, the program 
assumes that those impedances not given are zero. It is important to note 
that these impedance values should not Include the effect of the parasitic 
junction capacitance . 

(iii) UNIT NO. 8: The DC SIS I-V curve Is input in the form of the 

number of points (format 15) in the file followed by the digitized I-V 
curve (format 2E15.8). The ohmic part, will have been subtracted out 
by program GETCHI (See Section 3.2) which should be used to provide the 
DC I-V curve data. 

(iv) UNIT NO. 4: The CHI values for the SIS tunnel junction device 
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with the following format: the number of points in the file (15), 

followed by the values of CHI (2E15.8). The points for this file can be 
generated by using the program GETCHI (See Section 3.2). 

(v) UNIT NO. 11: This file contains all the SIS device parameters 
together with the various program option values which affect the operating 
and performance characteristics of program. The parameters and their 
function are listed sequentially below (all integers are input in the 
format 15, all real values are input in the format El 5. 8): 

(a) FREQ: frequency of the local oscillator (LO) in Hertz; 

(b) FIF: intermediate frequency (IF) in Hertz; 

(c) CAP: parasitic junction capacitance of the SIS device in Farads; 

(d) RN: normal resistance of the SIS device in ohms. This value is 

supplied by the program GETCHI; 

(e) TEMP: operating temperature of the junction in Kelvins; 

(f) IBASE2: binary logarithm of NOIV; 

(g) NOIV: number of points in the I-V data file (this must be a 

power of 2). The value of NOIV must be '■he value used in the program 
GETCHI. Assistance in selecting the value of NOIV is given when program 
GETCHI is run (See Section 3.2); 

(h) IBASE1 : binary logarithm of NOPNTS; 

(i) NOPNTS: number of points into which each LO cycle is subdivided 

(this must be a power of 2). This value must be obtained from the program 
GETCHI. 

(j) MMAX: number of sidebands to be considered (e.g. 9 implies 4 

upper, 4 lower and the IF sidebands are to be included in the analysis). 
MMAX must be an odd number. Typically, a value of 15 has been* found to 
give satisfactory results. Its value is always limited to a maximum of 
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N0PNTS/2 - 1; 

(k) NMAX: number of photon points to be considered in the small 

signal analysis. This establishes the number of terms summed in Eq. 16 
to calculate the admittance matrix elements. At typical LO drive levels, 
at least 6 photon points on each side of the DC bias point are required 
to ensure reasonable accuracy. Therefore, the minimum for NMAX is 

2*6 + 1 = 13. The maximum value for NMAX is NOPNTS/2 - 1; 

(l) ISIG: indicates the type of mixer analysis to be performed: 

1 = three port model, 2 = general mixer analysis using voltage update 
method, 3 = general mixer analysis using multiple reflection technique; 

(m) NOPER: maximum number of iterations allowed in the large signal 

general mixer analysis. This parameter acts as an emergency brake in case 
the nonlinear analysis does not converge. Typically, its value should be 
100. The value of this parameter is irrelevant if the three port 
approximation is used (ISIG = 1); 

(n) NOCHI: number of points required in the CHI (x) file. The value 

of this parameter is obtained from the program GETCHI; 

(o) NODAT: number of different DC bias voltage cases to be analyzed 

(see also DELV and VSTART); 

(p) TOL: the maximum relative error allowed between successive 

iterations in the large signal analysis before the solution has 
converged. Since the computational accuracy of each computer varies from 
machine to machine, the value assigned to TOL must therefore reflect the 
respective rounding error of each individual machine. Typically, TOL 
equal to 0.00001 has been found to be effective on IBM 32 bit machines; 

(q) EPS: in a vector of values, EPS is the fraction of the maximum 

component below which all values are considered to be insignificant and 
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therefore unduly affected by rounding error and the like. Since the 
computational accuracy of each computer varies from machine to machine, 
the value assigned to EPS must therefore reflect the respective rounding 
error of each individual machine. Typically, EPS equal to 0.00001 has 
been found to be effective on IBM 32 bit machines; 

(r) DELV: increment in DC bias voltage between the different bias 

cases considered for analysis (in volts). See also NODAT and VSTART; 

(s) VSTART: first of the NODAT values of DC bias voltage (in volts). 

See also NODAT and DELV; 

(t) VLOSTR: first of the local oscillator peak voltage values (in 

volts). See also DELVLO and NLO; 

(u) DELVLO: increment in local oscillator peak voltage value (in 

volts). See also VLOSTR and NLO; 

(v) NLO: number of different LO voltage cases to be analyzed. See 

also VLOSTR and NLO; 

(w) P: convergence parameter used in the voltage update large signal 

analysis; its value lies in the range 0 < P < 1. Typically, its value is 
unity; 

(x) RID: value of the identity resistance (in ohms) to be used in 

the voltage update large signal analysis method. Typically, R^p set equal 
to Rfj/2 has been found to give satisfactory results; 

(y) Z0: value of the characteristic impedance (in ohms) of the 

transmission line used in the multiple reflection large signal analysis 
method. Typically, Zq set equal to Rjj has been found to give satisfactory 
results. 

Output is received from program SISCAP as follows: 

(a) UNIT NO. 6: the console output. 
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3.5 Verification Tests 

The program has been verified for the case of large junction 
capacitance by comparison with the predictions of the three-frequency 
analysis. In a number of cases (not restricted to a large junction 
capacitance) , the large signal waveforms predicted by this program were 
found to be consistent with those calculated by an independently developed 
SIS nonlinear analysis program which was used to study chaos [23] . The 
mixer temperatures predicted by this program have not been independently 
verified but do appear reasonable. This project was terminated before the 
program was extensively employed, and it is possible that errors may occur 
in various untested parameter ranges. 
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4. Example 

In this section, the mixer analysis program is used to study the 
effects of parasitic junction capacitance on the performance of an SIS 
mixer. The experimentally measured I-V curve, shown in Fig. 9, of an SIS 
element [24] for which % = 72 ft was assumed to be in a mixer circuit with 
a 113.9 GHz local oscillator. The intermediate frequency was 1.0 MHz and 
the IF load resistance was 50 £2. In this low IF limit the signal, image, 
and LO source impedances must all be equal, and they were chosen to be 55 
+ i92 £2 (which includes Cj), the value which maximizes the conversion 
efficiency at a dc voltage in the center of the first photon step in the 
three-frequency approximation. This SIS mixer was analyzed for a wide 
range of junction capacitance. In each case the termination impedance at 
all higher harmonics and sidebands ( not including the LO and its side- 
bands) was the parallel combination of the capacitive reactance and a 
resistance of 72 £1 (arbitrarily chosen equal to %) • 

For each value of capacitance the LO power and the dc voltage were 
optimized for maximum conversion on the first photon step. Typical 
junction current and voltage waveforms resulting from the numerical 
computation in Section A. 3 are shown in Fig. 10. The overall results for 
varying junction capacitances are not plotted, for fear that relatively 
large interpolation errors may have made them unreliable. The trend of 
the results is however quite clear. 

The largest value of mixer conversion gain, a few dB, was found in 
the limit of large capacitance, and agreed with the three-frequency 
approximation. The conversion gain decreased to a minimum of about unity 
at an ojpR^Cj product of unity, but then increased again for smaller u)p%Cj 
values. The mixer noise temperatures corresponding to these results 
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showed an even less dramatic variation with capacitance, remaining between 
20 and 25 K. The mixer output impedance was negative at all capacitance 
values, while the input impedance was always positive. It must be 
emphasized that even if these conclusions are verified for this specific 
mixer circuit, they may be far from typical and should not be taken as 
general. 




JUNCTION CURRENT (/iA)-» JUNCTION VOLTAGE (mV) 



Fig. 10: Typical Junction Waveforms: (a) Voltage; (b) Current 
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5. Conclusions 

A FORTRAN computer program has been developed for analyzing SIS 
mixers with arbitrary embedding impedances at all LO harmonics and 
sidebands. This program has been verified using the three-frequency 
approximation and the multiple reflection algorithm. Typical current and 
voltage waveforms at the SIS -junction are plotted. Although a thorough 
study of the effect of varying junction capacitance has yet to be 
completed, in one specific example optimum mixer performance occurs for 
large or very small values of ojpRjjCj. Low values of this parameter, 
roughly 0.3 to 2.0, lead to a deterioration in mixer performance. 

Tests have also been performed on the relative speeds of the two 
nonlinear large signal analysis methods. For a typical SIS mixer 
situation, the voltage update method provides an order of magnitude 
improvement in total CPU time over that required by the multiple 
reflection technique. 
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APPENDIX A Supplementary Computer Program Information 
A. 1 Typical Run of Program SISJUN 

The program SISJUN generates, without input, a typical 
single-junction DC I-V curve in digitized form. It is output on unit 7 
and consists of 2049 lines of data. The first few lines of output are 
shown below. This data is to be used as input to GETCHI and is therefore 
given as a means for the verification of the correct functioning of 
programs GETCHI and SISCAP. Although SISJUN generates equally spaced data 
for input to program GETCHI, this is not a necessity for GETCHI input, as 
discussed in Section 3.2. 

A. 2 Typical Run of Program GETCHI 

Given in this section is the output of a typical run of the program 
GETCHI. It should be remembered that three files have to be defined 
corresponding to FORTRAN unit 7 (input I-V data), unit 8 (output CHI data) 
and unit 9 (output I-V data). In the example shown here, the program uses 
as input on unit 7 the output of the program SISJUN. 

The listing of the console output (given below) is self explanatory. 
It should be emphasized that, as a precaution, the two output files of CHI 
data and DC I-V data should be plotted to ensure there are no 
discontinuities in this data. Such a situation is possible when the input 
experimental data has been affected by noise during the measurement 
process. Such noise can destroy the monotonicity of the output DC I-V 
curve, a property relied upon by portions of the computer programs 


described in this report 
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FILE* SISJUN DOC A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


s I s Jun 

FI 07 DISK SISJUN DAT A1 { PERM RECFM FA LRECL 80 
EXEC FORTXCLG SISJUN XREF GOSTMT MAP 

FORTRAN H EXTENDED (ENHANCED) 

SOURCE STATEMENTS - 43, PROGRAM SIZE = 147G, SUBPROGRAM NAME - MAIN 

*** Files ’ SISJUN TEXT ’ and ’ SISJUN LISTING * have been created 
*** Loading ... file * LOAD MAP A * Is being created 
EXECUTION BEGINS. . . 

R; T=4.75/5.54 12:15:38 


FILE: SISJUN DATA A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


201,8 

0.0 0.0 
0 4S999999E-05 0. 1 6G66675E-08 
0. 9 999996 1 E-J0T5 0. 33333336E-08 
0. 14 99999SE-04 0 . 50000004 £-08 
0 19995992E-04 0 . 6 6 G 6£ 6 72E-08 
0.Z 5000001 E -04 0.82333340E-08 
0 2 999999GE -04 0 . U 000004 E-0 7 
0.34 9 99 990E-04 0. 1 1 G66C G8E-07 
0 . 3 9999 984E -04 0. 1 3333:34E-07 
0.44999 9 79E-04 0. 1 4 9999 5CE-07 
0. 4S99S973E-04 0 . 1 E 66 6 £ C 4E-07 
0.54999997 E -04 0 1 8 33 3 3 3 6 E -07 
0.59999991 E-04 0. 2<.h00< f LE-07 
0. £ 4999986E-04 0. 21 66666CE-07 
0. 69999 980E-04 0.2C333235C-,. 7 
0. 74999974E-04 0 2C 000! 0 2 E-07 
0.793S9998E-JC4 0.2lGG6£C<" 2-07 
0.84S55S92E-04 0. 2£ 3339 / 3C-07 
0.8 9 995 98 7E -04 0. 29 S5S995E-07 
0.94995581E-04 0 316C6E73E-07 
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FILE: GETCHI DOC A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


getch 1 

FI 07 DISK SISJUN DAT A1 { PERM RECFM FA LRECL 80 
FI 08 DISK CHI DAT A1 < PERM RECFM FA LRECL 80 
FI 09 DISK IV DAT A1 ( PERM RECFM FA LRECL 80 
EXEC FORTXCLG GETCHI XREF GOSTMT MAP 

FORTRAN H EXTENDED (ENHANCED) 


SOURCE 

STATEMENTS 

= 

31 , 

PROGRAM 

SIZE 

= 

67042, 

SUBPROGRAM 

NAME 

s 

MAIN 

SOURCE 

STATEMENTS 

= 

63, 

PROGRAM 

SIZE 

= 

2004, 

SUBPROGRAM 

NAME 

= 

EVEN 

SOURCE 

STATEMENTS 

= 

10, 

PROGRAM 

SIZE 

= 

494 , 

SUBPROGRAM 

NAME 

= 

REFINE 

SOURCE 

STATEMENTS 

= 

77, 

PROGRAM 

SIZE 

= 

2542, 

SUBPROGRAM 

NAME 

= 

INEXP 

SOURCE 

STATEMENTS 

= 

56, 

PROGRAM 

SIZE 

= 

1578, 

SUBPROGRAM 

NAME 

= 

I VTOCH 

SOURCE 

STATEMENTS 

= 

71, 

PROGRAM 

SIZE 

= 

2210, 

SL BPROGRAM 

NAME 

= 

FFT 

SOURCE 

STATEMENTS 

= 

15, 

PROGRAM 

SIZE 

= 

404, 

SUBPROGRAM 

NAME 

= 

JUXTA 

SOURCE 

STATEMENTS 

= 

15, 

PROGRAM 

SIZE 

= 

484 , 

SUBPROGRAM 

NAME 

= 

FACNO 

SOURCE 

STATEMENTS 

= 

63, 

PROGRAM 

SIZE 

= 

1978, 

SUBPROGRAM 

NAME 

= 

REEXP 


*** Files ’ GETCHI TEXT ’ and ’ GETCHI LISTING ’ have been created 
*** Loading ... file ’ LOAD MAP A ’ Is being created 
EXECUTION BEGINS. . . 

THIS PROGRAM, GETCHI, PERFORMS A PREPROCESSING 
FUNCTION ON THE INPUT DC I-V CURVE DATA PRIOR TO USING 
THE PROGRAM SISCAP. OUTPUT FROM GETCHI IS THE PRE- 
PROCESSED DC I-V DATA, CHI IN A!, APPROPRIATELY DIS- 
CRETIZED FORM AND VALUES FOR THE PARAMETERS: NOCHI, 

NOI V , NOP NTS AND RN 

INPUT THE GAP VOLTAGE (VGAP) IN VOLTS OF THE SINGLE 
JUNCTION EQUIVALENT DEVICE (FORMAT E15.8) 

0.003 

INPUTTING THE INPUT I-V CURVE DATA FROM UNIT 7 

INPUT THE NUMBER OF JUNCTIONS (NOJUN) IN DEVICE (15) 

1 

INPUT THE LO FREQUENCY IN HERTZ (FREQ) (E15.8) 

113. 9e9 

NUMBER OF LO CYCLES OF CHI IN SISCAP SHOULD BE > 2 

THE NUMBER OF POINTS PER LO CYCLE IN SISCAP SHOULD BE > 99 

INPUT NUMBER OF LO CYCLES (NOCYC) (15) 

INPUT NO OF POINTS PER LO CYCLE (NOPNTS) (15) 

4 

128 


FILE: GETCHI DOC A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


NOCHI = 512 

INPUT DELTAV (TYP 0.2000E-04 VOLTS) (E15.8) 

INPUT NUMBER OF POINTS REQUIRED (NOIV) ( PVR OF 2) (TYP 1024) (15) 
0.00002 
1024 


INFORMATION FOR TRANSFERRAL TO PROGRAM SISCAP: 
OUTPUTTING THE CHI FILE ON UNIT 8: NOCHI = 512 

OUTPUTTING THE IV FILE ON UNIT 9: NOIV = 1024 
NOPNTS = 128 

RN = 0. 72131729E+02 OHMS 

R; T=13. 93/15. 93 12:17-46 



58.03 


FILE: CHI 


DATA A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


512 

0.0 

0.68591 129E-13 
0. 13718226E-12 
0. 20577339E- 1 2 
0. 274364S2E- 1 2 
0. 34295564E- 1 2 
0.41154677E-12 
0.4801 3790E- 1 2 
0. 54 872 903E- 1 2 
0.61 73201 6E- 1 2 
0.68591 129E-12 
0. 75450242E- 1 2 
0.82309355E-12 
0.89168468E-12 
0. 96027526E- 1 2 
0. 1028E662E- 1 1 
0. 10974581E-1 1 
0. 1 1660490E-1 1 
0. 1 2346400E- 1 1 
0 1 30323 10E -11 


0.0 

-0.230 1 9856E+08 
-0.413349 1 2E+08 
-0.5221441 6E+08 
-0. 55727536E+08 
-0.54098384E+08 
-0.501 1 6304E+08 
-0. 45653440E+08 
-0.41 148400E+08 
-0.36080416E+08 
-0.2S899776E+08 
-0.22654000E+08 
-0. 1 4 9854 96 E+08 
-0.76890440E+07 
-0. 1 2803370E+07 
0 . 409808 90E +07 
0. 84673520E+07 
0. 1 1769160E+08 
0. 1381 2084E+08 
0. 1 4 4 265 68E + 08 


FILE: IV 

DATA 

A 

NASA/GISS 

1024 




0.0 

0.0 




0. 20000007E-04-0. 27060435E-06 
0. 4000001 4E-04-0. 541 2076 IE-06 
0. 60000020E-04-0. 8118 1088E-06 
0 . 8000002 7E -04 -0 . 10824 1 33E-05 
0. 1 0000003 E -03-0. 1 35301 7 1 E-05 
0. 12000004E-03-0. 1 6236208E-05 
0. 14000005E-03-0. 18942246E-05 
0 . 16000005E-03-0.21648284E-05 
0 . 18000006E-03-0.24354322E-05 
0. 20000 007 E -03-0. 2 70603 60E -05 
0 . 2 2000007 E -03-0 . 29 7 6 63 9 7E -05 
0. 24000008E-03-0. 3 24 724 4 4 E-05 
0. 259 999 8 8E -03-0. 351 7 84 4 6 E-05 
0.27999980E-03-0. 37884465E-05 
0.29999972E-03-0.40590494E -05 
0.3 1 999 96 4 E -03-0. 4329651 3 E-05 
0. 3 3 9999 56 E -03-0. 46002542E-05 
0.35999948E-03-0. 48708571 E-05 
0.37999940E-03-0. 51414 600E-05 


CONVERSATIONAL MONITOR SYSTEM SP302 
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FILE: SISCAP DOC A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


s I scap 

FI 11 DISK PARAM DAT A1 ( PERM RECFM FA LRECL 80 
FI 04 DISK CHI DAT A1 ( PERM RECFM FA LRECL 80 
FI 08 DISK IV DAT A1 ( PERM RECFM FA LRECL 80 
FI 14 DISK ZSSB DAT A1 ( PERM RECFM FA LRECL 80 
FI 15 DISK ZLO DAT A1 ( PERM RECFM FA LRECL 80 
EXEC FORTXCLG SISCAP GOSTMT XREF MAP 

FORTRAN H EXTENDED (ENHANCED) 


SOURCE 

STATEMENTS 

= 

170, 

PROGRAM 

SIZE 

= 

42222, SUBPROGRAM 

NAME 

= MAIN 

SOURCE 

STATEMENTS 

= 

47, 

PROGRAM 

SIZE 

= 

1810, SUBPROGRAM 

NAME 

=GETSSB 

SOURCE 

STATEMENTS 

= 

39, 

PROGRAM 

SIZE 

= 

1556, SUBPROGRAM 

NAME 

= GETZ 

SOURCE 

STATEMENTS 

= 

25, 

PROGRAM 

SIZE 

= 

2478, SUBPROGRAM 

NAME 

=LINEAR 

SOURCE 

STATEMENTS 

= 

68, 

PROGRAM 

SIZE 

= 

20228, SUBPROGRAM 

NAME 

= SISNL 

SOURCE 

STATEMENTS 

= 

27, 

PROGRAM 

SIZE 

= 

734, SUBPROGRAM 

NAME 

=COMPAR 

SOURCE 

STATEMENTS 

3 

57, 

PROGRAM 

SIZE 

= 

35034, SUBPROGRAM 

NAME 

= TOPS 

SOURCE 

STATEMENTS 

= 

26, 

PROGRAM 

SIZE 

= 

1054, SUBPROGRAM 

NAME 

“RESNOI 

SOURCE 

STATEMENTS 

= 

36, 

PROGRAM 

SIZE 

= 

1754, SUBPROGRAM 

NAME 

=PROCES 

SOURCE 

STATEMENTS 

= 

79, 

PROGRAM 

SIZE 

= 

4476, SUBPROGRAM 

NAME 

= SIGNL 

SOURCE 

STATEMENTS 

= 

31, 

PROGRAM 

SIZE 

= 

17608, SUBPROGRAM 

NAME 

“ GETO 

SOURCE 

STATEMENTS 

= 

16, 

PROGRAM 

SIZE 

- 

762, SUBPROGRAM 

NAME 

“SIGLIN 

SOURCE 

STATFMENTS 

= 

45, 

PROGRAM 

SIZE 

= 

2430, SUBPROGRAM 

NAME 

“LOSCAL 

SOURCE 

STATEMENTS 

= 

28, 

PROGRAM 

SIZE 

= 

866, SUBPROGRAM 

NAME 

= EVALJ 

SOURCE 

STATEMENTS 

= 

40, 

PROGRAM 

SIZE 

= 

1818, SUBPROGRAM 

NAME 

=S I SNO I 

SOURCE 

STATEMENTS 

= 

23, 

PROGRAM 

SIZE 

= 

998, SUBPROGRAM 

NAME 

= INVRT 

SOURCE 

STATEMENTS 

= 

25, 

PROGRAM 

SIZE 

3 

1238, SUBPROGRAM 

NAME 

“FACTOR 

SOURCE 

STATEMENTS 

= 

25, 

PROGRAM 

SIZE 

= 

1418, SUBPROGRAM 

NAME 

= SUBST 

SOURCE 

STATEMENTS 

= 

71, 

PROGRAM 

SIZE 

= 

2210, SUBPROGRAM 

NAME 

= F FT 

SOURCE 

STATEMENTS 

= 

15, 

PROGRAM 

SIZE 

= 

404, SUBPROGRAM 

NAME 

= JUXTA 

SOURCE 

STATEMENTS 

= 

15, 

PROGRAM 

SIZE 

= 

484, SUBPROGRAM 

NAME 

= FACNO 

SOURCE 

STATEMENTS 

= 

15, 

PROGRAM 

SIZE 

= 

576, SUBPROGRAM 

NAME 

= ZERO 

SOURCE 

STATEMENTS 

= 

140, 

PROGRAM 

SIZE 

= 

4574, SUBPROGRAM 

NAME 

“THEORY 
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FILE: SISCAP DOC A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


SOURCE 

STATEMENTS 

= 

26, 

PROGRAM 

SIZE 

= 

962, 

SUBPROGRAM 

NAME = DERI J 

SOURCE 

STATEMENTS 

= 

66, 

PROGRAM 

SIZE 

= 

1600, 

SUBPROGRAM 

NAME = BESO 

SOURCE 

STATEMENTS 

= 

12, 

PROGRAM 

SIZE 

= 

396, 

SUBPROGRAM 

NAME = COTH 

SOURCE 

STATEMENTS 

= 

56, 

PROGRAM 

SIZE 

= 

3042, 

SUBPROGRAM 

NAME =GEN3PT 

SOURCE 

STATEMENTS 


65, 

PROGRAM 

SIZE 

= 

1754, 

SUBPROGRAM 

NAME = TERP5 

SOURCE 

STATEMENTS 

= 

65 , 

'UROGRAM 

SIZE 

= 

6928, 

SUBPROGRAM 

NAME = MRT 

SOURCE 

STATEMENTS 

= 

20, 

PROGRAM 

SIZE 

= 

12252, 

SUBPROGRAM 

NAME =S ISCHA 

SOURCE 

STATEMENTS 

= 

1 10, 

PROGRAM 

SIZE 

= 

4450, 

SUBPROGRAM 

NAME = SIS 

SOURCE 

STATEMENTS 

= 

12, 

PROGRAM 

SIZE 

= 

430, 

SUBPROGRAM 

NAME =YFUNCT 

SOURCE 

STATEMENTS 

= 

38, 

PROGRAM 

SIZE 

= 

854, 

SUBPROGRAM 

NAME = COMP 2 

SOURCE 

STATEMENTS 

= 

8, 

PROGRAM 

SIZE 

= 

388, 

SUBPROGRAM 

NAME =F IXPHA 


*** Files ’ SISCAP TEXT ’ and ’ SISCAP LISTING ’ have been created 
*** Loading ... file ’ LOAD MAP A ’ Is being created 
EXECUTION BEGINS... 

FREQ = 0. 1 1390000E+12HZ FIF = 0. 10000000E+07HZ 

CAP = 0. 19371999E-13F RN = 0. 72 1 3 1 7 2 9E+02OHMS 

TEMP = 0.26999998E+01K 

IBASE2 = 10 NO IV = 1024 

I BASE 1 = 7 NOP NTS = 128 

MMAX = 31 NMAX = 31 

IS IG = 2 NOPER = 50 

NOCHI = 512 NODAT = 1 

TOL = 0. 10000000E-03 EPS = 0. 10000000E-03 

DELV = 0. 25000001 E-04VOLTS VSTART = 0. 27999999E-02VOLTS 

VLOSTR = 0. 35000001 E-02VOLTS DELVLO = 0. 10000000E-03VOLTS 

NLO =1 P = 0. 10000000E+01 

RID = 0 . 50000000E+02OHMS Z 0 = 0. 50000000E+02OHMS 

INPUTTING THE CHI FUNCTION 
DELTAT = 0 . 68597905E- 1 3 SECS 

NUMBER OF POINTS READ = 512 

INPUTTING THE SIS IV CHARACTERISTIC 
DELTAV = 0. 1 9997358E-04 VOLTS 

NUMBER OF POINTS READ = 1024 

CALCULATING THE VECTOR J 

INPUTTING THE LO HARMONIC IMPEDANCES 

THE FIRST 10 LO HARM IMPEDANCES (EXC CAPAC) ARE: 

0 0.0 +U 0.0 OHMS 

1 0 . 55000000E+02+U 0. 92000000E+02OHMS 

2 0.721 30005E+02+ J 0. 99999994E-10OHMS 
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FILE: SISCAP DOC 


A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


3 0.721 30005E+02+O 

4 0. 72130005E+02+O 

5 0.721 30005 E +02+0 

6 0.72130005E+0 2+0 

7 0. 72130005E+02+O 

8 0.721 3 00 0 5 E +02+0 

9 0.721 30005 E +02+0 
Y PUMP ( I NCL CAPAC ) = 
NUMBER OF POINTS READ = 


0.99999994E-10OHMS 
0.99999994E-10OHMS 
0.99999994E-10OHMS 
0. 99999994E-10OHMS 
0.99999994E-10OHMS 
0.99999994E-10OHMS 
0.99999994E- 10OHMS 
0. 4787 1 843E-02 +0 -0. 8007S605E-02S 
64 


INPUTTING SIDEBAND TERMINATION INFORMATION 
THE FIRST 10 SIDEBAND IMPEDANCES I EXC CAPAC): 

0 0.50000000E+02+O 0.0 OHMS 

1 0 . 5 5000000E+02 +0 0. 92000000E+02OHMS 

-1 0. 55000000E+02+O 0. 92000000E+02OHMS 

2 0.721 30005 E +02+0 0. 99999994E-10OHMS 

-2 0.721 3000 5 E +02+0 0 . 9 99 9 9 9 9 4 E - 10OHMS 

3 0.721 30005 E +02+0 0 . 9 9 9 9 99 9 4E - 10OHMS 

-3 0. 721300D 5E+02+O 0 . 9 9 9 9 9 9 9 4 E - 10OHMS 

4 0.72130005E+02+O 0 . 999 9 9 9 9 4E - 1 0OHMS 

-4 0.72130005E+02+O 0 . 9 9 9 9 9 9 9 4 E- 1 0OHMS 

5 0.721300C5E+02+O 0 . 999999 94E - 10OHMS 

Y SIGNAL (INCL CAPAC) = 0. 4787 1 843E-02 +0 -0. 80076605E-02S 

Y IMAGE (INCL CAPAC) = 0. 4787 1 843E-02 +0 0. 80076605E-02S 

NUMBER OF POINTS READ = 127 


********************************* 
*r*********£^C££ J*************** 


EMBEDDING THEVENIN VDC SOURCE = 0. 27999997E-02 VOLTS 

EMBEDDING THEVENIN VLO SOURCE (PEAK) = 0. 34999999E-02 VOLTS 

INCIDENT LOCAL OSCILLATOR POWER = 0. 27840887E-07 WATTS 


GENERAL NONLINEAR ANALYSIS USING VOLTAGE UPDATE METHOD 
BEGINNING THE ITERATION PHASE 


ITERATION NUMBER = 1 

ERR = 0. 10000000E+01POSN 

NO = 

1 

ITERATION NUMBER = 2 

ERR = 0. 13908285E+00POSN 

NO = 

74 

ITERATION NUMBER = 3 

ERR = 0. 78462720E-01POSN 

NO = 

85 

ITERATION NUMBER = 4 

ERR = 0 . 5975 3280E-01 POSN 

NO = 

92 

ITERATION NUMBER = 5 

ERR = 0. 16621690E-01POSN 

NO = 

99 
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FILEs SISCAP DOC A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


ITERATION NUMBER = 6 

ERR = 0. 53864606E-02POSN NO = 85 

ITERATION NUMBER = 7 

ERR = 0. 31008909E-02POSN NO = 93 

ITERATION NUMBER = 8 

ERR = 0.98307454E-03POSN NO = 97 

ITERATION NUMBER = 9 

ERR = 0.3879G710E-03POSN NO = 93 

ITERATION NUMBER = 10 

ERR = 0. 18293827E-03POSN NO = 94 


ITERATION NUMBER = 11 

ERR = 0 . G 1 9G22 1 5E-04POSN NO = 95 

LARGE SIGNAL ANALYSIS COMPLETED IN 11 ITNS 

POINT NO. = 1 VOLTAGE = 0 . 36563654E-02CURRENT 

POINT NO. = 2V0LTAGE = 0. 3G8 1 2625E-02CURRENT 

POINT NO. = 3V0LTAGE = 0 . 3 7026 438E-02CURRENT 

POINT NO. = 4V0LTAGE = 0. 37205843E-02CURRENT 

POINT NO. = 5V0LTAGE = 0. 3735 1 9 1 4E-02CURRENT 

POINT NO. = GVOLTAGE = 0. 3746597 1 E-02CURRENT 

POINT NO. = 7V0LTAGE = 0. 37549580E-02CURRENT 

POINT NO. = 8V0LTAGE = 0 . 37604398E-02CURRE NT 

POINT NO. = 9V0LTAGE = 0. 37632226E-02CURRENT 

POINT NO. = 10VOLTAGE = 0 . 37G348 1 3E-02CURRE NT 

POINT NO. = 1 1 VOLTAGE = 0.3761 39 4 6E-02CURRENT 

POINT NO. = 1 2V0LTAGE = 0 . 3757 1 308E-02CURR E NT 

POINT NO. = 1 3V0LTAGE = 0. 37508472E-02CURRENT 

POINT NO. = 1 4V0LTAGE = 0. 3742G832E-02CURRENT 

POINT NO. = 1 5V0LTAGE = 0. 37327G8 1 E-02CURRENT 

POINT NO. = 1 GVOLTAGE = 0.3721 2060E-02CURRENT 

POINT NO. = 1 7V0LTAGE = 0. 37080892E-02CURRENT 

POINT NO. = 1 8V0LTAGE = 0. 36934826E-02CURRENT 

POINT NO. = 1 9V0LTAGE = 0. 36774450E-02CURRENT 

POINT NO. = 20VOLTAGE = 0. 3GG001 23E-02CURRENT 

POINT NO. = 2 1 VOLTAGE = 0. 3G4 1 2 100E-02CURRENT 

POINT NO. = 22V0LTAGE = 0. 3G2 10504E-02CURRENT 

POINT NO. = 23VOLTAGE = 0. 3599540GE-02CURRENT 

POINT NO. = 2 4V0LTAGE = 0.357GG7G 1 E-02CURRENT 

POINT NO. = 2 5V0LTAGE = 0. 35524534E-02CURRENT 

POINT NO. = 26VOLTAGE = 0. 35268620E-02CURRENT 

POINT NO. = 27VOLTAGE = 0. 34998977E-02CURRENT 

POINT NO. = 28V0LTAGE = 0. 347 1 5536E-02CURRENT 

POINT NO. = 29VOLTAGE = 0. 344 1 8299E-02CURRENT 

POINT NO. = 30VOLTAGE = 0 . 34 10729 1 E-02CURR ENT 

POINT NO. = 3 1 VOLTAGE = 0 . 33 782658E-02CURRENT 

POINT NO. = 32VOLTAGE = 0 . 33 4 44576E-02CURRENT 

POINT NO. = 33VOLTAGE = 0 . 330933 G5E-02CURRENT 

POINT NO. = 34VOLTAGE = 0. 32729360E-02CURRENT 

POINT NO. = 35VOLTAGE = 0. 32353 105E-02CURRENT 


0.35102436E-04 
0.3G557831E-04 
0. 37954203E-04 
0. 39285034E-04 
0. 4053220GE-04 
0. 41695879E-04 
0. 427538G2E-04 
0. 437 1 565GE-04 
0.445G2665E-04 
0. 45304492E-04 
0. 45932538E-04 
0. 4645 1751 E-04 
0.4G863002E-04 
0.47166 1 92E-04 
0. 47370107E-04 
0.47475507E-04 
0. 4748G261E-04 
0.4741 5029 E -04 
0. 47253392E-04 
0. 47026 100E-04 
0. 46722096E-04 
0.463551 1 1 E -04 
0. 45924549E-04 
0.4 5438268E-04 
0. 44896689E-04 
0. 44308676E-04 
0.43670356E-04 
0. 4299 3066E-04 
0.42271 175E-04 
0.41 509833E-04 
0. 40709 1 72E-04 
0. 3987 1 244E-04 
0. 38996906E-04 
0.38083905E-04 
0. 37 1 37594 E-04 
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FILE: SISCAP DOC 


POINT 

NO. 

c 

POINT 

NO. 

= 

POINT 

NO. 

s 

POINT 

NO. 

£= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

- 

POINT 

NO. 


POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 


POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

S5 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 


POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 


POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

PO I N'T 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 


POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 

POINT 

NO. 

= 


36V0LTAGE 
37VOLTAGE 
38V0LTAGE 
39VOLTAGE 
40VOLTAGE 
4 1 VOLTAGE 
4 2V0LTAGE 
43V0LTAGE 
4 4V0LTAGE 
4 5V0LTAGE 
4 6V0LTAGE 
47VOLTAGE 
4 8V0LTAGE 
4 9V0LTAGE 
50VOLTAGE 
E1V0LTAGE 
52V0LTAGE 
53VOLTAGE 
54V0LTAGE 
55V0LTAGE 
56V0LTAGE 
57V0LTAGE 
58V0LTAGE 
59V0LTAGE 
60VOLTAGE 
G 1 VOLTAGE 
G2V0LTAGE 
63V0LTAGE 
64VOLTAGE 
65VOLTAGE 
66VOLTAGE 
67V0LTAGE 
G8V0LTAGE 
69V0LTAGE 
70VOLTAGE 
7 1 VOLTAGE 
72V0LTAGE 
73VOLTAGE 
74VOLTAGE 
75VOLTAGE 

7 GVOLTAGE 
77V0LTAGE 
78V0LTAGE 
79VOLTAGE 
80VOLTAGE 
8 1 VOLTAGE 
82V0LTAGE 
f 3V0LTAGE 
84VOLTAGE 
85VOLTAGE 
86VOLTAGE 

8 7V0LTAGE 
D8V0LTAGE 
69VOLTAGE 
S0VOLTAGE 
9 1V0LTAGE 
92VOLTAGE 
S3V0LTAGE 

9 4V0LTAGE 
95V0LTAGE 
96V0LTAGE 
97VOLTAGE 
98VOLTAGE 
S9V0LTAGE 

100VOLTAGE 

101 VOLTAGE 

102VOLTAGE 


NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


0.31 9G5207E-02CURRENT 
0.31 566401 E-02CURRENT 
0.31 157513E-02 CURRENT 
0. 30739547 E-02CURRE NT 
0.3031 35 74E-02CURRENT 
0.2988081 9 E-02CURRENT 
0.294 425 57E-02CURR ENT 
0. 2 9000205 E- 02 CUR RE NT 
0.28 55524 7 E-02CURR ENT 
0.281092G2E-02CURRENT 
0.27663852E-02CURRENT 
0.272 20668E-02CURRENT 
0.26781 39 1E-02CURRENT 
0. 26 3 47 6 95E-02CURRENT 
0. 25 92121 5E-02CURRENT 
0.25503556E-02CURRENT 
0. 2509G249E-02CURRENT 
0.24700745E-02CURRENT 
0.24318350E-02CURRENT 
0. 2395027 1E-02CURRENT 
0.23597558E-02CURRENT 
0.232G1076E-02CURRENT 
0.22941S44E-02CURRENT 
0.226394 5 1E-02CURRENT 
0. 22 3 55 143E-02CURRE N'T 
0.22088757E-02CURRENT 
0. 2 1 84024 1 E-02 CURRENT 
0.21 609 364E-02CURRENT 
0.21395762E-02CURRENT 
0. 21 198881E-02CURRENT 
0.210181 25E-02CURRENT 
0. 208527 18E-02CURR ENT 
0.20701 862E-02CURR ENT 
0. 205647 18E-02CURRENT 
0. 204 4042 1E-02CURRENT 
0.203281 18E-02CURR ENT 
0.2022701 8 E-02CURRENT 
0.201 3635 8 E-02CURRENT 
0 . 20055503E-02CURR ENT 
0. 1S983889E-02 CURRENT 
0. 1992 1 06 9E-02 CURRENT 
0. 19866740E-02CURRENT 
0. 19820742E-02CURRENT 
0. 1978303 1E-02CURRENT 
0. 19753748E-02CURRENT 
0. 19733123E-02CURRENT 
0. 197216 38E-02CURRENT 
0. 1S719815E-02CURRENT 
0. 1S728362E-02CURRENT 
0. 197481 22E-02CURRENT 
0. 19 7 8007 4E-02CURRENT 
0. 1S8252G6E-02 CURRENT 
0. 19 884 6 9 9E-02CURRENT 
0. 19960220E-02CURRENT 
0.20052569E-02CURRENT 
0. 201 6329 2E-02CURRENT 
0.2029 37 87E-02CURR ENT 
0.20445397E-02CURRENT 
0.2061 94 68E-02CURRENT 
0.2081 72 27E-02CURR ENT 
0.21039834E-02CURRENT 
0.21 2882 60E-02CURRENT 
0.2 15634 03 E-02CURRENT 
0.21865845E-02CURRENT 
0.22196025E-02CURRENT 
0. 2255408 1E-02CURRENT 
0. 2293995 1E-02CURRE NT 


0.361 54655E-04 
0.351 352 94E-04 
0. 3408 1 284E-04 
0.32992626E-04 
0.318671 6 6E-04 
0.30705502E-04 
0.29508359E-04 
0.28276278E-04 
0 . 27009 1 9 9E-04 
0. 25706700E-04 
0.24371 766E-04 
0. 2 3002 402 E -04 
0.21 6069 1 7E-04 
0 . 201 78020E-04 
0. 18730527E-04 
0. 1725S368E-04 
0. 15773418E-04 
0. 14274468E-04 
0. 12772449E-04 
0. 1 1269207E-04 
0 . 97750744E-05 
0. S2930055E-05 
0.6837581 2E-05 
0.54113771 E-05 
0. 4025991 IE-05 
0 . 2 685 104 8 E-05 
0. 1 4031393E-05 
0. 18249557E-06 
-0. 96604344E-06 
-0. 20427 105E-05 
-0. 30309020E-05 
-0.39313018E-05 
-0.47391368E-05 
-0.54530101 E-05 
-0. 60664S97E-05 
-0.6579E959E-05 
-0. 69922826E-05 
-0. 7308E 1 58E-05 
-0.75217 104 E-05 
-0. 7642E259E-05 
-0.7662441 9 E-05 
-0. 75992303E-05 
-0.74446871 E-05 
-0.721721 15E-05 
-0. 6 902 904 7 E-05 
-0.6531 7 727 E-05 
-0.60873426E-05 
-0. 5 59327 60E-05 
-0.5042S062E-05 
-0. 445383 19E-05 
-0. 38256440E-05 
-0.31720265E-05 
-0.24886249E-05 
-0. 1 8043793E-05 
-0. 1 09 98 625 E-05 
-0.39467704E-06 
0.31491 800E-06 
0. 10105286E-05 
0. 16979029E-05 
0.23706962E-05 
0.30339870E-05 
0. 367 1 5355E-05 
0. 42979227E-05 
0. 4908G666E-05 
0. 54999255E-05 
0. G0887687E-05 
0.66628563E-05 
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POINT NO. = 103VOLTAGE = 0. 23353205E-02CURRENT = 0. 7236 1 436E-05 

POINT NO. = 104VOLTAGE = 0. 23793 1 89E-02CURRENT = 0. 78079G 1 1 E-05 

POINT NO. = 1£ 5V0LTAGE = 0 . 24 25 8 88 1 E-02CURR ENT = 0. 83870545E-05 

POINT NO. = 106VOLTAGE = 0. 24749034E-02CURRENT = 0. 897G1 743E-05 

POINT NO. = 107VOLTAGE = 0 . 252G202 7E-02CURRENT •= 0. 958G5G89E-05 

POINT NO. = 108VOLTAGE = 0. 25795978E-02CURRENT = 0. 102 1 8384E-04 

POINT NO. = 1 jf 9V0LTAGE = 0. 2G34 8GG9E-02CURRENT = 0 . 10884 57 1 E-04 

POINT NO. = 1 10VOLTAGE = 0 . 2G9 1 7 G 65E-02CURR ENT = 0. 1 1 585973E-04 

POINT NO. = 111 VOLTAGE = 0. 2750022GE-02CURRENT = 0. 1 2333083E-04 

POINT NO. = 1 1 2V0LTAGE = 0. 2S093394E-02CURRENT = 0. 1 3 1 29320E -04 

POINT NO. = 1 1 3V0LTAGE = 0. 28693988E-02CURRENT *= 0 . 1 3978235E-04 

POINT NO. = 1 1 4V0LTAGE = 0. 29298 668E-02CURR ENT = 0. 1 4888705E-04 

POINT NO. = 1 1 5V0LTAGE = 0. 29903932E-02CURRENT = 0. 1 58G2883E-04 

POINT NO. = 1 1 6V0LTAGE = 0 . 3050G 1 90E-02CURR ENT - 0. 1 G903679E-04 

POINT NO. = 1 1 7V0LTAGE = 0.31 101 780E-02CURRENT = 0. 1 80108G2E-04 

POINT NO. = 1 1 8V0LTAGE = 0 . 3 1 687077E-02CURRENT * 0. 19187792E-04 

POINT NO. = 1 1 9V0LTAGE = 0. 32258469E-02CURRENT = 0. 20430G85E-04 

POINT NO. = 1 20VOLTAGE = 0 . 32 8 1 24 64E-02CURRENT = 0. 21738233E-04 

POINT NO. = 1 2 1 VOLTAGE = 0. 33345749E-02CURRENT = 0. 23 103901 E-04 

POINT NO. = 1 22VOLTAGE = 0 . 33 855 2 1 9E-02CURR ENT = 0. 24523295E-04 

POINT NO. = 1 23V0LTAGE = 0. 34 338073E-02CURRENT = 0. 25988586E-04 

POINT NO. = 1 24VOLTAGE = 0 . 3 4 79 1 802E-02C URR ENT = 0 . 2 7487 957E -04 

POINT NO. = 1 25VOLTAGE = 0 . 3 52 1 4 2 73E-02C URRENT = 0 . 2901 3652E-04 

POINT NO. = 1 2GV0LTAGE = 0. 35603768E-02CURRENT = 0. 30549549E-04 

POINT NO. = 1 27V0LTAGE = 0 . 35 9 59005E -02C URR E NT = 0. 32087453E-04 

POINT NO. = 1 28VOLTAGE = 0. 36279093E-02CURRENT = 0. 336 101 5 1 E-04 

DC VOLTAGE AT THE JUNCTION = 0. 27999 97 1 E-02 VOLTS 

LO VOLTAGE (PEAK) AT THE JUNCTION = 0. 92 1 97349E-03 VOLTS 

SMALL SIGNAL AND NOISE ANALYSIS SECTION 


DETERMINING SMALL SIGNAL MIXER PROPERTIES USING THE GENERAL ANALYSIS 

IF OUTPUT IMPEDANCE = 0. 548 1 8091 E+03 +J -0. 2379G079E-02 

IF VSWR REFERRED TO 50.00 OHMS = 0. 10963G 1 3E+02 

OUTPUT MISMATCH GAIN = 0. 30G4001 2E+00 ( -0. 5 1 37 1098E+01 DB ) 

INPUT SIGNAL IMPEDANCE = 0. 852G2G50E+02 +J 0. 93G24878E+01OHMS 

INPUT MISMATCH GAIN = 0. 82252783E+00 ( -0. 84849358E+00DB ) 

INPUT GAIN CALCULATION REFERRED TO = 0. 20889 105E+03 OHMS 

CONVERSION GAIN (USB) EQUALS -0. 3022932 1 E+01 DB 
CONVERSION GAIN (LSB) EQUALS -0. 302292G3E+01 DB 

SSB TEMP OF MIXER (USB) = 0. 23333755E+02 DEG K 
SSB TEMP OF MIXER (LSB) = 0. 233 33 7 10E+02 DEG K 
R; T=244. 59/252. 26 14:03:14 


FILE: PARAM DAT A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


2 


ISWIT 

1 


NODAT 


0.025E-3 

DEL V 


2 . 8000E-3 

VSTART 

3 

. 50000 E -3 

VLOSTR 

0 

. 10000E-3 

DELVLO 

1 


NLO 


1 .0E-4 

TOL 


1 .0E-4 

EPS 

128 


NOPNTS 

7 


I BASE 1 

1024 


NO I V 

10 


IBASE2 

512 


NCCHI 

50 


NCPER 

31 


m:;ax 

31 


NMAX 


72.131729 

RN 

19 

. 3720E- 1 5 

CAP 


2.7 

TEMP 


1 13.9E9 

FREQ 


1 .0E6 

F IF 


1 .0 

P 


50.000 

RID 


50.0 

Z0 
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FILE: ZSSB DATA A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 

127 



50.0 

0.0 

55 

. 00000 

92.00000 

55 

. 00000 

92.00000 


72.13 

1 .E-10 


72.13 

1 .E-10 


72.13 

1 .E-10 


72.13 

1 .E-10 


72.13 

1 .E-10 


72.13 

1 .E-10 


72.13 

1 .E-10 


72.13 

1 .E-10 


72.13 

1 .E-10 


72.13 

I .E-10 


72.13 

1 .E-10 


72.13 

1 .E-10 


72.13 

1 .E-10 


72. 13 

1 .E-10 


72.13 

1 .E-10 


72.13 

1 .E-10 


72.13 

1 .E-10 


FILE: ZLO DATA A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


0.0 

0.0 

55.00 

92.00 

72.13 

1 .E-10 

72.13 

1 .E-10 

72.13 

1 .E-10 

72.13 

1 .E-10 

72.13 

1 .E-10 

72.13 

1 .E-10 

72.13 

1 .E-10 

72.13 

1 .E-10 

72. 13 

1 .E-10 

72.13 

1 .E-10 

72. 13 

1 .E-10 

72.13 

1 .E-10 

72.13 

1 .E-10 

72.13 

1 .E-10 

72.13 

1 .E-10 

72.13 

1 .E-10 

72.13 

1 .E-10 

72.13 

1 .E-10 
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A. 3 Typical Run of Program SISCAP 

The following pages show a typical run of program SISCAP. It 
utilizes the CHI and DC I-V data produced from the run of GETCHI shown in 
the previous section. This program uses as input the following files: 

CHI data (unit no. 4) , processed DC I-V data (unit no. 8) , control 
parameters (unit no. 11), sideband termination impedances (unit no. 14) 
and LO harmonic impedances (unit no. 15). One output file is generated 
via the console. 

A. 4 Memory Allocation Requirements of GETCHI 

Memory is allocated in program GETCHI in accordance with the 
following rules. NUM is preset by subroutine EVEN to be 4096. NOCHI and 
NOIV are numbers which the user obtains by running program GETCHI. Since 
the user needs the value of these parameters prior to running the program 
(in order to allocate the correct amount of memory for the program), the 
program may have to be run twice, the first time with an excessive amount 
of memory in order to obtain the values of NOCHI and NOIV. Table 2 below 
indicates the correct memory allocation required and in which subroutine 
the memory space has to be declared. 


VARIABLE 

SUBROUTINE 

SIZE 

V0LT2 

MAIN 

NUM 

CUR2 

MAIN 

NUM 

V0LT3 

MAIN 

NOCHI 

CUR3 

MAIN 

NOCHI 

V0LT4 

MAIN 

NOIV 

CUR4 

MAIN 

NOIV 


TABLE 2 
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VARIABLE 

SUBROUTINE 

SIZE 

TIME 

MAIN 

NOCHI 

CHI 

MAIN 

NOCHI 

CCHI 

MAIN 

2*N0CHI 

FACTOR 

MAIN 

2*N0CHI 


TABLE 2 

A. 5 Memory Allocation Requirements of SISCAP 

This section outlines the memory allocation requirements of program 
SISCAP. There are 5 key parameter values which must be specified prior to 
allocating the appropriate memory, namely NOPNTS, NOIV, NOCHI, MMAX and 
NMAX. 

NOPNTS and NOCHI are determined by the discretization constraints on 
CHI. NOPNTS is the number of points per LO cycle and NOCHI is the number 
of CHI points input to SISCAP. Both are provided by using the program 
GETCHI. NOPNTS must be a power of 2 and NOCHI, in addition to being a 
power of 2, must be equal to or greater than NOPNTS. 

NOIV is the number of points in the I-V curve input to the program 
SISCAP. It must also be a power of 2 and must be large enough such that 
the discretization of the I-V curve is adequate. This value is obtained 
from program GETCHI. 

MMAX indicates the number of LO sidebands which are taken into 
account. MMAX must be odd and is limited to NOPNTS/2 - 1 in value. See 
Section 3.4 item (j) for a more detailed discussion. 
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NMAX (assumed odd) Indicates how many terms are significant in the 
series which are used to calculate the small signal conversion properties 
of the mixer. Its value should be as discussed in item (k) of Section 

3.4. 

Having established these values, the memory allocation is as per 
Table 4. Note that the symbol '*' indicates a multiplication whereas the 
symbol 'x' indicates a multidimensional array. 


VARIABLE 

SUBROUTINE 

SIZE 

ZIMPED 

MAIN 

NOPNTS/2 

VSOURG 

MAIN 

NOPNTS/2 

CVAL 

MAIN 

NOPNTS 

VOLT 

MAIN 

NOPNTS 

CUR 

MAIN 

NOPNTS 

OLDVOL 

MAIN 

NOPNTS 

CTEMP 

LINEAR 

NOPNTS 

W 

TOPS 

NOPNTS 

CVOL 

SIGNL 

NOPNTS 

CVOL 

SISNL 

NOPNTS 

FACTOR 

MAIN 

2*MAX OF (NOIV & NOCHI) 

J 

MAIN 

NOIV 

IVDAT 

MAIN 

NOIV 

CCHI 

SISNL 

N0CHI*2 

CTEMP 

GETJ 

NOIV* 2 


TABLE 4 
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VARIABLE 

ZSSB 

CHI 

UTERM 

YSIG 

NOISE 

ZSIG 

B 

STORE 

STORE 

YSI2 

NOI2 

ZSI2 

VAL 

STORE 

ESOURC 

ELEFT 

E RIGHT 

RHO 

CCUR 

ZTERM 

PHASE 

VOLT2 

UTERM 

CVOL 


SUBROUTINE 

MAIN 

MAIN 

SISNL 

TOPS 

TOPS 

TOPS 

TOPS 

TOPS 

PROCES 

TOPS 

TOPS 

TOPS 

THEORY 

TERP5 

MRT 

MRT 

MRT 

MRT 

SISCHA 

SISCHA 

SISCHA 

SISCHA 

SISCHA 

SISNL 


SIZE 

NOPNTS-1 
NOCHI 
NOCHI*2 
MMAX x MMAX 
MMAX x MMAX 
MMAX x MMAX 
MMAX 

MMAX x MMAX 

MMAX 

3x3 

3x3 

3x3 

NMAX+10 

4 

NOPNTS 

NOPNTS 

NOPNTS 

NOPNTS 

NOPNTS 

NOPNTS 

NOCHI 

NOCHI 

NOCHI 

NOPNTS 


TABLE 4 CONT 
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A. 6 Listing of Program SISJUN 

This section gives the listing of the test data program SISJUN. 

FILEs SISJUN FORTRAN A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


THIS PROGRAM PROVIDES A STANDARD TEST DC I-V CURVE 
IN WHICH TO VERIFY THE CORRECT FUNCTIONING OF THE 
PROGRAMS GETCHI AND SISCAP 

I-V CURVE FOR FUNCTION PNBS W2 C2736 D4 


VI = 0.0 
V2 = 2.44 
V3 = 5.99 
V4 = 9.20 
C 

NOPNTS = 2048 
DELTAV = 0.01 
C 

WR ITE ( 7 , 1 9 ) NOPNTS 
19 FORMAT (15) 

C 

VSPAN = FLOAT! NOP NTS -1 ) * DELTAV 
C 

11 = 1 

1 2M1 = V2 / DELTAV + 0.5 

12 = I 2M1 + 1 

I3M1 = V3 / DELTAV + 0.5 

13 = I3M1 + 1 

I 4M1 = V4 / DELTAV + 0.5 

14 = I4M1 + 1 
C 

DO 100 I = I1.I2M1 

VOUT = (FLOAT! I - 1 ) * DELTAV) 

CUROUT = (VOUT/3.0) * 1.0E-6 / 2.0 
VOUT = VOUT * 1.0E-3 / 2.0 
WRITE(7, 1 10) VOUT, CUROUT 
110 FORMAT ! 2E 1 5 . 8 ) 

100 CONTINUE 
C 

DO 200 I = I2.I3M1 

VCUT = l FLOAT! I - 1 > * DELTAV) 

CUROUT = ( 0 . 9*VOUT - 1.37 + EXP ! ! VOUT-4 . 98 > *3 . 0775 ) + 

1 EXP ! ( VOUT-5 . 749 )» 12. E 1 ) ) * 1 .0E-6 / 2.0 
VOUT = VOUT * 1.0E-3 / 2.0 
WR I TE ( 7 , 1 1 0 ) VOUT, CUROUT 
200 CONTINUE 
C 

DO 300 I = 13, I4M1 

VOUT » (FLOAT! 1-1) * DELTAV) 

CUROUT = ( 15.4325'VOUT - 25.33 - EXP ! ( 7 .08G-VOUT > *2 .01 25 ) 
1 EXP! (S.27-VOUT)*E.45) ) * 1 .0E-6 / 2.0 
VOUT = VOUT * 1.0E-3 / 2.0 
WP I TE ( 7 , 1 1 0 ) VOUT, CUROUT 
300 CONTINUE 
C 

DC 400 I = 14, NOPNTS 
VCUT = ( FLOAT! 1-1 )*DELTAV> 

CUROUT = (14.1 637*VOUT - 125.78/VOUT) * 1 .0E-6 / 2.0 
VOUT = VOUT * 1.0E-3 / 2.0 
WR ITE ( 7 , 1 10) VOUT, CUROUT 
400 CONTINUE 
C 

STOP 

END 


SIS00010 
S IS00020 
S IS00030 
SIS00040 
SIS00050 
SIS00060 
S IS00070 
SIS00080 
SIS00090 
SIS00100 
SIS001 10 
SIS00120 
SIS00130 
SIS00140 
SIS00150 
SIS00160 
SIS00170 
SIS00180 
SIS00190 
SIS00200 
SIS00210 
SIS00220 
SIS00230 
SIS00240 
S IS00250 
SIS00260 
SIS00270 
SIS00280 
S IS00290 
SIS00300 
SIS00310 
S IS00320 
S IS00330 
SIS00340 
SIS00350 
SIS00360 
S IS00370 
SIS00380 
S IS00390 
SIS00400 
SIS00410 
S IS00420 
SIS00430 
SIS00440 
SIS00450 
SIS004G0 
SIS00470 
SIS00480 
SIS00490 
SIS00500 
SIS00510 
SIS00520 
SIS00530 
SIS00540 
SIS00550 
S IS00560 
S IS00570 
SIS00580 
S IS00590 
SIS00600 
SIS00610 
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A. 7 Listing of Program GETCHI 


This section gives the listing of the preprocessing program GETCHI. 


FILE: GETCHI 


FORTRAN A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


REAL VOLT2 ( 4096 ) , VOLT3(512>, VOLT4U024) 

REAL CUR2U096), CUR3(512>, CUR4U024) 

REAL TIME ( 51 2 ) , CHK512) 

COMPLEX CCHIU024), FACTOR(1024) 

COMMON /GAP INF/ VGAP 

SUBROUTINE GETCHI PRODUCES THE CHI FUNCTION WITH THE 
APPROPRIATE SPECIFIED DISCRETIZATION. THIS PROGRAM 
REQUIRES AS INPUT THE I-V CURVE OF THE SIS JUNCTION . 

OUTPUT IS TWO-FOLD: 

(I) THE TABULATED CHI FUNCTION {UNIT 8), AND 

(II) A COMPATIBLE TABULATED DISCRETIZATION OF THE 
I-V DATA (UNIT 9). 

THE INPUT I-V DATA WILL BE CHECKED FOR THE FOLLOWING REQUIREMENTS 

(A) ThE FIRST DATA POINT MUST BE THE ORIGIN; 

(B) THE I-V FUNCTION MUST BE MONOTONIC; 

<C> THE RANGE OF INPUT DATA MUST BE SUCH THAT THE VOLTAGE 

RANGE EXTENDS FROM ZERO TO VOLTAGES WELL INTO THE LINEAR 
RANGE. THE INPUT DATA DOES NOT HAVE TO BE EVENLY DISCRETIZED 


C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

WR ITE ( 6 , 2000 ) 

2000 FORMAT! /IX , ’THIS PROGRAM, GETCHI, PERFORMS A PREPROCESSING’, 
1/1X, ’FUNCTION ON THE INPUT DC I-V CURVE DATA PRIOR TO USING’, 
2/ 1 X , 'THE PROGRAM SISCAP. OUTPUT FROM GETCHI IS THE PRE-’, 
3/1X, 'PROCESSED DC I-V DATA, CHI IN AN APPROPRIATELY DIS-’, 
4/1X, ’CRETIZED FORM AND VALUES FOR THE PARAMETERS: NOCHI,’, 
5/1X, ’NOIV, NOP NTS AND RN’) 


C 

20 


25 

C 


C 

C 

c 

c 

c 

10 

c 

40 

30 


WRITE (6, 20) 

F0RMAT(/1X, 'INPUT THE GAP VOLTAGE (VGAP) IN VOLTS OF THE SINGLE' 
1/10X, ’JUNCTION EQUIVALENT DEVICE (FORMAT E15.8)’) 

READ (5,25) VGAP 
FORMAT ( E 1 5 . 8 ) 

CALL EVEN(VOLT2,CUR2,NUM,DELTAV) 

RSTART = 1.0E9 

CALL REFINE(VOLT2,CUR2,NUM, RSTART, RFINIS) 

CALL INEXP(VOLT2,CUR2,NUM,VOLT3,CUR3, NOCHI, NOPNTS) 

RSTART = RFINIS 

CALL REF I NE(VOLT3,CUR3, NOCHI .RSTART, RFINIS) 

CALL IVTOCH(VOLT3,CURC,NOCHI , DELTAV, TIME, CHI ,CCHI .FACTOR) 

CALL REEXP (VOLT2.CUR2, NUM.VOLT4 ,CUR4, NOIV, DELTAV) 

WR ITE ( 8 , 10 ) NOCHI 
FORMAT! 15) 

DO 40 I = 1, NOCHI 

WR I TE ( 8 , 30 ) TIME(I), CHI(I) 

F ORMAT ( 2E 1 5 . 8 > 


GET00010 
GET00020 
GET00030 
GET00040 
GET00050 
GET00050 
GET00070 
GET! 0080 
C ET r 0090 
GETC0100 
GET001 10 
GET001 20 
GET00130 
GET001 40 
: GET00 1 50 
GET001S0 
GET001 70 
GET001 80 
GET00190 
. GET00200 
GET00210 
GET00220 
GET80230 
GET00240 
GET00250 
GET002S0 
GET00270 
GET00280 
GET00290 
GETX0300 
GET003 10 
GET00320 
GET00330 
GET00340 
GET00350 
GET00360 
GET! 0370 
GET00380 
GETL0390 
GET0040 v 
GETL04 a. 
GETL 0420 
GET! 0430 
GET00440 
GET00450 
GET00460 
GET00470 
GET0 04 80 
GETP0490 
GET> 0500 
GETC0510 
GET00520 
GETC0530 
GET00540 
GET00550 
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FILE: GETCHI 


FORTRAN A 


NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


C 

C 

50 

C 

60 


C 


C 

c 

c 

c 

c 

c 


c 

10 

2 

C 

c 

C5 

C 

C 

C 

C 

c 

20 

30 

C 

C 

C 

C 


C 


WR ITE ( 9 , 10) NOIV 

DO 50 I = 1 , NOIV 

WR ITE ( 9 , 30 ) VOLT4 ( I > , CUR4(I) 

WRITE (6, 60) NOCHI, NOIV, NOPNTS, RFINIS 

F0RMAT(//1X, ’ INFORMATION FOR TRANSFERRAL TO PROGRAM SISCAP:’ 

1/1X, ’OUTPUTTING THE CHI FILE ON UNIT 8: NOCHI = ’,15/ 

21X, ’OUTPUTTING THE IV FILE ON UNIT 9: NOIV = ’,15/ 

3 1 X , ’NOPNTS = ’ , 15/ 

4 IX, ’RN = ’ ,E15.8, ’ OHMS’ ) 

STOP 

END 

SUBROUTINE EVEN ( VOLTV ,CURV, NUM , DELTAV ) 

REAL CURV(NUM), VOLTV(NUM) 

LOGICAL FLAG 
COMMON /GAPINF/ VGAP 

SUBROUTINE EVEN: THIS SUBROUTINE PROVIDES AN EVENLY SPACED 
VERSION OF THE EXPERIMENTALLY OBTAINED INPUT I-V DATA. 

LINEAR INTERPOLATION IS USED TO PPOVIDE THE OUTPUT POINTS 
THE EXPERIMENTAL DATA IS SCALED TO A SINGLE JUNCTION EQUIVALENT 

DELTAV = 0.00005 
NUM = 4096 

READ ( 7 , 10, END=2000 ) NTMP 
FORMAT ( 15) 

WP I TEC 6, 2) 

FORMAT? /IX, ’ INPUTTING THE INPUT I-V CURVE DATA FROM UNIT 7’) 

WR I TE ( 6 , 5 ) 

F0RMAT(/1X, ’SUBROUTINE EVEN: THIS SUBROUTINE PROVIDES AN EVENLY’ 
1 1 X , 'SPACED VERS I ON ’ / 1 8X , ’OF THE EXPERIMENTALLY OBTAINED INPUT I- 
2,’ DATAV18X, ’LINEAR INTERPOLATION IS USED TO PROVIDE THE POINTS 
3 / 1 8X , ’ON OUTPUT, THE EXPERIMENTAL DATA IS SCAL ED ’ / 1 8X , ’ TO A’, 

4’ SINGLE JUNCTION EQUIVALENT’) 

WR ITE ( 6 , 20 ) 

FORMAT1/1X, ’ INPUT THE NUMBER OF JUNCTIONS (NOJUN) IN DEVICE (15)’ 
READ ( 5,30) NOJUN 
FORMAT! 15) 

ADJUST DELTAV TO CORRESPOND WITH A MULT I - JUNCTI ON DEVICE 

DELTAV = FLOAT ( NOJUN ) * DELTAV 

VOLT = 0.0 
CUR = 0.0 
VI = 0.0 
Cl = 0.0 
FLAG = .FALSE. 


GET00560 
GET00570 
GET00580 
GET00590 
GET00600 
GET006 10 
GET00620 
GET00630 
GET00640 
GET00650 
GET00660 
GET00670 
GET00680 
GET00690 
GET00700 
GET007 10 
GET00720 
GET00730 
GET00740 
GET00750 
GET00760 
GET00770 
GET007 80 
GET007 90 
GET00800 
GET008 10 
GET00820 
GET00830 
GET00840 
GET00850 
GET00860 
GET00870 
GET00880 
GET00890 
, GET00900 
VGET009 10 
’GET00920 
GET00930 
GET00940 
GET00950 
GET00960 
>GET00970 
GET0098C 
GET00990 
GET01000 
GET01010 
GET01020 
GET01030 
GET01040 
GET01050 
GET01060 
GET01070 
GET01080 
GETB 1090 
GET01 100 
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FILE: GETCHI 


FORTRAN A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 



READ (7,40) V2 , C2 


GET0U10 

40 

F ORMAT { 2E 1 5 . 8 ) 


GET01 120 

C 



GET01 130 


I F ( V2 . N£ .0.0 .OR. C2.NE.0.0) GO TO 3000 


GET01 140 

C 



GET01 150 


DO 100 I = 1 , NUM 


GET01 160 


VOLTV(I) = VOLT / FLOAT { NOOUN ) 


GET01 170 


CURV(I) = CUR / FLOAT { NOJUN ) 


GET01 180 


VOLT = VOLT + DELTAV 


GET01 190 

50 

I F ( VOLT . LT . V2 ) GO TO 60 


GET01200 


IF1C2.LT. Cl .OR. V2.LT.V1) GO TO 4000 


GET01210 


VI = V2 


GET01220 


Cl = C2 


GET01230 


I F ( FLAG ) GO TO 70 


GET01240 


READ(7,40,END=70) V2 , C2 


GET01250 


GO TO 50 


GET01260 

60 

CUR = Cl + (C2 - Cl) / ( V2 - VI) * (VOLT - VI) 


GET01270 


GO TO 100 


GET01280 

70 

RTST = V2 / C2 


GET01290 


I F ( V2 . LT . 2 .0*VGAP * FLOAT (NOOUN ) ) GO TO 6000 


GET01300 


CUR = VOLT /RTST 


GET01310 


FLAG = .TRUE. 


GET01320 

100 

CONTINUE 


GET01330 

C 



GET01340 

1000 

RETURN 


GET01350 

C 



GET01360 

2000 

WR I TE ( 6 , 2001 ) 


GET01370 

2001 

FORMAT( IX, ’ERROR DETECTED IN SUBROUTINE EVEN’/ 


GET01380 


llX.’NO DATA FOUND IN IV FILE’) 


GET01390 


STOP 


GET01400 

3000 

WR ITE ( 6 , 3001 ) 


GET01410 

3001 

FORMAT! IX, ’ERROR DETECTED IN SUBROUTINE EVEN’/ 


GET01420 


11X, ’FIRST POINT IN IV EXPERIMENTAL DATA SHOULD BE 

0.0, 0.0’ ) 

GET01430 


STOP 


GET01440 

4000 

WR I TE ( 6 , 4001 ) VI, V 2, Cl, C2 


GET01450 

4001 

F ORMAT ( IX , ’ERROR DETECTED IN SUBROUTINE EVEN’/ 


GET01460 


UX, ’POINTS ARE NOT MONOTON I CAL LY INCREASING’/ 


GET01470 


21X, ’VOLTAGE 1 = ’ , E 1 5 . 8 , 5X , ’ VOLTAGE 2 = ’.E15.8/ 


GET01480 


31X, ’CURRENT 1 = ’ , E 1 5 . 8 , 5X , ’ CURRENT 2 = ’,E15.8> 


GET01490 


STOP 


GET01500 

5000 

WRITE! 6, 5001 ) 


GET01510 

5001 

FORMAT! IX, ’ERROR DETECTED IN SUBROUTINE EVEN’/ 


GET01520 


1 IX, ’RANGE OF GIVEN I-V DATA POINTS DOES NOT COVER 

REQD SPACE’ ! 

) GET01530 


STOP 


GET01540 

6000 

WRITE (6, 6001 ) 


GET01550 

6001 

FORMAT! IX, ’ERROR DETECTED IN SUBROUTINE EVEN’/ 


GET01560 


1 1 X , ’INPUT I-V CURVE DATA DOES NOT EXTEND FAR ENOUGH »/ 

GET01570 


21X , ’ INTO THE LINEAR RANGE’) 


GET01580 


STOP 


GETf 1590 


END 


GET01600 


SUBROUTINE REF I NE ! VOL , CUR , NUM , RBEFOR , RAFTER ) 


GET01610 


DIMENSION VOL ( NUM ) , CUR(NUM) 


GET01620 

C 



GET01630 

C 

SUBROUTINE REFINE REMOVES THE OHMIC PORTION FROM 

THE I-V DATA 

SUCHGET01 640 

C 

THAT THE LAST DATA POINT HAS ZERO CURRENT 


GET01650 
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FILEs GETCHI FORTRAN A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


C 

C WRITEIG.90) NUM 

C90 F0RMATI/1X, ’SUBROUTINE REFINE: REMOVING THE OHMIC PORTION FROM’ 
C 1/20X,’I-V DATA SUCH THAT THE LAST POINT HAS ZERO CURRENT’/ 

C 220X , ’ NUMEER OF POINTS PROCESSED = ’,15) 

C 

GADD = CUR ( NUM > / VOLiNUM) 

C 

DO 40 I = 1 , NUM 

40 CUR ( I ) = CUR C I > - GADD * VOL ( I ) 

C 

GBEFOR = l./RBEFOR 
GAFTEP. = GBEFOR + GADD 
RAFTER = 1 . / GAFTER 
C 

RETURN 

C 

END 

SUBROUTINE I NEXP C T I MI N , CH I , NUM , TOUT , CHOUT , NTOTAL , NOPNTS > 

REAL TOUT ( NTOTAL ) , CHOUT ( NTOTAL > 

DIMENSION TIMIN(NUM), CHI(NUM) 

LOGICAL FLAG 
COMMON /GAP INF/ VGAP 

DATA PI, HBAR , E/3 . 1 4 1 592S , 1 . 055E-34 , 1 . S02E- 1 9/ 

DATA TOL 1 , TOL2/1 .E-7, 1 .E-3/ 

C 

C INTERPOLATION PROGRAM 

C 

C SUBROUTINE INEXP INTERPOLATES THE I-V DATA SUCH THAT 

C THE SPACING AND NUMBER OF POINTS CORRESPOND TO THE 
C CORRECT SPACING FOR THE CHI CURVE 

C 

C THE INTERPOLATION FORMULA USED IS A + B*X + C*EXP(D*X) 

C EQUAL SPACED POINTS MUST BE USED 

C 

FLAG = .FALSE. 

C 

C WR ITE ( S , 700) 

C700 FORMATI/1X, ’SUBROUTINE INEXP: INTERPOLATES THE I-V DATA SUCH’, 

C 1 IX , ’ THAT THE SPAC I NG ’ / 1 9X , ’ AND NUMBER OF POINTS CORRESPOND’/ 

C 2 1 9X , ’TO THE CORRECT SPACING FOR THE CHI CURVE’/) 

C 

WR ITE ( G , 720) 

720 FCRMATI/1X, ’ INPUT THE LO FREQUENCY IN HERTZ < FREQ) (E15.8)’) 

READ { 5,740) FPEQ 
740 FORMAT { E 1 5 . 8 ) 

DELTA = E * VGAP / 2.0 

N 1 = 2.0 * DELTA / (0. 064*FREQ*2 .0* P I *HBAR ) 

N2 = 260 / N 1 

WR I TE ( 6 , 7E0 ) N2 , N1 

750 F0RMAT(/1X, ’NUMBER OF LO CYCLES OF CHI IN SISCAP SHOULD BE’, 

iix ’ > ’ is/ 

2 1 X ] ’THE NUMBER OF POINTS PER LO CYCLE IN SISCAP ’, 

3’SHOULD BE > ’ , 155 
WR ITE ( 6 , 755 ) 


GET01GG0 

GET01670 

GET01G80 

GET01G90 

GET01700 

GET01710 

GET01720 

GET01730 

GET01740 

GET01750 

GET01760 

GET01770 

GET01780 

GET01790 

GET0180Q 

GET01810 

GET01820 

GET01830 

GET0184Q 

GET01850 

GET01860 

GET01870 

GET01880 

GET01890 

GET01900 

GET01910 

GET0192P 

GET01930 

GET01940 

GET01950 

GET01960 

GET01970 

GET01980 

GET01990 

GET02000 

GET02010 

GET02020 

GET02030 

GET02040 

GET02050 

GET02060 

GET02070 

GET02080 

GET02090 

GET02100 

GET02110 

GET02120 

GET02130 

GET02140 

GET02150 

GET02160 

GET02170 

GET02180 

GET02190 

GET* 2200 
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FILE: 

GETCHI FORTRAN A NASA/GISS CONVERSATIONAL MONITOR 

SYSTEM 

SP302 

755 

FORMATI/1X, ’INPUT NUMBER OF LO CYCLES (NOCYC) <I5)’/1X, 

’INPUT 

NO’ .GET02210 


1* OF POINTS PER LO CYCLE (NOPNTS) (15)’) 


GET02220 


READ (5,760) NOCYC, NOPNTS 


GET02230 

760 

FORMAT ( I 5 > 


GET02240 

C 



GET02250 


DELTAT = 1. / FREO / F LOAT ( NOP NTS ) 


GET02260 


DELTAW = 1 .0/FLOAT!NOPNTS)/FLOAT!NOCYC)/DELTAT*2.0*PI 


GET02270 


DELTAV = HBAR * DELTAW / E / 2.0 


GET02280 


NTOTAL = NOCYC * NOPNTS 


GET02290 


VTOTAL = FLOAT( NTOTAL ) * DELTAV 


GET02300 


WR I TE ( 6 , 55 4 ) NTOTAL 


GET023 10 

554 

FORMAT! IX, ’NOCHI = ’,15) 


GET02320 

C 



GET02330 

C 

WR ITE(6,555> DELTAV, VTOTAL 


GET02340 

C555 

FORMAT! IX, ’DELTAV = ’.E15.8,’ VOLTS’ , IX, ’TOTAL VOLTAGE 

REQD ■= 

‘ GET02350 

C 

1 E 1 5 . 8 , ’ VOLTS’) 


GET02360 

c 



GET02370 


TIME = 0. 


GET02380 


IPNT = 2 


GET02390 


VAL = 0. 


GET02400 


TOUT! 1 > = TIME 


GET024 10 


CHOUT(l) = VAL 


GET02420 

c 



GET02430 


NLIMIT = NTOTAL - I 


GET02440 


DO 100 I = 1, NLIMIT 


GET02450 


TIME = TIKE + DELTAV 


GET02460 


I F ( FLAG ) GO TO 410 


GET02470 

60 

IF (TIM IN! IPNT+1 ) . GT . TIME ) GO TO 70 


GET02480 


IPNT = IPNT + 1 


GET02490 


IF! IPNT+2.GT.NUM) GO TO 410 


GET02500 


GO TO 60 


GET025 10 

C 



GET02520 

70 

YM1 = CHI! I P NT- 1 ) 


GET02530 


Y 0 = CHI ( IPNT) 


GET02540 


Y1 = CHI ( IPNT+1 ) 


GET02550 


Y2 = CHI ( IPNT+2) 


GET02560 


CURV1 = Y2 - 2 . 0*Y 1 + Y 0 


GET02570 


CURV2 = Y1 - 2.0*Y0 + YM1 


GET02580 


IF (CURV1*CURV2 .LT. 0.0) GO TO 470 


GET02590 


I F ! ABS ( CURV1 ) . LE . TOL 1 .OR. ABS ! CUR V2 ) . LE .TOL 1 )GOTO470 


GET02600 


D = AL0G(CURV1/CURV2> 


GET02610 


IF !ABS(D > .LT.TOL2 > GO TO 470 


GET02620 


C = CURV1 / ( ( EXP C D ) — 1 . 0) **2 ) 


GET02630 


A = Y0 - C 


GET02640 


B = Y1 - A - C*EXP ( D > 


GET02650 


X = (TIME-TIMIN! IPNT) )/(TIMIN! IPNT+1 >-TIMIN( IPNT) ) 


GET02660 


VAL = A + B *X + C* EXP ( D*X ) 


GET02670 


GO TO 400 


GET02680 

C 



GET02690 

470 

VAL = (TIME-TIMIN! IPNT) )/(TIMIN! IPNT+1 )-TIMIN( IPNT) ) 


GET02700 


1 * (CHI ( IPNT+1 ) - CHI ( IPNT) > + CHI ( IPNT) 


GET02710 

C 



GET02720 


GO TO 400 


GET02730 

410 

FLAG = .TRUE. 


GET02740 


VAL = 0.0 


GET02750 
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FILEt 


C 

400 

C 

100 

c 

c 


210 


C 


C 

c 

c 

c 

c 

C300 

c 


777 


799 

778 

223 

224 

225 

228 

C 

C 


GETCHI FORTRAN A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


TOUTIl+l) = TIME 
CHOUT { I + 1 > = VAL 

CONTINUE 

I F ( .NOT. FLAG) RETURN 

PERCNT = TIMIN(NUM) / VTOTAL * 100.0 
WRITE!6,210) PERCNT 

FORMAT< IX, ’THE GIVEN POINTS WERE NOT SUFFICIENT TO COVER THE’ 

1, ’ REOD INTERP REG ION ’ / IX , 'THE FUNCTION WAS ASSUMED’, 

2, ’ TO BE 2ER0 FOR THE EXTRAPOLATED PORTION’,/ 

3 1 X , ’THE INTERPOLATED REGION EQUALS ’ ,F8.2,’X OF THE TOTAL RANGE’) 
RETURN 

END 

SUBROUTINE I VTOCH < VOL , CURNEW , NOCH I , DELTAV , T I ME , CH I , CCH I , FACTOR ) 
REAL VOL ! NOCHI ) , CURNEW( NOCH I ) , TIME(NOCHI), CHI(NOCHI) 

REAL LENGTH 

COMPLEX CCH I ( NOCHI ) , FACTOR C NOCH I ) 

THIS SUBROUTINE TRANSFORMS THE DISCRETIZED I-V DATA 
INTO THE EQUIVALENT DISCRETIZED CHI DATA 

WR I TE ( 8 , 300 > 

F0RMATI/1X. ’SUBROUTINE IVTOCH: TRANSFORMING I-V DATA TO CHI DATA 

IBASE = -9 
HEAR = 1 .055E-34 
E = 1.602E-19 


NDIV2 = 

NOCH 

I 




N = 

NDIV2 * 

2 




I F ( N 

.EQ. 

512) 


GO TO 

222 

I F ( N 

.EQ. 

1024 

) 

GO 

TO 

223 

I F ( N 

.EQ. 

2048 

) 

GO 

TO 

224 

I F I N 

.EQ. 

4096 

) 

GO 

TO 

225 

I F ( N 

.EQ. 

8192 

) 

GO 

TO 

226 

WR ITE ( 8 , 

777 ) 


N 



FORMAT! IX, ’N 


DOES 

; NOT El 


1/1X, ’ IT EQUALS ’ , 15) 
WR ITE( 6 , 799 ) 
FORMATUX, ’INPUT THE 
READ 15,778 ) IBASE 
FORMAT! 15 > 


EQUAL 512,1024,2048,4096,8192’ 


BINARY LEOGAR ITHM OF THIS NUMBER’) 


GO 

TO 

222 


IBASE 

= - 

10 

GO 

TO 

222 


IBASE 

= - 

1 1 

GO 

TO 

222 


IBASE 

= - 

12 

GO 

TO 

222 


IBASE 

= - 

13 


SET UP DC VALUE AT ZERO ON I-V CURVE 


GET027G0 
GET02770 
GET02780 
GET02790 
GET02800 
GET028 1 0 
G ETC 2 8 20 
GET02830 
GET£ 2840 
GET02850 
GETC2860 
GET02870 
GET02880 
GET02890 
G ETC 2 900 
GET029 10 
GET02920 
GET02930 
GET02940 
GET02950 
GET02960 
GET0297C 
GET02980 
GETC 2990 
GET03000 
GET03010 
’ GET03020 
GET03030 
GET03040 
GET03050 
GET03060 
GETf 3070 
GET03080 
GET03090 
GET03100 
GET031 10 
GET03120 
GET03130 
GET03140 
GET03150 
GET03160 
GET03170 
GET03180 
GET03190 
GET03200 
GETC 32 10 
GET03220 
GET03230 
GET03240 
GET03250 
GET03260 
GET03270 
GET03280 
GET03290 
GET03300 
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C 

222 TMP = CURNEW! 1 ) 

VOL ( 1 ) = 0. 

CCHI ( 1 > = CMPLX (0. ,0. ) 

C 

DO 100 I = 2.NDIV2 
TMP = CURNEWI I) 

DELTAV = VOL I I > - VHLD 
VHLD = VOL I I ) 

I NEXT = N + 2 - I 

CCHI ( INEXT) = CMPLXI0. ,TMP) 

100 CCH I ( I ) = CMPLX (0. , -TMP ) 

110 F ORMAT (E15.8.E15.8) 

CCHIINDIV2 + 1) = CMPLX(0.,0.) 

C 

C WR I TE I 6 , 4 46 > DELTAV, NDIV2 

C446 FORMAT! 19X, ’DELTAV INPUT IS ’,E15.8,’ VOLTS’/ 

C 1 1 9 X , ’NO OF POINTS INPUT = ’,15) 

DELTAF = DELTAV / HBAR / 2. / 3.1415926 * E 
FRMAX = DELTAF * FLOAT I ND I V2 > 

C 

C NYQUIST SAMPLING IS TWICE THAT OF MAX FREQUENCY IN W FORM 
C 

DELTAT = 1 . / (FRMAX * 2.) 

LENGTH = FRMAX * 2. 

C 

CALL FFTI I BASE , CCH I , FACTOR ) 

C 

DO 200 I = 1.NDIV2 
TIME! I ) = FLOAT ( I - 1 ) * DELTAT 
C 

C THE VALUE OF L ENGTH/F L OAT I N ) EQUAL DELTAF 

C 

CHIII) = 2. * REAL I CCHI I I > ) * LENGTH / FLOATIN) 

200 CONTINUE 
C 

C WRITEI6.310) DELTAT, NDIV2 

C310 FORMAT! 19X, ’DELTAT ON OUTPUT = ’.E15.8,’ SECS’/ 

C 1 1 9X , ’ NO OF POINTS OUTPUT = ’,15) 

C 

RETURN 

END 

SUBROUTINE FFT! NUMBER , DATA , FACTOR > 

INTEGER CRPSZE.POSN.POINTR , ELMTNO , FACPTR .ABSNUM 
INTEGER GPPTR 
LOGICAL INVERT 

COMPLEX TEMP ,W, FACTOR ( NUMBER ) .DATA! NUMBER > .WSCALE, TEMPI , 
1TEMP2 

DATA PI/3.1415926/ 

C 

C AUTHOR: R.G. HICKS, UNIV. OF QLD . , APRIL 1977. 

C 

C AIM: THIS SUBROUTINE PERFORMS A FAST FOURIER TRANSFORM ON THE 

C INPUT DATA USING THE COOL EY-TURKE Y ALGORITHM 

C THERE MUST EE 2**NUMBER POINTS GIVEN 


GET033 10 
GET03320 
GET03330 
GET03340 
GET03350 
GET03360 
GET03370 
GET03380 
GET0339C 
GET0340C 
GET03410 
GET03420 
GET03430 
GET03440 
GET03450 
GET03460 
GET03470 
GET03480 
GET03490 
GET03500 
GET03510 
GET03520 
GET03530 
GET03540 
GETL3550 
GET03560 
GET03570 
GET03580 
GETD3590 
GET03600 
GET03610 
GET03620 
GET03630 
GET03640 
GET03650 
GET0366P 
GET03670 
GET03680 
GET03690 
GET03700 
GET03710 
GET03720 
GET03730 
GET03740 
GET03750 
GET03760 
GET03770 
GETf 3780 
GET03790 
GET03800 
GETL3810 
GETL3820 
GET03830 
GET03840 
GET03850 
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FILE: 

GETCHI FORTRAN A NASA/GISS 

CONVERSATIONAL 

MONITOR SYSTEM 

SP302 

C 

AND THE GIVEN POINTS ARE 

DESTROYED BY THE 

SUBROUTINE. 

GET03860 

c 

IF NUMBER IS NEGATIVE, AN 

INVERSE FOURIER 

TRANSFORM IS 

DONE.GET03870 

c 




GET03880 


INVERT = .FALSE. 



GET03890 


ABSNUM = NUMBER 



GET03900 


I F ( NUMBER . GT .0 ) GO TO 10 



GET039 10 


ABSNUM = -NUMEER 



GET03920 


INVERT = .TRUE. 



GET03930 

10 

N = 2 ** ABSNUM 



GET03940 


ARG = 2. * PI / FLOAT < N ) 



GET03950 


TEMP = CMPLX(0. ,ARG) 



GET03960 


W = CEXP ( TEMP ) 



GET03970 

c 




GET03980 


FACTOR < 1 ) = W 



GET03990 


NDIV2 = N/2 



GET04000 


DO 15 INDEX = 2.NDIV2 



GET04010 

15 

FACTOR! INDEX) = FACTOR! INDEX - 

1 >*W 


GET04020 

c 




GET04030 


IF! INVERT) GO TO 20 



GET04040 


DO 17 POSN = 1 ,N 



GET04050 

17 

DATA! POSN ) = CONOG ( DATA! POSN ) ) 



GET04060 

c 




GET04070 

20 

GRPSZE = N 



GET04080 


NOGRP = 1 



GET04090 


ITERNO = 0 



GET04100 

25 

LIMIT = GRPSZE/2 



GET041 10 


ITERNO = ITERNO + 1 



GET04120 


POINTR = GRPSZE/2 + 1 



GET04130 


POSN = 1 



GET04140 

c 




GET04150 


DO 60 GRPTR = 1 , NOGRP 



GET04160 


DO 50 ELMTNO = 1 .LIMIT 



GET04170 


I POSN = POSN - 1 



GET04180 


CALL FACNO! IPOSN.N, ITERNO, FACPTR) 


GET04 1 90 


I F ( FACPTR . EQ .0) GO TO 30 



GET04200 


WSCALE = FACTOR I FACPTR > 



GETI 4210 


TEMPI = DATA! POINTR ) * WSCALE + 

DATA! POSN) 


GET04220 


WSCALE = -WSCALE 



GET! 4230 


TEMP2 = DATA! POINTR ) * WSCALE + 

DATA! POSN ) 


GET0424C 


GO TO 40 



GET04250 

30 

TEMPI = DATA! POI NTR ) + DATA! POSN > 


GET04260 


TEMP2 = DATA! POSN ) - DATA! POI NTR ) 


GET04270 

40 

DATA! POSN) = TEMPI 



GET04280 


DATA! POI NTR ) = TEMP2 



GET04290 


POSN = POSN + 1 



GET04300 

50 

POINTR = POINTR + 1 



GET04310 

c 




GET04320 


POSN = POSN + GRPSZE/2 



GET04330 

60 

POINTR = POINTR + GRPSZE/2 



GET04340 

c 




GET04350 


GRPSZE = GRPSZE/2 



GET04360 


NOGRP = NOGRP * 2 



GET04370 


I F ( GRPSZE .GT. 1 ) GO TO 25 



GET04380 

c 




GET04390 


DO 70 POSN = 1 , N 



GET04400 
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70 

C 


80 

C 


C 

c 

c 

c 


c 


10 

20 


C 

c 

c 

c 


10 


IPOSN = PCSN - 1 

CALL OUXTAIN, IPOSN, OUXPOS) 

OUXPOS = OUXPOS + 1 
IFIOUXPOS.LT. FOSN) GO TO 70 
TEMP = DATA! POSN) 

DATA(POSN) = DATA! OUXPOS > 

DATA! OUXPOS ) = TEMP 
CONTINUE 

IF(INVERT) RETURN 

DO 80 POSN = 1 ,N 

TEMP = CONOG ( DATA! POSN ) ) 

DATA(POSN) = TEMP / N 

RETURN 

END 

SUBROUTINE OUXTAIN, IPOSN, OUXPOS ) 

THIS SUBROUTINE DETERMINES THE CORRECT MULTIPL 
IN THE BUTTERFLY ALGORITHM OF THE FAST FOURIER 

NN = N 
NUM = IPOSN 
OUXPOS = 0 

DO 10 1=1 ,N 

NN = NN/2 

MEM = NUM 

NUM = NUM/2 

IREM = MEM - 2*NUM 

OUXPOS = 2*OUXPOS + IREM 

IFINN.EQ. 1 ) GO TO 20 

CONTINUE 

RETURN 

END 

SUBROUTINE FACNOI I POSN , N , ITERNO, FACPTR ) 

INTEGER FACPTR 

THIS SUBROUTINE DETERMINES THE CORRECT FACTOR 
CALCULATION OF THE BUTTERFLY ALGORITHM IN THE 


I CATION SEQUENCING 
TRANSFORM 


GET044 10 
GET04420 
GET04430 
GET04440 
GET04450 
GET044G0 
GET04470 
GET04480 
GET04490 
GET04500 
GET04510 
GET04520 
GET04530 
GET04540 
GET04550 
GET04560 
GET04570 
GET04580 
GET04590 
GET£ 4600 
GET04610 
GET04620 
GET04630 
GET046 40 
GET04650 
GET04660 
GET04670 
GET04G80 
GET04690 
GET04700 
GETf 4710 
GET04720 
GET04730 
GET04740 
GET04750 
GET047G0 
GET04770 
GET04780 

FOR CORRECT GET04790 

FAST FOURIER TRANSFOGET04800 

GET04810 


FACPTR = 0 
I D I V = N/ ( 2** ITERNO) 

ITEMP = IPOSN 

ITEMP = ITEMP/IDIV 

MUL = N/2 

DO 10 1=1 , ITERNO 

IREM = ITEMP - ITEMP/2*2 

FACPTR = IREM*MUL + FACPTR 

ITEMP = I TEMP / 2 

MUL = MUL / 2 

CONTINUE 

RETURN 

END 

SUBROUTINE REEXP < TI MI N , CH I , NUM , TNEW, CHNEW, NOI V , DELTAV ) 


GET04820 
GET04830 
GET04840 
GET04850 
GET04860 
GET04870 
GET04880 
GET04890 
GET04900 
GET049 10 
GET04920 
GET04930 
GET04940 
GETC4950 



non noooonoooooo 
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C30 

C 

40 


50 

C 


C 


C 


G0 


C 

70 


REAL TNEW(NOIV), CHNEW(NOIV) 

DIMENSION TIMIN(NUM), CHI(NUM) 

LOGICAL FLAG 

COMMON /GAP INF/ VGAP 

DATA TOL 1 , TOL2/1 .E-7, 1 .E-3/ 

INTERPOLATION PROGRAM 

THIS SUBROUTINE INTERPOLATES THE DATA SUCH THAT THE VECTOR 
J MAY BE CALCULATED WITH REASONABLE ACCURACY. 

TYPICALLY DELTAV SHOULD BE ( VGAP/3 .0E-3 )*2 .0E-5 VOLTS 
SHOULD BE USED. 

THE INTERPOLATION FORMULA USED IS A + B*X + C*EXP(D*X> 

THE POINTS USED MUST BE EQUALLY SPACED 
FLAG = .FALSE. 

VSUG = (VGAP/3. 0E-3) * 0.00002 


WRITE (6, 30) 

FORMAT( /IX , ’SUBROUTINE REEXP: PROVIDES THE OUTPUT IV DATA’, 

1 , ’ FOR PROGRAM SISCAP ’ ) 

WR ITE ( G , 40) VSUG 

F0RMAT(/1X, ’INPUT DELTAV ( TYP ’.Eli. 4,’ VOLTS) (E15.8)’/ 

1 1 X , ’INPUT NUMEER OF POINTS REQUIRED (NOIV) (PWR OF 2) (TYP 1024) 
2’ ( 15) ’ ) 

READ ( 5 , 50) DELTAV, NOIV 
FORMAT ( E 1 5 . 8/ I 5 ) 

NTOTAL = NOIV 

VTOTAL = FLOAT( NTOTAL) * DELTAV 

TIME = 0. 

IPNT = 2 
VAL = 0. 

TNEW( 1 ) = TIME 
CHNEW( 1 ) = VAL 

NLIMIT = NTOTAL - 1 

DO 100 I = 1, NLIMIT 

TIME = TIME + DELTAV 

I F ( FLAG ) GO TO 410 

IF (TIMIN( IPNT+1 ) .GT.TIME) GO TO 70 

IPNT = IPNT + 1 

IF ( IPNT+2.GT.NUM) GO TO 410 

GO TO G0 

YM1 = CHI ( IPNT-1 ) 

Y0 = CHI ( IPNT) 

Y 1 = CHI ( IPNT+1 ) 

Y 2 = CH I ( IPNT+2) 

CURV1 = Y2 - 2 . 0*Y 1 + Y0 
CURV2 = Y1 - 2.0*Y0 + YM1 


GET049G0 
GET04970 
GET04980 
GET04990 
GET05000 
GET05010 
GET05020 
GET05030 
GET05040 
GET05050 
GET05060 
G ETC 5070 
GET05080 
GET05090 
GET05100 
GET051 10 
GET05120 
GET05130 
GET05140 
GET05150 
GET05160 
GET05170 
GET05180 
GET05190 
GET05200 
GET052 1 0 
, GET05220 
GET05230 
GET05240 
GET05250 
GET052G0 
GET05270 
GET05280 
GET05290 
GETC 5300 
GET053 10 
GET05320 
GET05330 
GET05340 
GETC 5350 
GET053G0 
GET05370 
GET05380 
GETC5390 
GETC 5 400 
GETC 5410 
GETC 5420 
GET05430 
GET05440 
GET-j5450 
GETC5460 
GETC5470 
GET0548C 
GETC 5490 
GET05500 



64.10 


FILE: GETCHI FORTRAN A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


IF (CUPV1*CURV2 .LT. 0.0) GO TO 470 

IF (ABSICURV1 ) .LE.T0L1 .OR. ABS (CURV2) .LE.T0L1 )GOTO470 
D = AL0G(CURV1/CURV2) 

IF (ABS(D) .LT.T0L2) GO TO 470 
C = CURV1 / ( (EXP(D)-l .0)**2> 

A = Y0 - C 

B = Y1 - A - C*EXP ( D ) 

X = (TIME-TIMINI IPNT) )/(TIMIN( IPNT+1 )— TIMINt IP NT) ) 

VAL = A + B*X + C*EXP ( D*X ) 

GO TO 400 

VAL = (TIME-TIMIN'! IPNT) >/{TIMIN< IPNT+1 )-TIMIN( IPNT) ) 
1 * (CHKIPNT+1) - CHI ( IPNT) ) + CHI ( IPNT) 

GO TO 400 
FLAG = .TRUE. 

VAL =0.0 

TNEWf 1+1 ) = TIME 
CHNEWt 1+1 > = VAL 
CONTINUE 


I F ( .NOT. FLAG) RETURN 


PERCNT = TIMIN(NUM) / VTOTAL * 100.0 
WR ITE (6,210) PERCNT 

FORMAT! IX. ’THE GIVEN POINTS WERE NOT SUFFICIENT TO COVER THE’ 

1, ’ REOD INTERP REG I ON ’ / IX , ’ THE FUNCTION WAS ASSUMED’, 

2, ’ TO BE ZERO FOR THE EXTRAPOLATED PORTION’,/ 

3 1 X , ’THE INTERPOLATED REGION EQUALS *,F8.2,’X OF THE TOTAL RANGE’) 
RETURN 


100.0 


GET05510 
GET05520 
GET05530 
GET05540 
GET05550 
GET055G0 
GET05570 
GETL5580 
GET0 5590 
GET05S00 
GET05S10 
GET05S20 
GET! 5630 
GET05640 
GET05650 
GET05660 
GET05670 
GET05630 
GET05690 
GET05700 
GET057 10 
GET05720 
GETC5730 
GET05740 
GET05750 
GET05760 
GET05770 
GET0 5780 
GET05790 
GET05800 
GET058 10 
GET05820 
GET05830 
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A. 8 Listing of Program SISCAP 

This section gives the listing of the SIS mixer analysis program 
SISCAP. 


FILE: SISCAP 


FORTRAN A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP30Z 


CUROC, CURPMP . 
PI/3.1415926/ 


RLOSS , RTEMP 


C 

C 

C 

C 


1 
C 

C 

C 

C4001 

C 


COMPLEX ZIMPEDt 128) , VSOURC (128) 

COMPLEX C VAL (128), FACTOR 1 204 8 ) , 0( 10241 , ZSSBU27) 

COMPLEX Y0 , Yl, Y1M, FIDDLE, ZLO, CHOLD 
REAL VOLT ( 128), CHK512), CURU28J, OLDVOLU28) 

REAL I VDAT ( 1024 ) 

LOGICAL FLAG 

COMMON /PER I P H/ IN, IOUT, IPLOT 
COMMON /ERRALL/ TOL , EPS 
COMMON /RBOUT/ VDC , VLO, VODC, VOLO 
DATA E/1 .602E-19/, HBAR/1 .055E-34/, 

IN = 5 
IOUT = 6 
IPLOT = 20 

OBTAIN THE VARIOUS INPUT PARAMETERS, BOTH JUNCTION AND ALGORITHM 
PARAMETERS 

READ (11,1) IS I G , NODAT, DELV , VSTART, VLOSTR , DELVLO , NLO , TOL, 

1 EPS, NOPNTS , I BASE 1 , NOIV, IBASE2, NOCHI, NOPER, MMAX , NMAX , 

2 RN, CAP, TEMP, FREQ, FIF, P, RID, Z0 

FORMAT ( I5/I5,4(/E15.8) ,/I5,2(/E15.8) ,8(/I5),8(/E15.8) ) 

NOCOL = 4 


WRITE! IPLOT, 4001 ) 
F ORMAT (215) 


NODAT, NOCOL 


700 


WRITE! 6,700) FREQ, FIF, CAP, RN, TEMP , I BASE2 , NOI V , I BASE 1 , NOP NTS , MMAX , 
1 NMAX , ISIG, NOPER, NOCHI, NODAT, TOL, EPS, DELV, VSTART, 

2VLOSTR , DELVLO, NLO, P, RID, Z0 
FORMAT! /IX, ’FREQ = ’ , E 1 5 . 8 , ’ HZ ’ , 1 X , ’ F I F = ’ , E 1 5 . 8 , ’ HZ ’ / 


SIS00010 
SIS00020 
SIS00030 
SIS00040 
SIS00050 
SIS00060 
SIS00070 
SIS00080 
SIS00090 
SIS00100 
SIS001 10 
SIS00120 
SIS00130 
SIS00140 
SIS00150 
SIS00160 
SIS00170 
SIS00180 
SIS00190 
SIS00200 
SIS00210 
SIS00220 
S IS00230 
SIS00240 
SIS00250 
SIS00260 
SIS00270 
SIS00280 
S IS00290 
SIS00300 
SIS00310 
SIS00320 


1 IX, 

’CAP = ’ , E 1 

5.8, ’F ’, IX, ’RN = ’ ,E15.8, ’OHMS’/ 

SIS00330 

2 1 X , 

’TEMP = ’ , E 1 

5.8, ’K ’/ 

SIS00340 

3 1 X , 

’ I BASE2 = ’, 

15, IX, ’NOIV = ’,15/ 

SIS00350 

4 1 X , 

’ I EASE 1 = ’, 

15, IX, ’NOPNTS = ’ , 15/ 

SIS00360 

5 1 X , 

'MMAX = ’, 

15. IX, ’NMAX = ’,15/ 

SIS00370 

6 1 X , 

'ISIG = ’ , 

15, IX, ’NOPER = ’,15/ 

S I S00380 

7 1 X , 

’NOCHI = ’, 

15, IX, ’NODAT = ’,15/ 

S IS00390 

8 1 X , 

’TOL = ’ , E 1 5 

,8,9X, ’EPS = ’ ,E15.8/ 

SIS00400 

9 1 X , 

’DELV = ’ , E 1 

5.8, ’VOLTS’ ,3X, ’VSTART = ’ , E 1 5 . 8 , ’ VOLTS ’ / 

SIS00410 

1 IX , 

’VLOSTR = ’, 

E15. 8, ’VOLTS’ , IX, ’DELVLO = ’ , E 1 5 . 8 , ’ VOLTS ’ / 

SIS00420 

1 1 X , 

’NLO = ’ , 15, 

1 9X , ’ P = ’ , E 1 5 . 8 / 

SIS00430 

2 1 X , 

’RID = ’ ,E15 

.8, ’OHMS’ ,5X, ’Z0 = ’ ,E15.8, ’OHMS’ ) 

SIS00440 

SIS00450 

CHECK THE VALID 

ITY OF CONTROL PARAMETERS 

SIS00460 




S IS00470 

IF(FREQ.LE.FIF) 

GO TO 8002 

SIS00480 

IF(FIF.LE.0 S3) 

GO TO 8002 

SIS00490 

I F ! 2** I BASE 2 . NE 

. NC I V ) GO TO 8002 

SIS00500 


IF(2**IBASE1.NE. NOPNTS) GO TO 8002 
I F < MMAX .GT.NOPNTS/2-1 > GO TO 8002 
IF (N MAX. GT. NOP NTS/2-1 > GO TO 8002 
IF(NOCHI.LT. NOPNTS) GO TO 8002 
IF! ISIG.LE.0.OR. ISIG.GT.3) GO TO 8002 


SIS00510 
S IS00520 
SIS00530 
SIS00540 
S IS00550 
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c 



S IS00560 

c 

OBTAIN THE CHI FUNCTION BY READING IT FROM DEVICE 4 

SIS00570 

c 



SIS00580 


REWIND 4 


SIS00590 


WR ITE ( 6 , 601 ) 


SIS00600 

601 

F0RMATI/1X, ’INPUTTING THE CHI FUNCTION’) 


SIS00610 


READ ( 4 , 50 ) NFINAL 


S IS00620 

50 

FORMAT (15) 


S IS00630 


IF ( NF INAL. NE. NOCHI) GO TO 8000 


SIS00640 


DUMHLD = 0.0 


SIS00650 


DO 70 I = 1 .NFINAL 


SIS00660 


READ I 4 , 80) TMP , CHIII) 


SIS00670 

80 

F ORMAT ( 2E 1 5 . 8 ) 


SIS00680 


DELTAT = TMP - DUMHLD 


SIS00690 


DUMHLD = TMP 


SIS00700 

70 

CONTINUE 


SIS00710 


WR I TE 1 6 , 602 ) DELTAT, NFINAL 


SIS00720 

602 

FORMAT! IX, ’DELTAT = ’.E15.8,’ SECS’/ 


SIS00730 


I IX , ’ NUMBEP OF POINTS READ = ’,15) 


SIS00740 

C 



SIS00750 

C 

OBTAIN THE SIS I-V CHARACTERISTIC BY READING IT FROM DEVICE 8 

SIS00760 

c 



S IS00770 


REWIND 8 


SIS00780 

498 

WR I TE < 6 , 50 1 > 


SIS00790 

501 

F0RMATI/1X ,’ INPUTTING THE SIS IV CHARACTERISTIC’) 

SIS00800 


READ ( 8 , 100 ) NDUM 


SIS00810 

100 

FORMAT < I 5 > 


SIS00820 


IF(NDUM.NE.NOIV) GO TO 7000 


SIS00830 


DUMHLD = 0.0 


SIS00840 


DO 110 II = 1 , NDUM 


S IS00850 


READ (8,1 20 > DUMMY, I VDAT (II) 


SIS00860 

120 

FORMAT { 2E 1 5 . 8 ) 


SIS00870 


DELTAV = DUMMY - DUMHLD 


SIS00880 


DUMHLD = DUMMY 


SIS00890 

110 

CONTINUE 


SIS00900 


WR I TE ( 6 , 507 ) DELTAV, NDUM 


SIS00910 

507 

FORMAT( IX, 'DELTAV = ’.E15.8,’ VOLTS’/ 


SIS00920 


11X, ’NUMBEP OF POINTS READ = ’,15) 


SIS00930 

C 



SIS00940 

c 

CALCULATE THE VECTOR J 


SIS00950 

c 



SIS00960 


WRITE (6, 502) 


SIS00970 

502 

F0RMAT(/1X, ’CALCULATING THE VECTOR O’) 


SIS00980 


CALL GETO ( I VDAT . NO I V . FREQ , I BASE2 , 0 , FACTOR , RN , DELTAV ) 

S IS00990 

C 



SIS01000 

C 

INITIALIZE TO ZERO THE JUNCTION VOLTAGE 

VECTOR 

SIS01010 

c 



S IS01020 


DO 10 1 = 1 , NOP NTS 


SIS01030 


OLDVOL ( I ) = 0. 


SIS01040 

10 

VOLT ( I ) = 0. 


SIS01050 

c 



SIS01060 

c 

INPUT THE HARMONIC IMPEDANCES FROM UNIT 

NO. 15 

SIS01070 

c 



SIS01080 


WR I TE ( 6 , 603 ) 


S IS01090 

603 

F0RMATI/1X, ’ INPUTTING THE LO HARMONIC IMPEDANCES’) 

SIS01 100 
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C 

C 

c 

499 


CALL GETZ (Z IMPED, NOP NTS, FREQ, CAP .VSOURC) 

ZLO = Z I MPED ( 2 ) 

GET THE SIDEBAND IMPEDANCES FROM UNIT 14 
WR ITE ( 6 , 499 ) 

F0RMAT(/1X, • INPUTTING SIDEBAND TERMINATION INFORMATION’) 
CALL GETSSB ( ZSSB, NOP NTS, FREQ, FIF , CAP ) 

Y0 = CMP LX ( 1 .0,0.0) / ZSSB ( 1 > 

Y 1 = CMPLXI 1 .0,0.0) / ZSSB ( 2 ) 

Y1M = CMPLXC 1 .0,0.0) / CONOG ( ZSSB ( 3 ) ) 

VLO = VLOSTR - DEL VLO 
DO 5000 LIM = 1 ,NLO 
VLO = VLO + DELVLO 


VDC = VSTART 


DEL V 


RUN THE LARGE SIGNAL ANALYSIS PROGRAM ’NODAT’ TIMES 
FOR EACH DC VOLTAGE BIAS POINT TO BE CONSIDERED 


DO 5000 KIM = 1, NODAT 
KKIMM = KIM + (LIM - 1) * NODAT 
WRITE (6,6000) KKIMM 

6000 FORMAT! /IX , <***********************•*********'/ 
1 IX , ’ ******** * *CASE ’ , 13, ’***************'/ 

2 ix >*********************************’/) 

VDC = VDC + DELV 


C 

607 

C 


C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

9000 


WR 1 TE ( 6 , 607 ) VDC, VLO 

F0RMAT(/1X, ’EMBEDDING THEVENIN VDC SOURCE ’, 8X ,’ = ’,E15.8,’ 
11X, ’EMBEDDING THEVENIN VLO SOURCE (PEAK) = ’,E15.8,’ VOLTS' 

VJ DC = 0.0 
VJLO = 0.0 
CURDC = 0.0 
CURPMP = 0.0 
R LOSS = 0.0 
RTEMP = 0.0 

DETERMINE WHICH METHOD OF ANALYSIS IS REQUIRED 

1 = THREE FREQUENCY MODEL ASSUMPTION 

2 = VOLTAGE UPDATE GENERAL ANALYSIS 

3 = MULTIPLE REFLECTION GENERAL ANALYSIS TECHNIQUE 

GO TO (9000,4000,4000), ISIG 

FOR THE 3-PORT ANALYSIS METHOD 

VODC = VDC 
VOLO = VLO 

ARG = 2.0 * PI / FLOAT ( NOP NTS ) 

DO 7777 MNZ = l.NOPNTS 


IS01U0 
IS01 120 
IS01130 
IS01140 
IS01 150 
IS01160 
IS01 170 
IS01 180 
IS01 190 
IS01200 
IS01210 
IS01220 
IS01230 
IS01240 
IS01250 
IS01260 
IS01270 
IS01280 
IS01290 
IS01300 
1S01310 
IS01320 
IS01330 
1S01340 
IS01350 
IS01360 
IS01370 
IS01380 
IS01390 
IS01400 
VOLTS ’/SIS01410 
IS01420 
IS01430 
IS01440 
IS01450 
IS01460 
IS01470 
IS01480 
IS01490 
IS01500 
IS01510 
IS0 1520 
IS0 1530 
IS01540 
IS01550 
I SC 1 560 
IS01570 
IS01580 
ISC 1590 
ISf 1600 
ISf 1610 
ISC 1 620 
IS01630 
IS01640 
XS01650 


SI 

SI 

SI 

SI 

S) 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 


SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

SI 

S] 

SI 

SI 

SI 

SI 

SI 

SI 

SI 
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SYSTEM SP 302 


PHASE = ARG * FLOATIMNZ -1) 


SIS01660 

7777 

VOLT(MNZ) = VDC - VLO * COS(PHASE) 


SIS01670 

C 



SIS01680 


WRITE <6, 8888) 


SIS01690 

8888 

F0RMAT(/1X, ’PERFORMING A 3-PORT ANALYSIS METHOD’) 


SIS01700 


GO TO 3000 


SIS01710 

C 



SIS01720 

C 

FOR THE GENERAL ANALYSIS METHOD 


S I S£ 1730 

C 



SISjl 1740 

4000 

VSOURC ( 1 ) = VDC 


SIS01750 


V50URC { 2 ) = VLO/2 . 


SIS01760 

C 



SIS01770 

C 

BCAP = 2. * PI * FREQ * CAP 


SIS01780 

c 

FIDDLE = ZLO * CMP LX ( 0 . , BCAP ) + CMPLX ( 1 . 0 , 0 . 0 ) 


SIS01790 

c 

VSOURC ( 2 ) = VSOURC < 2 ) / FIDDLE 


SIS01800 

c 



SIS018U 


RSOURC = REAL ( ZLO ) 


SIS01820 


PLO = ( CABS (VSOURC(2))**2) / 2. / RSOURC 


SIS01830 


WR I TE ( 6 , 4002 ) PLO 


S I S01 840 

4002 

F0RMAT(/1X, ’INCIDENT LOCAL OSCILLATOR POWER = ’.E15.8,’ 

WATTS ’ / ) 

SIS01850 

c 



SIS018 6X' 


IF( ISIG.E0.2) GO TO 6003 


SIS01870 

c 



SIS01880 

c 

FOR THE MULTIPLE REFLECTION TECHNIQUE 


SIS01890 

c 



SIS01900 


CALL MRT (VOLT, CUR, NOP NTS, I BASE 1 , CHI, NOCHI ,RN, 


SIS01910 


1 FACTOR .FREQ, Z IMP ED, VSOURC, OLDVOL .FLAG, Z0, NOP ER ) 


SIS01920 

c 



S ISC 1 930 


I F ( FLAG ) GO TO 5001 


S ISC 1 940 

c 



SIS01950 


DO 4321 MGM = 1 ,NOPNTS 


SIS0196r 

4321 

WR ITE ( 6 , 4322 ) MGM. VOLT(MGM) , CUR ( MGM ) 


SIS01970 

4322 

FORMAT! IX. ’POINT NO. = ’ , I 5 , ’ VOLTAGE = ’ , E 1 5 . 8 , ’ CURRENT 

= ’ , E 1 5 . 

8 > S I S01 980 

C 



SISC1990 


GO TO 6001 


SIS02000 

C 



SIS02010 

C 

FOR THE VOLTAGE UPDATE METHOD 


S I SC 2020 

c 



S IS02030 

6003 

WR ITE ( 6 , 604 ) 


S I SC 2040 

604 

F0RMAT(/1X, ’GENERAL NONLINEAR ANALYSIS USING VOLTAGE UPDATE’, 

SIS02050 


1’ METHOD’/IX, ’BEGINNING THE ITERATION PHASE’) 


SIS02060 

C 



SIS02070 


DO 1000 I = 1 , NOP ER 


SIS02080 

C 



SIS02090 


WRITE (6, 400) I 


SIS02100 

400 

FORMAT! /8X ,’ ITERATION NUMBER = ’,15) 


SIS02110 

C 



SIS02120 


CALL S I SNL! VOLT, CUR, NOP NTS, I BASE 1 , CHI , NOCH I , RN , FACTOR , FR EQ , R I D ) 

SIS02130 


CALL L I NEAR (VOLT, CUR, NOP NTS, I BASE 1 .FACTOR, Z IMPED, VSOURC 

, R I D ) 

SIS02140 


CALL COMP AR( VOLT, CL DVOL , NOP NTS .FLAG) 


SIS0215C 

C 



SIS02160 


DO 500 00 = 1 , NOP NTS 


SIS02170 

500 

OLDVOL ( 00 ) = P * VOLT ! 00 > + (1.0 - P) * OLDVOL(OO) 


SIS02180 

C 



SIS02190 


IF!. NOT. FLAG) GO TO 2000 


S IS02200 
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C 

1000 

C 

60S 

C 

2000 

2002 

C 

C 

C 

c 

G002 

9002 

C 

4323 

C 

C 

C 

C 

6001 

2010 


CONTINUE 

WR ITE( 6 , 606 ) NOPER 

FORMATI/8X, ’LARGE SIGNAL ANALYSIS DID NOT CONVERGE IN ’,15,’ITNS’ 
GO TO 5001 


,15,’ ITNS ’ / ) 


C 

C 

C 

C 

C 

3000 

C 

C5001 

C4500 

C 

5001 

C 

5000 

C 

c 

7000 

1001 

C 

8000 


WRITE ( 6 , 2002) I 

FORMATI/1X, ’LARGE SIGNAL ANALYSIS COMPLETED IN 

CORRECT FOR THE IDENTITY ELEMENT IN THE VOLTAGE 
UPDATE METHOD 

DO 9002 IM = 1 .NOPNTS 

CUR ( I M > = CUR(IM) + VOLT ( I M ) / RID 

DO 4323 MGM = 1 .NOPNTS 

WR ITE ( 6 , 4322 ) MGM, VOLT(MGM), CUR I MGM ) 

DETERMINE THE NONLINEAR OPERATING POINT PARAMETERS 
FROM THE RESULTS OF THE GENERAL ANALYSIS 

DO 2010 IM = 1, NOPNTS 
CVALUM) = VOLT(IM) 

CONTINUE 

CALL FFT! IBASE1 , CVAL, FACTOR) 

VUDC = REAL ( CVAL ( 1 ) ) 

VJLO = CAES ! CVAL < 2 ) ) * 2.0 

VB IAS = VO DC 

VNORM = E * VBIAS / HBAR 

CALL EVALO! VNORM, CHOL D , 0 , NOI V, DELTAV, RN ) 

CURDC = AIMAG(CHOLD) 

SMALL SIGNAL ANALYSIS SECTION 

DO THE SMALL SIGNAL CONVERSION AND NOISE ANALYSIS 

CALL TOPS ( FREQ, F IF, I BASE2, NOI V, I BASE 1 .NOPNTS, MMAX, CAP, RN, TEMP, 
1 NMAX, VOLT, Y0, Y1 ,Y1M,ZL0, I VDAT.O .FACTOR, DELTAV, IS IG.ZSSB) 

WRITE! I PLOT , 4500) VODC, CURDC, CURPMP, RLOSS, RTEMP 
F ORMAT < 5E 1 5 . 8 > 

CONTINUE 

CONTINUE 

STOP 

WRITE ( 6 , 1001 ) 

FORMAT! IX, ’PROBLEM WITH NUMBER OF POINTS IN IVDAT FILE’) 

STOP 

WR ITE I 6 , 8001 ) 


SIS02210 
SIS02220 
S IS02230 
SIS02240 
>S I S02250 
SIS02260 
SIS02270 
SIS02280 
SIS02290 
SIS02300 
SIS02310 
SIS02320 
S IS0233B 
SIS02340 
SIS02350 
SIS02360 
S IS02370 
SIS02380 
SIS02390 
S IS02400 
SIS02410 
SIS02420 
SIS02430 
S IS02440 
S IS02450 
SIS02460 
SIS02470 
SIS02480 
SIS02490 
SIS02500 
SIS02510 
SIS02520 
SIS02530 
SIS02540 
SIS02550 
S IS02560 
S IS02570 
SIS0258C 
S IS02590 
S IS02600 
SIS02610 
S IS02620 
SIS02630 
S IS02640 
S IS02650 
SIS02660 
S I SD 2670 
SIS02680 
SIS02690 
S IS02700 
SIS02710 
SIS02720 
SIS02730 
SIS02740 
SIS02750 
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8001 FORMAT ( IX, ’PROBLEM WITH NUMBER OF POINTS IN NOCHI FILE’) 

STOP 

C 

8002 WR I TE { S , 8003 ) 

8003 F OR MAT ( IX, ’ERROR DETECTED IN CONTROL PARAMETERS’) 

STOP 

C 

END 

SUBROUTINE GETSSB ! ZSSB , NOP NTS , FREQ , F I F , CAP ) 

COMPLEX ZSSB { NOPNTS ) , YCAP 
COMPLEX YSIG, YIM 
DATA PI/3. 1415926/ 

DATA NOPRT/ 10/ 

C 

C READ IN THE SIDEBAND IMPEDANCES FROM UNIT NUMBER 

C 14. THE ORDER OF THE IMPEDANCE DATA IS IF, 

C USB NO . 1 , LSB NO . 1 , USB NO. 2, LSB NO. 2, AND SO ON 

C 

READ (14,1101) NCHK 
1101 FORMAT (15) 

C IF { NChK.NE.NOPNTS-1 ) GO TO 8000 

NOMINI = NOPNTS - 1 
READ ( 14,1120) ZSSB ( 1 > 

1120 F ORMAT ( 2E 1 5 . 8 ) 

WR I TE ( 6 , 300 ) NOPRT, ZSSB(l) 

300 FORMAT! IX, ’THE FIRST ’,13,’ SIDEBAND IMPEDANCES !EXC CAPAC ) s ’, 

1/1X,’ 0 ’ .E15.8, ’+0 ’ ,E15.8, ’OHMS’ ) 

DO 801 II = 2.NOMIN1 
ZSSB (II) = CMFLX(0. 0,0.0) 

ISDBD =11/2 

IF! ISDBD*2.NE. I I .AND. II.NE.l) ISDBD = -ISDBD 
IF! II .GT. NCHK) GO TO 2000 
READ (14,11 20 > ZSSB! II ) 

2000 IF! I I .LE. NOPRT) WRITE<6,301) ISDBD, ZSSB(II) 

301 FORMAT! 1X,I5,1X,E15.8,’+0’,E15.8, ’OHMS’ ) 

IF! 1 1/2*2. NE. I I ) ZSSB (II) = CONUG ( ZSSE (II)) 

801 CONTINUE 
C 

C ADD IN THE JUNCTION PARASITIC CAPACITANCE CONTRIBUTION 

C 

CON = 2. * PI * CAP 
DO 398 II = l.NCHF 

I F ! CAES ( ZSSB (II)) .LT. l.E-10) GO TO 398 
M = II / 2 

I F ( M . EQ . 1 .OR. M.EO.0) GO TO 398 
IF (M*2 .NE. I I .AND. I I .NE. 1 ) M = -M 
BCAP = ! FLOAT ( M ) *FREQ + FIF) * CON 
YCAP = CMP LX ( 0 . , BCAP > 

ZSSB ! I I ) = 1 . /! 1 . /ZSSB! I I ) + YCAP) 

398 CONTINUE 
C 

YSIG = CMPLX! 1 .0,0.0) / ZSSB ! 2 ) 

YIM = CMPLX! 1 ,0,0.0) / ZSSB ( 3 ) 

WR I TE ( G , 7000 ) YSIG, YIM 

7000 F ORMAT ( 1 X , ’ Y SIGNAL (INCL CAPAC) = ’,E15.8,’ +J ’ , E 1 5 . 8 , ’ S ’ / 1 X , 


SIS02760 
SIS02770 
S IS02780 
S IS02790 
SIS02800 
SIS02810 
S IS02820 
SIS02830 
S IS02840 
SIS02850 
S IS028S0 
SIS02870 
SIS02880 
SIS02890 
S IS02900 
SIS02910 
SIS0292B 
SIS02930 
SIS02940 
S IS02950 
SIS029B0 
SIS02970 
SIS02980 
SIS02990 
SIS03000 
SIS03010 
SIS03020 
SIS03030 
SIS03040 
SIS03050 
SIS03060 
SIS03070 
SIS03080 
SIS03090 
SIS03100 
SISB31 10 
SIS03120 
SIS03130 
SIS03140 
SIS03150 
SIS03160 
SIS03170 
SIS03180 
SIS03190 
SIS03200 
SIS03210 
SIS03220 
S IS03230 
SIS03240 
SIS03250 
SIS032G0 
SIS03270 
S IS03280 
SIS03290 
SIS03300 
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FILE: 

SISCAP FORTRAN A NASA/GISS CONVERSATIONAL 

MONITOR SYSTEM 

SP 302 


l’Y IMAGE ( I NCL CAPAC) = ’.E15.8,’ +J ’.E15.8, 

’S’ ) 


SIS03310 

c 




SIS03320 


WRITE ( 6 , 200 ) NCHK 



SIS03330 

200 

FORMAT ( IX, ’NUMBER OF POINTS READ = ’,15) 



S I SO 3340 


RETURN 



S IS03350 

C 




SIS03360 

C8000 

WR ITE ( 6 , 8001 ) 



S I SO 3370 

C8001 

FORMAT! IX. 'PROBEM WITH NUMBER OF POINTS IN ZSSB FILE’) 

S IS03380 

C 

STOP 



SIS03390 


END 



SIS03400 


SUBROUTINE GETZ ! Z I MPED , NOP NTS , FREQ , CAP , VSOURC ) 



SISP341P 


COMPLEX ZIMPED(NOPNTS) , ZCAP , VSOURC ( NOP NTS ) 



S ISP 3420 


COMPLEX YPMP 



S IS03430 


DATA PI/3 .1415926/ 



S I S0 3 4 40 


DATA NOPRT /10/ 



SIS03450 

c 




SIS03460 

c 

READ IN FROM UNIT NO. 15 THE EMBEDDING IMPEDANCES 


S IS03470 

c 

AT THE LO AND ITS HARMONICS 



SIS03480 

c 




S IS03490 


CONST = 2. * PI * FREQ * CAP 



S IS03500 


NODIV2 = NOP NTS/ 2 



SIS03510 

c 




S ISt 3520 


READ (15, 100) NCHK 



SIS03530 

100 

FORMAT! 15 > 



SIS03540 

c 

IF (NCHK.NE.NODIV2) GO TO 1000 



SISC3550 

c 




SIS03560 


READ ! 1 5 , 1 50 ) Z IMPED ( 1 ) 



S IS03570 


WR I TE ( G , 300 > NOPRT. ZIMPEDI1) 



SIS03580 

300 

FORMAT! IX, ’THE FIRST ’,13,’ LO HARM IMPEDANCES 

(EXC 

CAPAC) ARE 

: ’ S IS03590 


1/1X,’ 0 ’ ,E15.8, ’+0’ ,E15.8, ’OHMS’ ) 



S IS03600 


DO 10 I = 2 , NOD I V2 



SIS03610 


ZIMPED! I ) = CMPLX!0. 0,0.0) 



SIS03620 


IMIN1 =1-1 



S IS03630 


IFII.GT.NCHK) GO TO 2000 



SIS03G40 


READ <15,1 50 ) ZIMPED! I ) 



SIS03G50 

150 

FORMAT ! 2E 1 5 . 8 ) 



SIS03G60 

2000 

IF! I .LE. NOPRT) WRITE(G,301) IMIN1, ZIMPED(I) 



SIS03G70 

301 

FORMAT! 1X,I5,1X,E15.8,’+0’,E15.8, ’OHMS’ ) 



SIS03680 


IFU.EQ.2) GO TO 10 



SIS03G90 

c 




SIS03700 

c 

ADD IN THE JUNCTION PARASITIC CAPACITANCE CONTRIBUTION 

SIS03710 

c 




S I SP 37 20 


XCAP = -1. / CONST / FLOAT! I - 1 ) 



S IS03730 


ZCAP = CMPLX<0. ,XCAP > 



SIS03740 


IF (CABS (ZIMPED! I)).LT.1.0E-10) GO TO 10 



SIS03750 


ZIMPED! I ) = 1. / ( 1 ./ZIMPED! I > + l./ZCAP) 



S IS037G0 

10 

CONTINUE 



S IS03770 


YPMP = CMPLX! 1 .0,0.0) / Z IMP ED ( 2 ) 



S IS03780 


WR I TE ( 6 , 1 1 0 ) YPMP 



SIS03790 

1 10 

FORMAT! IX, ’Y PUMP (INCL CAPAC) = *,E15.8,’ +0 

’ ,E15 

.8, ’S’ > 

SIS03800 

c 




SIS03810 

c 

ZERO OUT THE SOURCE VOLTAGE VECTOR 



SIS03820 

c 




SIS03830 


DO 20 I = 1 , NOD I V2 



SIS03840 

20 

VSOURC! I ) = CMP LX ( 0. , 0 . ) 



SIS03850 



oooo o rs) oooo non o — oooo oooooo 


65.07 
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C 

200 


WR I TE < G , 200 ) NCHK 

F OR MAT ( IX, ’NUMBER OF POINTS READ = 
RETURN 


,15) 


1000 WRITEC6, 1001 ) 

1001 FORMAT! IX , ’ERROR IN NUMBER OF POINTS IN ZLO FILE’) 

STOP 

END 

SUBROUTINE L I NEAR { VOLT , CUR , NOP NTS , I BAS E , FACTOR , Z IMP ED , VSOURC , R I D > 
REAL VOLT! NOP NTS ) , CUR(NOPNTS) 

COMPLEX FACTOR (NOPNTS), Z I MP ED { NOPNTS ) , CTEMP(128>, VSOURC { NOPNTS 
COMMON /RBOUT/ VDC , VLO, VJDC, VULO, CURDC, CURPMP , RLOSS, RTEMP 

THIS IS THE LINEAR EMBEDDING NETWORK CALCULATION PART 
OF THE VOLTAGE UPDATE METHOD 

DO 10 I = 1 , NOPNTS 
0 CTEMP ( I ) = -CUR(I) 

CALL FFT( IBASE .CTEMP .FACTOR) 

OBTAIN THE PUMPED I -V CURVE CURRENT VALUE ON THE WAY THROUGH 
CURPMP = -REAL ( CTEMP ! 1 ) ) 

DON’T FORGET TO TAKE INTO ACCOUNT THE EFFECT OF 
THE IDENTITY ELEMENT 

N0DIV2 = NOP NTS/2 
DO 20 I = 1 , NOD I V2 

CTEMP ( I ) = CTEMP! I >*(ZIMPED( I ) *R I D ) / ! Z IMP ED ( I )+RID) + 

1 VSOURC ( I ) * RID / (ZIMPED(I) + RID) 

IF! I ,EQ. 1 ) GO TO 20 

ICORR = NOPNTS -1+2 

CTEMP ! I COR R ) = CONJG ( CTEMP ( I ) ) 

CONTINUE 

CTEMP(NODIV2+l ) = CMPLX(0.,0.) 

CALL ZEROiCTEMP, NOPNTS) 

COMPENSATE FOR THE IDENTITY ELEMENT WHEN DETERMINING 
THE PUMPED I-V CURVE CURRENT VALUE 

CURPMP = CURPMP + REAL (CTEMP! 1 ) ) / RID 

INTER = -IBASE 

C/LL FFTONTER, CTEMP, FACTOR) 

DO 30 I = 1. NOPNTS 
VOLT! I ) = REAL (CTEMP! I ) ) 

RETURN 


0 


30 

C 


SIS038G0 
SIS03870 
S IS03880 
SIS03890 
S IS03900 
SIS03910 
S IS03920 
SIS03930 
SIS03940 
SIS03950 
S IS039G0 
SIS03970 
SIS03980 
)S IS03990 
S IS04000 
SIS04010 
SIS04020 
S IS04030 
SIS04040 
SIS04050 
SIS04060 
SIS04070 
S1SC4080 
SIS04090 
SIS04100 
SIS041 10 
SIS04I20 
SIS04130 
SIS04140 
SIS04150 
SIS04160 
SIS04170 
SIS04180 
SIS04190 
SISD4200 
SIS04210 
S I SP 4 220 
SIS04230 
SIS04240 
SIS04250 
SIS04260 
S IS04270 
SIS04280 
S IS04290 
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SIS04320 
S I S0 4 330 
S IS04340 
SIS04350 
SIS04360 
SIS04370 
SISC 4380 
SIS04390 
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65.08 
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C 


C 

C 

C 

C 


C 


C 

C 

C 


C 

C 

c 

1 00 


c 

c 

c 


200 


c 

c 

c 

c 

c 

c 


END 

SUBROUTINE S ISNL ( VOLT ,CUR , NOPNTS , I BASE 1 , CH I , NOCH I , RN , FACTOR , FREQ, 
1RID) 

REAL VOLT! NOPNTS > , CHKNOCHI), CUR(NOPNTS) 

COMPLEX UTERM! 1024 ) , FACTOR ( NOP NTS > , CVOL(128), CSUM, CTEMP 
COMPLEX CCH I ( 1024 ) 

DATA PI.E.HBAR/3. 1415926. 1.602 E- 19, 1.05 5E-34/ 

THIS IS THE SIS NONLINEAR NETWORK CALCULATION SUBROUTINE 
SUBROUTINE PART OF THE MULTIPLE REFLECTION TECHNIQUE 

ITST1 = 10 

IFINOCHI .EQ. 1024) ITST1 = 11 
ITST2 = -ITST1 

IF (NOCHI .GT. 1024) GO TO 1000 

I F ( NOCH 1 . NE . 1024 .AND- NOCH I . NE . 5 1 2 > GO TO 1002 

SET UP A FEW PARAMETERS 

OMEGA = 2. * PI * FREO 
N0DIV2 = NOPNTS/2 
INTER = -I BASE 1 

CONST = -1. * E / HBAP. / OMEGA 
HSTEP = 1. / FREQ / F L OAT < NOP NTS ) 

CALCULATE THE DC SEPARATELY 

DO 100 I = 1, NOPNTS 
CVOL ( I ) = VOLT < I ) 

CALL FFT( IBASE1 , CVOL, FACTOR) 

RAMP = REA LI CVOL (1)) * (-1.) * E / HBAR 

INTEGRATE THE NON-DC COMPONENTS 

DO 200 I = 2 , NOD I V2 
ICORR = NOPNTS -1+2 

CVOL I I ) = CVOL I I ) / CMP LX ( 0 . , 1 . ) / FLOAT! I - 1 > * CONST 
CVOL I I CORR ) = CONOG I CVOL I I > ) 

CONTINUE 

CVOL ( 1 ) = CMPLXI0 . ,0. ) 

CVOL I NOD I V2+ 1 > = CMPLXI0. ,0. ) 

CALL FFTI INTER, CVOL , FACTOR ) 

NTIME = NOCHI / NOPNTS + 1 

CALCULATE U 

NTIME = 2 * NOCHI / NOPNTS 
DO 300 0 = 1 .NTIME 
DO 300 I = 1, NOPNTS 
I PTR = NOPNTS * IJ-1) + I 

TEMP = REAL! CVOL < I ) ) + FLOAT! I PTR- 1 ) * HSTEP * RAMP 


SIS04410 

S IS04420 

SIS04430 

SIS04440 

SIS04450 

S IS04460 

SIS04470 

S IS04480 

SIS04490 

SIS04500 

SIS04510 

SIS04520 

SIS04530 

S IS04540 

SIS04550 

SIS04560 
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IF IABSITEKP ) . L E . 2 . *PI ) GO TO 250 
IREVS = TEMP/2. /PI 

TEMP = TEMP - FLOAT! IREVS) * 2. * PI 
250 CTEMP = CMPLX < 0. , TEMP ) 

UTERM(IPTP) = CEXP ! CTEMP > 

300 CONTINUE 
C 

DO 5000 IMF = 1, NOCHI 
CCHI ( IMF ) = CHI! IMF) 

CCHI ! IMF+NOCHI ) = 0.0 
5000 CONTINUE 
C 

C DO THE CONVOLUTION CALCULATIONS IN THE FREQUENCY DOMAIN 

C 

CALL FFT! ITST1 , CCHI .FACTOR ) 

CALL FFT! ITST1 ,UTERM, FACTOR) 

C 

NOCHM2 = 2 * NOCHI 
DC 6000 IMF = 1 , NOCHM2 

CCHI ! IMF ) = CCHIIIMF) * UTERM! IMF ) * FLOAT! NOCHM2 ) 

6000 CONTINUE 
C 

CALL FFT! ITST2, CCHI .FACTOR) 

CALL FFT! ITST2, UTERM, FACTOR) 

C 

C PERFORM THE FINAL CALCULATIONS AND CORRECT FOR THE 

C IDENTITY ELEMENT VALUE 

C 

DO 600 I = l.NOPNTS 
N = 2*NOCh I - NOP NTS + I 
C CSUM = Ct PLX(0. .0. ) 

C N = NOCHI + I 

C NOMINI = NOCHI - 1 

C NOMIN3 = NOCHI - 3 

C DO 500 K = 1 .NOMIN3.2 

C CSUM = (CHI! K >*UTERM( N-K+l )+4 .0*CHI(K+I >*UTERM! N-K>+ 

C 1CHI (K+2)*UTERM(N-K-1 ) )/3.0 + CSUM 

C500 CONTINUE 

C CSUM = CSUM+!CHI ( NOM INI )*UTERM!N-NOCKI + 2> + 4 . 0*CH I {NOCHI ) 

C 1 *UTERM ( N-NOCH I + 1 ) ) /3 0 

CSUM = CCHI ( N > 

CSUM = CSIM * CONUG ( U' ERrl! N > ) * HSTEP 
CUR ( I ) = VOLT! I > / RN + AIMAG!CSUM) - VOLT!I) / RID 
600 CONTINUE 
C 

RETURN 

C 

1000 WRITE! 6 , 1001 ) NOCHI 

1001 FORMAT! IX. ’NOCHI = ’,15,’ AND IS TOO BIG FOR THE AVAIL MEMORY’) 
STOP 

C 

1002 WRITE! 6, 1003) NOCHI 

1003 FORMAT! IX. ’NOCHI = ’,15,’ AND IS AN INCORRECT VALUE’) 

STOP 

END 
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SUBROUTINE COMPAR I LATEST , OL D , NOPNTS , FRE I TR > 

REAL LATEST! NOPNTS ) , OLD(NOPNTS) 

LOGICAL FREITR 
COMMON /PERIPH/IN, IOUT 
COMMON /ERRALL/ TOL , EPS 
C 

C AUTHOR: R.G. HICKS. UNIV. OF OLD . , APRIL 1977 
C 

C AIM: THIS SUBROUTINE COMPARES TWO VECTORS ’LATEST’ 

C AND ’OLD’ AND PROVIDES A LOGIC ANSWER ’FREITR’ 

C AS TO WHETHER THE TWO VECTORS ARE SUFFICIENTLY EQUAL. 

C 

FREITR = .FALSE. 

ERR = 0.0 
C 

TOP = 0. 

DO 5 I = 1 .NOPNTS 

IF !ABS!LATEST! I ) ) . GT . TOP ) TOP = ABS ( LATEST ( I > ) 

5 CONTINUE 

C 

DO 10 1 = 1 , NOF NTS 

IF !ABS(LA7EST! I ) ) .LE . EPS*TOP > GO TO 10 
TEST = ( LATEST ( I ) - OL D ( I ) ) /LATEST { I > 

IF (ABS(TEST) .GT.TOL) FREITR = .TRUE. 

I F (ABS ! TEST ) .LE. ERR) GO TO 10 
ERR = ABS ( TEST ) 

MEM = I 

10 CONTINUE 

C 

WRITE! IOUT. 20) ERR, MEM 

20 FORMAT! 8X ,’ ERP = ’ , E 1 5 . 8 , ' POSN NO= ’,15) 

C 

RETURN 

END 

SUBROUTINE TOP S ( FREO , F I F , I BASE 2 . NOl V , I BASE 1 , NOPNTS , MMAX , CAP , 

1RN , TEMP ,NMAX, VOLT, Y0,Y1 .Y1M.ZLO, I VDAT , 0 , FACTOR , DELTAV, ISIG.ZSSB) 
REAL I VDAT ! NO 1 V ) , VOLTINOPNTS) 

COMPLEX YSIG131.31), NO I SE ( 3 1 , 3 1 ) , ZSIGI31.31), ZOUT 
COMPLEX B(31>, STORE (31,31) 

COMPLEX FACTOR ( NOIV) , Y0, Y1.Y1M, Z0, ZSSBIMMAX) 

COMPLEX J(NOIV) , W! 128) 

COMPLEX YS12(3, 3), N012(3,3), ZSI2(3,3) 

COMPLEX YCAP , ZLO 

COMMON /PERIPH/ IN, IOUT, IPLOT 

COMMON /RBOUT / VDC , VLO , VODC, VJLO, CURDC, CURPMP, RLOSS, RTEMP 
C 

C THIS IS THE MAIN SMALL SIGNAL CONVERSION AND NOISE 
C CONTROLLING CALCULATION ROUTINE 

C 

WFITE<6,709) V J DC , VJLO 

709 FORMAT! IX , ’DC VOLTAGE AT THE J UNCT I ON ’ , 8X , ’ = ’.E15.8,’ VOLTS’/ 
llX.’LC VOLTAGE (PEAK) AT THE JUNCTION = ’.E15.8,’ VOLTS’) 

C 

WP ITECS, 1 ) 

1 FORMAT! /IX , ’SMALL SIGNAL AND NOISE ANALYSIS SECTION’/) 
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C 

IF(ISIG.EQ.l) GO TO 5500 
C 

C FOR THE GENERAL ANALYSIS METHOD 

C 

DO 10 I = 1 , MMAX 
B < I ) = CMP LX (0 . ,0. ) 

DO 10 II = 1 , MMAX 

ZSIG< I , I I > = CMPLX (0. .0. ) 

STORE < 1,11) = CMPLX ( 0 . ,0. ) 

NOISE! 1 ,11) = CMPLXC0. ,0. ) 

YSIG! I ,11) = CMPLXC0. ,0. ) 

10 CONTINUE 

C 

WR ITE ( S , 504 > 

504 FORMAT! /IX, ’ DETERMINING SMALL SIGNAL MIXER PROPERTIES USING’, 

1* THE GENERAL ANALYSIS’) 

C 

CALL S 1 GNL ! I VDAT ,NOIV,FREQ,I BASE 1 , I BASE 2 , VOLT , NOP NTS , YS IG , MMAX , 
IF IF, DELTAV.O.W, FACTOR, NMAX.RN) 

C 

CALL SIGL I N ( FREQ , YS IG , MMAX ,FIF,CAP,ZSSB ) 

C 

CALL I NVRT ! MMAX ,YSIG,ZSIG,B, STOR E ) 

C 

CALL LOSCAL ( MMAX ,YSIG,ZSIG,Y0,Y1 ,Y1M> 

C 

CALL SISNOI ( W, NMAX , MMAX ,NOIV,FREQ,FIF ,VUDC, 

10 .DELTAV, TEMP .NOISE , NOP NTS ,RN> 

CALL P ESNOI! TEMP, NOISE, MMAX, FREQ, F IF, CAP, ZSSB) 

CALL PR OCES!ZSIG, NOISE , MMAX , Y1 ,Y1M> 

C 

GO TO 5000 
C 

C USING THE APPROXIMATE THREE-PORT MODEL 

C 

5500 DO 15 I = 1,3 
DO 15 II = 1,3 
NO 1 2 ( I , I I > = CMPLX ! 0 . ,0 . ) 

ZS 1 2 ! I , I I ) = CMPLX !0. ,0. > 

YSI2! I , II ) = CMPLX ( 0 . ,0. ) 

15 CONTINUE 

DO 17 I = l.MMAX 

B! I > = CMPLX !0. ,0. ) 

DO 17 II = l.MMAX 

STORE ( I , 1 1 ) = CMPLX (0. ,0. ) 

17 CONTINUE 

C 

WR I TE ! 6 , 505 ) 

505 F0RMAT!/1X, ’THREE PORT SMALL SIGNAL ANALYSIS USING FORMULAE’) 

C 

CALL THEORY! VOL 0, VO DC , NMAX , I VDAT ,NOIV,0,FREQ,YSI2, NOI 2 .TEMP , 

1 DELTAV, RN.Y0.Y1 .Y1M.ZL0) 

C 

4999 WR I TE ( 6 , 506 ) 
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506 F0RMAT(/1X, ’THREE PORT ANALYSIS USING MATRIX TECHNIQUES’) 

C 

CALL SIGLIN(FREQ,YSI2,3,FIF .CAP ,ZSSB> 

CALL INVRT(3, YSI2.ZSI2.B, STORE) 

CALL LOSCAL(3,YSI2,ZSI2, :0,Y1 ,Y1M) 

CALL RESNOI (TEMP, NOI2, 3, FREQ, FIF , CAP ,ZSSB) 

CALL PROCES(ZSI2,NOI2,3,Yl ,Y1M) 

C 

5000 RETURN 
C 

END 

SUBROUTINE RESNOI ( TEMP , NOl SE , MMAX , FREQ, F I F , CAP , ZSSB > 

REAL K 

COMPLEX NOISE ( MMAX , MMAX ) , YPAP. , ZSSB(MMAX) 

DATA HBAR/1 .055E-34/ 

DATA PI/3.1415926/ 

DATA K/l . 380622E-23/ 

C 

C DETERMINE THE THERMAL NOISE COMPONENT EMANATING FROM 

C THE RESISTANCES USING THE QUANTUM CORRECTED FORMULA 

C NO NOISE FROM THE EMBEDDING IMPEDANCES AT THE SIGNAL 

C AND IMAGE IS CONSIDERED IN THE DEFINITION 

C 

DO 100 I = 1 , MMAX 

M = (MMAX/ 2) - I + 1 

I PTR = I APS ( M ) * 2 

IF ( M . L E . 0 > I PTR = I PTR + 1 

I F ( CABS (ZSSB( IPTR) > .LT. 1 .E-10) GO TO 500 

YPAR = 1 . / ZSSB ( IPTR) 

GO TO 600 

500 YPAR = CMPLX(0.0, 1 .E10) 

600 I F ( IABS ( M ) . EQ . 1 ) YPAR = CMPLX (0.0, 0.0) 

I F ( M . EQ . 0 ) YPAR = CMP LX ( 0 . 0 ,0 . 0 > 

OMEGA!', * ( FLOAT ( M * FREQ + FIF) * 2.0 * PI 
VAL = HBAR * CMEGaM / K /TEMP 

NOISE(I.I) = f CISC ( I , I ) + 4 .0*HBAR*OMEGAM*REAL (YPAR > 

1 / (EXF(VAL) - 1.0) 

100 CONTINUE 
RETURN 
END 

SUBROUTINE PROCES ( ZS IG , NOI SE , MMAX , Y 1 , Y1 M ) 

REAL K 

COMPLEX STORE ( 3 1 ) . SUM 

COMPLEX Yl, Y1M, NOISE ( MMAX .MMAX ) , ZS I G ( MMAX , MMAX ) 

COMPLEX ZIN 
INTEGER SGPORT 

COMMON /R BGU1 / VDC, VLO, VODC, VOLO, CURDC, CURPMP , RLOSS, TEMPI 
DATA K/l 380622E-23/, TOL/l.E-4/ 

C 

C DETERMINE THE NOISE TEMPEPATURE USING THE IMPEDANCE 

C MATRIX AND THE NOISE CONTRIBUTIONS MATRIX 

C 

I F PORT = MMAX / 2 + 1 
SGPORT = MMAX / 2 
C 
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DO 20 I = l.MMAX 
SUM = CMPLXI0. ,0. ) 

DO 10 0 = l.MMAX 

SUM = SUM + NOISE(I.J) * CONJG ( ZS I G < I F PORT , 0 ) > 

10 CONTINUE 

STORE! I ) = SUM 
20 CONTINUE 
C 

SUM = CMPLX (0. ,0. ) 

C 

DO 30 I = l.MMAX 

SUM = SUM + ZSIG( IFPORT, I ) * STORE(I) 

30 CONTINUE 

C 

VREAL = REAL (SUM > 

VIMAG = A I MAG (SUM) 

I F ( ABS ( VI MAG ) . GT . TOL *ABS { VREAL ) ) GO TO 100 
C 

VAL1 = CAES(ZSIG( IFPORT, SGPORT) ) 

VAL2 = r.A?S(ZSIG< IFPORT, IFPORT+1 ) ) 

C 

C POWER = REAL (SUM) / 50.0 

C WR I TE ( 6 , 50 ) PO'.EF 

C50 F OR MAT { EX , ’ I F OUTPUT POWER INTO 50 OHM LOAD = ’ .E15.B) 

C 

TEMPI = REAL { SUM ) / (4.*K*REAL(Y1 ) *VAL 1**2) 

TEMP 2 = REAL ( SUM ) / < 4 . *K*REAL ( Y 1 M ) * VAL 2* * 2 ) 

WP I TE ( 6 , 40 ) TEMPI, TEMP2 

40 F CiR MAT ( /8X , ’ SS B TEMP OF MIXER (USB) = ’,£15.8,’ DEG K’/ 

1 8X , ’ SSB TEMP OF MIXER ( LSB ) » ’,E15.8,’ DEG K’) 

C 

RETURN 

C 

100 WR ITE ( S , 1 10) VREAL, VIMAG 

110 FORMAT( IX, ’PROBLEM IN TEMPERATURE CALCULATIONS:’/ 

1 IX, ’TEMP = ’,E15.G,’ +0 ’ , E 1 5 . 8 > 

RETURN 

END 

SUBROUTINE S I GNL ( I VDAT , NO I V , FR EQ , I BASE 1 , I BASE2 , VOLT , NOPNTS , YS IG , 
1 MMAX.F IF, DELTA V.U.W, FACTOR, NMAX.RN) 

REAL I VDAT ( NO I V ) , VOLT(NOPNTS) 

COMPLEX C VOL (128). J(NOIV), W( NOPNTS), YS I G { MMAX , MMAX ) 

COMPLEX C VAL 1 , CVAL2 , CVAL3, CVAL4, CSUM 
COMPLEX FACTOR(NOIV) 

COMMON /RBOUT/ VDC, VLO, VJDC, V0LO, CURDC, CURPMP, RLOSS, RTEMP 
DATA P I , E , HBAR/3 .1415928,1 . S02E- 1 9 , 1 .055E-34/ 

C 

C DETERMINE THE ADMITTANCE CONVERSION COMPONENTS DUE TO 

C THE SIS TUNNEL JUNCTION AS PER TUCKERS THEORY 

C 

OMEGA = 2. * PI * FREQ 

N0DIV2 = NOPNTS/2 

INTER 1 = - I BASE 1 

INTER2 = -IBASE2 

CONST = -1. * E / HBAP / 'EGA 
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HSTEP = 1. / FREQ / FLOAT< NOPNTS > 

WIF = 2. * PI * FIF 
C 

DO 100 I = 1, NOPNTS 
100 CVOLU ) = VOLT ( I ) 

CALL FFTII BASE 1 , C VOL , FACTOR ) 

C 

VB IAS = REAL ( CVOL ( 1 ) ) 

C 

DO 200 I = 2 , NOD I V2 
ICORR = '' PNTS -1+2 

CVOLU) * - VOL ( I ) / CMPLX (0. , 1 . ) / FLOAT! I - 1 ) * CONST 
CVOL ( I COR) ' = CONOG { CVOL ! I > 

200 CONTINUE 
C 

CVOLU ) = CMPLX (. ,0. > 

CVOL ( NOD I V2+ 1 ) = C.MP LX ! 0. , 0 . ) 

C 

CALL FFTUNTER1, CVOL, FACTOR) 

C 

DO 300 I = 1, NOPNTS 
TEMP = REAL ( CVOL < I > ) 

CVOLU ) = CMP LX ( 0. , TEMP > 

WU ) = CEXP < CVOL ( I ) ) 

300 CONTINUE 
C 

CALL FFTUBASE1 ,W, FACTOR) 

C 

VNORM = E * VBIAS / HBAR 
C 

DO 1000 Ml = 1 , MMAX 
DO 1000 M2 = 1 , MMAX 
M = MKAX/2 - Ml + 1 
MDASH = MMAX/2 - M2 + 1 
I D I F = M - MDASH 

WMDASH = FLOAT { MDASH > * OMEGA + WIF 
CSUM = CMPLX ( 0 . ,0. ) 

DO 500 N 1 = 1 , NMAX 

NDASH = NMAX/ 2 - N1 + 1 

N = NDASH - IDIF 

IFUABS(N) .GT.NMAX/2) GO TO 500 

ARG 1 = FLOAT < N > * OMEGA + VNORM 

ARG2 = FLOAT! N) * OMEGA - WMDASH + VNORM 

ARG3 = FLOATINDASH) * OMEGA + WMDASH + VNORM 

ARG4 = FLOAT ! NDASH ) * OMEGA + VNORM 

CALL EVALO " ARG 1 , CVAL 1,0,1. 'IV, DELTAV ,RN > 

CALL EVALO ( AkG 2 , CVAL 2 , 0 ,'.v I V , DELTAV , R N ) 

CALL EVALO (ARC 3, CVAL3,0.N‘- IV, DELTAV, RN> 

CALL EVALO ( ARC 4 .CVAL 4.0 ,NOIV, DELTAV, RN) 

IF ( ABS ( WMT ASH ) . GT 1 .0E9'2 .0*P I ) GO TO 4355 
CALL DERIO (ARG1 ,C\AL 1 , 0 , NO I V , DE LTAV , R N > 

CALL DERI0(ARG4,CVAL4,0,NOIV, DELTAV ,RN ) 

CVAL4 = Cl PLX(-1 .X ,0.0) * CVAL 4 
CVAL2 = Cl. PLX<0. 0,0.0) 

CVAL 3 = Cl, PLX(0. 0,0.0) 
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4356 NPTR1 = N + 1 

IF< IABS(N) .GE.NOPNTS/2) GO TO 5 00 
IFINPTRl . LE . J0T) NPTR1 = NPTR1 + NOPNTS 
NPTR2 = NDASH + 1 

I F ( IAES(NDASH) .GE.NOPNTS/2) GO TO 500 
IF ! NPTR2. LE.0) NPTR2 = NPTR2 + NOPNTS 

CSUM = W(NPTRl) * CONOG ( W( NPTR2 ) ) * ( CVAL I -CVAL2-CONOG ( CVAL 3 > 
t 1 CONJG ( CVAL4 ) ) + CSUM 
500 CONTINUE 

YS I G ( Ml , M2 ) = CMPLX(0. ,-l . )*E/2./HBAR/WMDASH*CSUM 
IF lABS(WMDASH) .GT. 1 .0E9"2.0*P1 ) GO TO 1000 
Yo IG ( Ml , M2 ) = CMPLX(0 0,-1. 0) / 2.0 * CSUM 
C 1 963 WR ITE ( 6 , 2P3 1 ) Ml, M2, YSIG(M1,M2> 

C2831 FORMAT! 1X,’M1=’,I5,’M2=’,I5,’YSIG=’,2E15.8> 

1000 CONTINUE 
RETURN 
END 

SUBROUTINE GETU! I VDAT , NO I V , FREQ , I BASE 2 , 0 , F ACTOR , RN , DELTAV ) 
COMPLEX CTEMP ( 2048 ) , FACTOR < NO I V > , O(NOIV) 

REAL I VDAT ( NO I V ) 

C 

C OBTAIN THF VECTOR 0 BY PERFORMING A KRONERS KRONIG 

C TRANSFORM ON THE REDUCED I-V CURVE AS PER THE REPORT 

C 

IBASE3 = IBASE2 + 1 
NOUSE = 2 * NO I V 
INTER3 = -IBASE3 
C 

DO 100 I = 1 , N 0 1 V 
CTEMP! I) = I VDAT ( I ) 

IF! I .EQ. 1 ) GO TO 100 
ICORR = NOUSE -1+2 
CTEMP ( I COP R ) = - I VDAT ! I ) 

100 CONTINUE 

CTEMP ( NO I V+ 1 ) = CMP LX ! 0. , 0 . ) 

C 

CALL FFT! IBASE3, CTEMP .FACTOR) 

C 

DO 200 I = 2, NO IV 

CTEMP ( I ) = CMPLX(0.,1.) * CTEMP ! I > 

ICORR = NOUSE -1+2 

CTEMP ( ICORR ) = CMPLXI0. ,-l . ) * CTEMPIICORR) 

200 CONTINUE 
C 

CTEMP! 1 ) = CMP LX ! 0 . ,0. ) 

CTEMP ( NOIV+1 ) = CMPLX(0. ,0. ) 

C 

CALL FFT! INTER3, CTEMP, FACTOR) 

C 

DIFF = REA L (CTEMP ( 1 ) ) 

C 

DO 300 I = 1 , NO IV 

TEMP = REA L ( CTEMP ( I ) ) - DIFF 

HOLD = I VDAT! I ) 

0(1) = CMPLXtTEMP .HOLD) 
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300 CONTINUE 
C 

RETURN 

C 

END 

SUBROUTINE S IGL I N ( FREQ . YS IG ,MMAX , F IF , CAP ,ZSSB ) 

COMPLEX YS I G ( MMAX , MMAX ) , ZSSB(MMAX), CVAL 
C 

C OETAIN THE CONTRIBUTIONS TO THE SIGNAL CONVERSION MATRIX 

C FROM THE EMBEDDING IMPEDANCES 

C 

DO 2000 Ml = 1 , MMAX 
M = MMAX / 2 - Ml + 1 
C 

I PTR = IAES(M) * 2 
I F ( M . L E . 0 ) I PTR = I PTR + 1 
IF (CABS(ZSSB( IPTR ) ) .LT. 1 .E-10) GO TO 500 
CVAL = 1 . / ZSSB ( I PTR > 

GO TO 600 

500 CVAL = CMFLX( 1 .E10, 1 .E10) 

600 YSIG(Ml.Ml) = YS IG ( Ml , M 1 ) + CVAL 
2000 CONTINUE 
RETURN 
END 

SUBROUTINE LOSCAL ( MMAX , YS IG , ZS IG , Y0, Y1 , Y1M ) 

C 

INTEGER SGPORT 
REAL LOSS, LOSSL 
REAL MSMCHI , MSMCHO 

COMPLEX Zll, ZIN, CTMP1, ZSOURC, ZOUT 

COMPLEX Y£ , Y1 , YS IG ( MMAX , MMAX ) , ZS IG ( MMAX , MMAX ) , Z01 , Z00, RHO 
COMPLEX Y1M, Z0M 1 

COMMON /REOUT/ VDC, VLO, VODC, VJLO, CURDC, CURPMP , LOSS, RTEMP 
C 

C CALCULATE THE CONVERSION LOSS FROM THE SMALL SIGNAL 

C IMPEDANCE CONVERSION MATRIX 

C 

RSOURC = 1 .0/REAL ( Y 1 ) 

R0 = 1 .0 / REAL ( Y0) 

C 

I F PORT = MMAX/2 + 1 
SGPORT = MMAX/2 
IMPORT = MMAX/2 + 2 
Z01 = ZSIGC IFPORT, SGPORT) 

Z00 = ZSIGC IFPORT, IFPORT) 

Z0M1 = ZSIGC IFPORT, IMPORT) 

Zll = ZSIGCSGPORT, SGPORT) 

C 

ZOUT = 1.0 / (1.0/Z00 - 1.0/R0) 

RHO = (ZOUT - R0) / (ZOUT + R0) 

VSWR = C1.0 + CABS (RHO) ) / (1.0 - CABS(RHO)) 

WR ITE ( 6 , 555 > ZOUT, R0, VSWR 

555 F ORMAT ( / 8X , ’ I F OUTPUT IMPEDANCE = ’,E15.8,’ +U ’,E15.8/ 

1 8X , ’ I F VSWR REFERRED TO ’,F5.2,’ OHMS = ’,E15.8) 

C 


SIS08810 
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ZIN = 1.0 / (1.0/Z11 - 1.0/RSOURC) 

CTMP1 = RSOURC + ZIN 

MSMCHI = 4 . 0*RSOUP,C*REAL ( Z I N ) / ( CABS I CTMP 1 ) **2 ) 

DB1 = ABS ( MSMCH I ) 

DB1 = 10.0 * ALOG10IDB1 ) 

IF (MSMCHI .LT.0.0) DB1 = 0.0 
C 

CTMP 1 * R0 + ZOUT 

MSMCHO = 4 . 0*R0*R EAL ( ZOUT > / ( CABS (CTMP1)**2) 

DE2 = ABS ( MSMCHO > 

DE2 = 10.0*ALOG10(DB2) 

IF (MSMCHO. LT.0.0) DB2 = 0.0 
C 

WR I TE ( 6 , 4000 ) MSMCHO, DB2 , ZIN, MSMCHI, DB1, RSOURC 
4000 FORMAT(8X, ’OUTPUT MISMATCH GAIN = ’.E15.8,’ I ’ , E 1 5 . 8 , ’ DB ) ’ // 
18X, ’INPUT SIGNAL IMPEDANCE = ’,E15.8,’ +0 ’ , E 1 5 . 8 , ’ OHMS ’ / 
28X, ’INPUT MISMATCH GAIN = ’,E15.8,’ ( ’ , E 1 5 . 8 , ’ DB ) ’ / 

38X, ’INPUT GAIN CALCULATION REFERRED TO = ’.E15.8,’ OHMS’) 

C 

LOSS = 0. 25/CABS ( Z01 ) **2/REAL ( Y 1 ) /REAL ( Y0) 

LOSS = -10. * ALOG 10( LOSS ) 

LOSSL = 0.25/CABS(Z0Ml)**2/REAL(YlM)/REAL(Y0) 

LOSSL = -10.0 * ALOG 10( LOSSL ) 

WRITERS, 3000) LOSS, LOSSL 

3000 FCRMAT(/8X, ’CONVERSION GAIN (USB) EQUALS ’.E15.8,’ DB’/ 

18X , ’CONVERSION GAIN < LSB > EQUALS ’,E15.8,’ DB’) 

RETURN 
END 

SUBROUTINE EVALO ( DUM , C VAL , 0 , NO I V , DELTAV , RN > 

COMPLEX OVAL , J(NOIV), TERP , TERP2 , TERP3, TERP5 
LOGICAL FLAG 

DATA E, HEAR/1 .602E-19, 1 .055E-34/ 

EVALUATE THE REQUIRED VALUE OF J BY INTERPOLATING 
THROUGH THE ARRAY OF CALCULATED 0 VALUES. 

FOR NEGATIVE ARGUMENTS, THE CONJUGATE VALUE OF J 
MUST BE USED. 

FLAG = .FALSE. 

ARG = DUM 

IF(ARG) 100,400,300 
100 FLAG = .TRUE. 

ARG = -DUM 

300 PTR = ARG / DELTAV / E * HBAR 
I PTR = PTP + 1 . 

IF ( IPTR .GE . NOIV-1 ) GO TO 200 
RATIO = AMOD ( PTR , 1 . ) 

C VAL = TEPP5(J ,NOIV, IPTR, RATIO, .FALSE. ) 

TEMP = PTF * TELTAV / RN 
C VAL = CV/ L + r.4P LX (0 . , TEMP > 

IF ( FLAG 5CVAL = CONOG(CVAL) 

RETURN 

200 TEMP = PTR * DELTAV / RN 

C VAL = J(NOIV) + CMPLX(0. ,TEMF > 

I F ( FLAG ) C VAL = CONOG(CVAL) 
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RETURN 

400 CVAL = CMPLX (0. ,0. ) 

RETURN 

END 

SUBROUTINE SISNOI ( W, NMAX , MMAX , NO I V , FREQ , F I F , VB I AS , 0 , DELTAV , 
1TEMP, NOISE, NOPNTS.RN) 

COMPLEX W(NOPNTS), NOISE (MMAX ,MMAX ) , J(NOIV>, CVAL3 , CVAL4 
DATA PI , E , HBAR/3. 1415926, 1 . 602E-19 , 1 .055E-34/ 

DATA BOLTZ/1 .380622E-23/ 

EVALUATE THE NOISE COMPONENT EMANATING FROM THE 
SHOT NOISE IN THE TUNNEL JUNCTION AS PER TUCKERS 
THEORY. 


BETA = 1 . / BOLTZ / TEMP 
OMEGA = 2. * PI * FREO 
WIF = 2. * PI * FIF 
CONST = E / HBAR 
C 

DO 1000 Ml = 1 , MMAX 

DO 1000 M2 = 1 , MMAX 

M = MhAX/2 - Ml + 1 

MDASH = Mt'AX/2 - M2 + 1 

WMDASH = FLOAT {MDASH ) * OMEGA + WIF 

I ADJ = M - MD/SH 

DO 500 N 1 = :.NMA> 

N = NMAX/2 - N1 + 1 
NDASH = N + I A D J 

ARG1 = BETA*! E * VB I AS + F LOAT ( NDASH > *HBAR*OMEGA + HBAR*WMDASH ) /2 . 
ARG2 = BETA*(E*VBIAS+FLOAT(N)*HBAR*OMEGA-HBAR*WMDASH>/2. 

ARG3 = VBIAS+FLOAT(NDASH)*HBAR*OMEGA/E + HBAR*WMDASH/E 
ARG4 = VB I AS + FLOAT ( N ) *HEAR*OMEGA/£ - HBAR*WMDASH/E 
ARG3 = CONST * ARG3 
APG4 = COT ST * ARG4 

CALL EVALJ (ARG3 .CVAL3.0 , NO I V , DE LTAV , RN ) 

CALL EVALJ (ARG4 , CVAL 4 ,0 , NO I V , DELTAV , RN ) 

VAL3 = AIMAG (CVAL3 ) 

VAL4 = AIMAG (CVAL 4) 

VAL1 = COTH ( ARG 1 ) 

VAL2 = COTH ( ARG2 ) 

NPTR1 = N + 1 

I F ( NPTR 1 . L E . 0 > NPTR1 = NPTR1 + NOPNTS 
NPTR2 = NDASH + 1 

IF ( NPTR2 . LE .0) NPTR2 = NPTR 2 + NOPNTS 

NO I SE ( Ml , M2 > = NOI SE ( Ml , M2 ) + E *W< NPTR 1 > *CONOG ( W( NPTR2 ) > * ( VAL 1 
1 * VAL 3 + VAL 2 * VAL 4 ) 

500 CONTINUE 
1000 CONTINUE 
RETURN 
END 

SUBROUTINE I NVRT ( NOS IG , YS I G , ZS I G , B , STORE ) 

C 

COMPLEX YSIG(NOSIG,NOSIG) , ZSIG(NOSIG), B(NOSIG) 

COMPLEX STORE ( NOS IG , NOS IG > 

C 
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SIS09920 
SIS09930 
S IS09940 
SIS09950 
SIS09960 
SIS09970 
SIS09980 
S IS09990 
SIS 10000 
SIS 10010 
SIS 10020 
SIS 10030 
SIS 10040 
SIS 10050 
SIS 10060 
SIS 10070 
S I S 10080 
SIS 10090 
SIS10100 
SIS101 10 
S IS 101 20 
S IS 101 30 
S I S 1 0 1 40 
S I S 1 0 1 50 
S I S 1 01 60 
S IS 101 70 
S I S 1 0 1 80 
S IS 101 90 
S I S 1 020£ 
SIS 102 1 0 
SIS 10220 
SIS 10230 
SIS 10240 
SIS 10250 
S I S 1 0260 
SIS 10270 
SIS 10280 
SIS 10290 
SIS 10302 
SIS 103 10 
S IS 1 0320 
SIS 10330 
SIS 10340 
S I SJ 0350 
S I S 1 0360 
SIS 10370 
S I S 1 0380 
S I S 1 03 90 
S I S 1 0400 
S IS 1 04 10 
S I S 1 0420 
SIS 10430 
SIS 10440 
S I S 10450 




65.19 


FILE: 

SISCAP FORTRAN A NASA/GISS CONVERSATIONAL MONITOR 

SYSTEM SP 302 


COMMON /PER I PH/ IN, IOUT 

SIS 10460 

C 


SIS 10470 

c 

THIS ROUTINE INVERTS THE SQUARE MATRIX 

SIS10480 

c 


S I S 104 90 


NSQ = NOSIG * NOSIG 

S I S 1 0500 

c 


S I S 1 05 10 


DO 5 I = 1, NOSIG 

SIS 10520 


DO 5 J = 1 .NOSIG 

SIS 10530 

5 

STORE ( I , J ) = YSIG(I.O) 

SIS 10540 

c 


S I S 1 0550 


CALL FACTOR (STORE, STORE, NOSIG, I FLAG) 

SIS 10560 

c 


S I S 1 0570 


GO TO ( 20 , 1 1 ) , I FLAG 

SIS 10580 

c 


SIS 10590 

1 1 

WRITE! IOUT, 611 ) 

SIS 10600 

611 

FORMAT ( IX , ’ THE MATRIX IS SINGULAR’ > 

S I S 1 06 1 0 


RETURN 

S I S 1 0621 

C 


S I S 1 06 30 

2J0f 

DO 21 I = 1 .NOSIG 

SISI0640 

21 

B ( I > = CMPLX (0. ,0. ) 

SIS. 06 50 

C 


SIS. 0660 


IBEG = 1 

SIS 10670 


DO 30 0 = 1 .NOSIG 

S I S 1 06 80 


B ( J ) = CMPLX ( 1 . ,0. ) 

S I S 1 06 90 


CALL SUBST(STORE,E,ZSIG( IBEG) .NOSIG) 

SIS 10700 


B ( 0 ) = 0. 

SIS 107 10 

30 

I EEG = IBEG + NOSIG 

S I S 1 0726 

C 


SIS:0730 


RETURN 

SIS 0740 


END 

SIS'0750 


SUBROUTINE FACTOR ( A , W , N , I F LAG > 

S I S 1 07 60 

C 


SIS 10770 


COMPLEX RATIO 

SIS 10780 


COMPLEX A ( N , N > , W(N,N) 

SIS 10790 


DATA TOL/1 .E-7/ 

SIS 10806 

c 


SIS 108 1 0 

c 

THIS SUBROUTINE FACTORS A GIVEN MATRIX INTO 

S I S 1 08 20 

c 

LU FORM USING GAUSSIAN ELIMINATION. 

SIS 10830 

c 

SEE CONTE AND DE EOOR : ’ELEMENTARY NUMERICAL ANALYSIS’ 

SIS 10840 

c 


S I S 108 50 


I FLAG = 1 

SIS 10860 

c 


S I S 1 0870 

c 

INITIALIZE W 

SIS 10880 

c 


S I S 108 90 


DO 5 I = 1 ,N 

SIS 10900 


DO 5 0 = 1 ,N 

S I S 1 09 10 

5 

W( I ,0 ) = A( I ,0 > 

S I S 1 0920 

c 


S I S 1 09 30 

c 

GAUSS ELIMINATION WITH SCALED PARTIAL PIVOTING 

S I S 1 0940 

c 


S I S 1 09 50 


NMIN1 = N - 1 

S I S 1 09 60 


DO 100 I = 1.NMIN1 

S IS 10970 


IF(CAES(W( I , I ) ) .LT.TOL ) GO TO 999 

S I S 10980 

c 


SIS 10990 


IPLU1 =1+1 

SIS' 1000 
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DO 80 0 = IPLU1 ,N 

SIS11010 


wc o , i > = wij.i) / wn.n 

SIS1 1020 


RATIO = -W(0,I> 

SIS1 1030 


DO 40 K = IPLU1 ,N 

SIS1 1040 


W(J.K) = W(J,K> + RATIO * U( I , K ) 

SIS1 1050 

40 

CONTINUE 

SIS1 1060 

80 

CONTINUE 

SIS1 1070 

100 

CONTINUE 

SIS1 1080 


RETURN 

SIS1 1090 

C 


SIS1 1 100 

999 

IFLAG = 2 

SIS1 1 1 10 


RETURN 

SIS1 1120 


END 

SIS1 1 130 


SUBROUTINE SUBST( W, B ,X , N ) 

SIS1 1 140 

c 


SIS1 1 150 


COMPLEX SUM 

SIS1 1 160 


COMPLEX W(N,N), B ( N ) , X ( N ) 

SIS1 1 170 

C 


SIS1 1 180 

C 

THIS SUBROUTINE SOLVES THE MATRIX EQUATION AX=B 

SIS1 1 190 

C 

WITH A IN LU FORM. SEE CONTE AND DE BOOR: 

S IS 1 1200 

C 

’ELEMENTARY NUMERICAL ANALYSIS’ 

SIS1 1210 

C 


SIS1 1220 


IF(N.GT.l) GO TO 10 

SIS1 1230 


X ( 1 > = B(l> / W< 1 , 1 ) 

SIS1 1240 


RETURN 

SIS1 1250 

C 


SIS1 1260 

10 

X( 1 ) = B(l) 

SIS1 1270 


DO 15 K = 2 , N 

SIS1 1280 


KM1 = K - 1 

SIS1 1290 


SUM = 0. 

SIS1 1300 


DO 14 J = 1 , KM1 

SIS1 1310 

14 

SUM = W( K , U ) *X ( 0 ) + SUM 

SIS1 1320 

15 

X ( K ) = B ( K ) - SUM 

SIS1 1330 

C 


SIS1 1340 


X ( N ) = X ( N ) / W(N,N) 

SIS1 1350 


K = N 

SIS1 1360 


DO 20 NP1MK = 2 , N 

SIS1 1370 


KP1 = K 

SIS1 1380 


K = K - 1 

SIS1 1390 


SUM = 0. 

SIS1 1400 


DO 19 0 = KPl.N 

SIS1 1410 

19 

SUM = W(K,0 ) * X<0) + SUM 

SIS1 1420 

20 

X ( K ) = ( X ( K ) - SUM) / W(K,K) 

SIS1 1430 

c 


SIS1 1440 


RETURN 

SIS1 1450 


END 

SIS1 1460 


SUBROUTINE F FT ( NUMBER , DATA , FACTOR ) 

SIS1 1470 


INTEGER CRPSZE,POSN,POINTR,ELMTNO,FACPTR,ABSNUM 

SIS1 1480 


INTEGER GRPTR 

SIS1 1490 


LOGICAL INVERT 

S IS 1 1500 


COMPLEX TEMP ,W, FACTOR ( NUMBER > , DATA! NUMBER ) .WSCALE , TEMP 1 , 

SIS1 1510 


1TEMP2 

SIS1 1520 


DATA P 1/3. 141 5926/ 

SIS1 1530 

C 


SIS1 1540 

C AUTHOR • R.G. HICKS, UNIV. OF OLD., APRIL 1977. 

SIS1 1550 
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C 

C 

C 

C 

C 

C 

C 


10 


C 


15 

C 


17 

C 

20 


25 


C 


30 

40 


50 

C 

60 

C 


S 

AIM: THIS SUBROUTINE PERFORMS A FAST FOURIER TRANSFORM ON THE S 
INPUT DATA USING THE COOL EY-TURKEY ALGORITHM S 

THERE MUST BE 2**NUMBER POINTS GIVEN AND THE GIVEN S 

AND THE GIVEN POINTS ARE DESTROYED BY THE SUBROUTINE S 

IF NUMBER IS NEGATIVE, AN INVERSE FOURIER TRANSFORM IS DONE.S 

S 

INVERT = .FALSE. S 

ABSNUK = NUMBER S 

I F ( NUMBER . GT .0 ) GO TO 10 S 

AESNUM = -NUMEER S 

INVERT = .TRUE. S 

N = 2 ** ABSNUM S 

ARG = 2. * PI / FLOAT ( N ) S 

TEMP = CMPLXI0. , ARG ) S 

W = CEXP ( TEMP ) S 

S 

FACTOR ( 1 ) = W S 

NODIV2 = N/2 S 

DO 15 INDEX = 2.NODIV2 S 

FACTOR! INDEX) = FACTOR! INDEX - 1 > *W S 

S 

IF! INVERT) GO TO 20 S 

DO 17 POSN = 1 ,N S 

DATA! POSN ) = CONOG ! DATA! POSN ) ) S 

S 

GRPSZE = N S 

NOGRP = 1 S 

ITERNO =0 S 

LIMIT = GRPSZE/2 S 

ITERNO = ITERNO + 1 S 

POINTR = GRPSZE/2 + 1 S 

POSN = 1 S 

S 

DO 60 GRPTR = 1 , NOGRP S 

DO 50 ELMTNO = 1 .LIMIT S 

I POSN = POSN - 1 S 

CALL FACNO! IPOSN.N, ITERNO, FACPTR) S 

IF ( FACPTR . EQ.0) GO TO 30 S 

WSCALE = FACTOR (FACPTR) S 

TEMPI = DATA! POINTR ) * WSCALE + DATA! POSN ) S 

WSCALE = -WSCALE S 

TEMP2 = DATA! PO I NTR > * WSCALE + DATA! POSN ) S 

GO TO 40 S 

TEMPI = DATA! POI NTR ) + DATA { POSN ) S 

TEMP2 = DATA! POSN ) - DATA! POI NTR > S 

DATA! POSN) = TEMPI S 

DATA! POI NTR ) = TEMP2 S 

POSN = POSN + 1 S 

POINTR = POINTR + 1 S 

S 

POSN = POSN + GRPSZE/2 S 

POINTR = POINTR + GRPSZE/2 S 

S 

s 


IS 

11560 

IS 

11570 

IS 

11580 

IS 

11590 

IS 

11600 

IS 

11610 

IS 

11620 

IS 

11630 

IS 

11640 

IS 

11650 

IS 

11660 

IS 

11670 

IS 

1 1680 

IS 

11696 

IS 

11700 

IS 

11710 

IS 

11720 

IS 

11730 

IS 

11740 

IS 

11750 

IS 

11760 

IS 

11770 

IS 

11780 

IS 

11790 

IS 

11800 

IS 

11810 

IS 

11820 

IS 

11836 

IS 

11840 

IS 

11850 

IS 

1 1860 

IS 

11870 

IS 

11880 

IS 

11890 

IS 

1 1900 

IS 

11910 

IS 

1 1920 

IS 

11930 

IS 

1 1940 

IS 

11950 

IS 

11960 

IS 

11970 

IS 

11980 

IS 

11990 

IS 

12000 

IS 

12010 

IS 

12020 

IS 

12030 

IS 

12040 

IS 

12050 

IS 

12060 

IS 

12070 

IS 

12080 

IS 

12090 

IS 

12100 


GRPSZE 


GRPSZE/2 
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NOGRP = NOGRP * 2 
IF ( GRPSZE .GT. 1 ) GO TO 25 
C 

DO 70 POSN = 1 , N 
IPOSN = POSN - 1 
CALL JUXTA1N, IPOSN, JUXPOS) 

JUXPOS = JUXPOS + 1 
IFIJUXPOS.LT. POSN) GO TO 70 
TEMP - DATA! POSN) 

DATA! POSN ) = DATAIJUXPOS) 

DATA( JUXPOS ) = TEMP 
70 CONTINUE 
C 

IF(INVERT) RETURN 

DO 80 POSN = 1 ,N 

TEMP = CONJG < DATA! POSN ) ) 

80 DATA(POSN) = TEMP / N 
C 

RETURN 

END 

SUBROUTINE J UXTAl N , I POSN , OUXPOS ) 

NN = N 
NUM = IPOSN 
OUXPOS = 0 
C 

DO 10 1=1 ,N 
NN = NN/2 
MEM = NUM 
NUM = NUM/ 2 
IREM = MEM - 2*NUM 
OUXPOS = 2*JUXP0S + IREM 
IFINN.EQ. 1 ) GO TO 20 
10 CONTINUE 
20 RETURN 
END 

SUBROUTINE FACNOII POSN , N , ITERNO , FACPTR ) 

INTEGER FACPTR 

THIS ROUTINE IS PART OF THE BUTTERFLY CALCULATION 
PROCEDURE IN PERFORMING A FAST FOURIER TRANSFORM 

FACPTR = 0 

I D I V = N/ ( 2** ITERNO) 

ITEMP = IPOSN 
ITEMP = ITEMP / I D I V 
MUL = N/2 
DO 10 1=1, ITERNO 
IP EM = ITEMP - ITEMP /2*2 
FACPTR = IREM*MUL + FACPTR 
ITEMP = ITEMP / 2 
MUL = MUL/2 
10 CONTINUE 
RETURN 
END 

SUBROUTINE ZERO( POI NTS , NUMBER ) 


SIS121 10 
SIS12120 
S IS 1 2 1 30 
SIS12140 
SIS12150 
SIS12160 
SIS12170 
SIS12180 
SIS12190 
SIS12200 
SIS12210 
S IS 1 2220 
S IS 12230 
S I S 1 2240 
S I S 1 2250 
SIS122G0 
S IS 1 2270 
S IS 1 2280 
S I S 1 2290 
S I S 1 2300 
S IS 1 23 10 
S I S 1 2320 
S I S 1 233£ 
S IS 1 2340 
S IS 1 235C 
S I S 1 2360 
S IS 1 2370 
S I S 1 2380 
S I S 1 23 90 
S IS 1 2400 
S IS 1 24 10 
S IS 1 2420 
S IS 1 2430 
S I S 1 24 40 
S I S 1 2450 
S I S 1 2460 
SIS12470 
S I S 1 2480 
S I S 1 24 90 
SIS 12500 
SIS12510 
S I S 1 2520 
S I S 1 2530 
S IS 1 2540 
S I S 1 2550 
S I S 1 2560 
S I S 1 2570 
S I S 1 2580 
S I S 1 2590 
SIS 12600 
SIS12610 
S I S 1 2620 
S I S 1 2630 
S IS 1 2640 
S I S 1 26 50 



n o n o r> o o •— nm oooooooo 


65.23 


FILE: SISCAP 


FORTRAN A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


COMP 

REAL 

DATA 

LEX POI NTS ( NUMBER ) 
MAX 

TOL/1 .E-5/ 


AUTHOR: 

R.G. HICKS, UNIV. OF OLD., 

APRIL 1977. 

AIM: 

THIS SUBROUTINE AUTOMATICALLY ZEROS ANY 
QUANTITI I ES WHICH ARE LESS THAN TOL 
TIMES THE LARGEST QUANTITY IN THE SECTION. 

MAX 

= CABS ( PO I NTS ( 1 > ) 



DO 5 J=2, NUMBER 

IFIMAX.LT. CABSIPOINTSfU > > > MAX = CABS ( POI NTS {0 ) ) 

CONTINUE 

DO 10 0 = 1 .NUMBER 

IF (CABS( POINTSU ) ) .LT.TOL*MAX) POINTStO) = 0.0 
CONTINUE 

RETURN 
END 

SUBROUTINE THEORY ( VLO , VB IAS , NMAX , I VDAT , NOI V , J , FREO , YS IG , 
INCISE. TEMP , DEL TAV , RN , Y0, Y1 .Y1M.ZLO) 

REAL I VDAT { NOI V ) , VALU1) 

COMPLEX J(NOIV), YS I G { 3 , 3 ) , NOISE<3,3>, CVAL 1 , CVAL2 
COMPLEX YjC, Y1.Y1M, ZLO, YLO 

COMMON /PEOUT/ VX , VY , VJDC, VOLO, CURDC, CURPMP, RLOSS, RTEMP 
DATA E.P I .HBAP/l .G02E-19, 3. 1415926,1 .055E-34/ 

DATA EOLTZ, TOL / 1 . 380622E-23 , 1 . E-5/ 

THIS SUBROUTINE SETS UP THE VALUES FOR CALCULATION 
OF THE SIGNAL CONVERSION AND NOISE TEMPERATURE FOR THE 
THREE FREQUENCY APPROXIMATION 

OMEGA = 2. * PI * FREQ 
ALPHA = E * VLO / HBAR / OMEGA 
C 

WRITE(G.lll) VB IAS , VLO, ALPHA 
111 F ORMAT ( IX. ’ VB I AS = ’ , E 1 5 . 8 , 1 X , ’ VLO = ’.E15.8, 

1 IX, ’ALPHA = ’ ,E15.8) 

C 

NMPL4 = NMAX + 10 
NR EQ = NMFL4/2 + 1 
C 

DO 100 I = l.NREQ 
N = NR.EQ - I 

CALL BESOfALPh'A.N, VAL { I ) ,TOL , IER ) 

IF { IER.GT.0) GO TO 150 
ICORR = NMPL4 - I + 1 
VAL < ICORR ) = VAL ( I ) 

IF (N/2*2.EQ.N) GO TO 100 
VAL (ICORR) = -VAL(I) 

100 CONTINUE 


S 

IS 

12660 

S 

IS 

12G70 

S 

IS 

12G80 

S 

IS 

12G90 

S 

IS 

12700 

S 

IS 

12710 

S 

IS 

12720 

S 

IS 

12730 

S 

IS 

12740 

S 

IS 

12750 

S 

IS 

127G0 

s 

IS 

12770 

s 

IS 

12780 

s 

IS 

12790 

s 

IS 

12800 

s 

IS 

12810 

s 

IS 

12820 

s 

IS 

12830 

s 

IS 

12840 

s 

IS 

12850 

s 

IS 

128G0 

s 

IS 

12870 

s 

IS 

12880 

s 

IS 

12890 

s 

IS 

12900 

s 

IS 

12910 

s 

IS 

12920 

s 

IS 

12930 

s 

IS 

12940 

s 

IS 

12950 

s 

IS 

129G0 

s 

IS 

12970 

s 

IS 

12980 

s 

IS 

12990 

s 

IS 

13000 

s 

IS 

13010 

s 

IS 

13020 

s 

IS 

13030 

s 

IS 

13040 

s 

IS 

13050 

s 

IS 

13060 

s 

IS 

1 3070 

s 

IS 

13080 

s 

IS 

13090 

s 

IS 

13100 

s 

IS 

13110 

s 

IS 

13120 

s 

IS 

13130 

s 

IS 

13140 

s 

IS 

13150 

s 

IS 

13160 

s 

IS 

13170 

s 

IS 

13180 

s 

IS 

13190 

s 

IS 

13200 
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c 


SIS13210 


GO TO 200 

SIS 13220 

c 


S IS 1 3230 

150 

WRITE (6, 160) 

S I S 1 32 40 

1G0 

FORMAT! IX. ’PROBLEM IN EVALUATING THE BESSEL FUNCTION’) 

S I S 1 3250 


STOP 

S IS 1 3250 

C 


S IS 1 3270 

200 

CGNST = HBAR * OMEGA / E 

S I S 1 3280 

c 


S IS 1 3290 


GSUM00 = 0. 

S IS 1 3300 


GSUM10 = 0. 

SIS13310 


GSUM01 = 0. 

S I S 1 3320 


GSUM11 = 0. 

S I S 1 3 3 30 


GSUM1M = 0. 

S IS 1 3340 

c 


S I S 1 3 350 


BSUM00 = 0. 

S IS 1 3360 


BSUM10 = 0. 

S IS 1 3370 


BSUM01 = 0. 

SIS13380 


BSUM11 = 0. 

S I S 1 3390 


BSUM1M = 0. 

S IS 1 3400 

c 


SIS13410 


H00 = 0. 

S I S 1 3 420 


H 10 = 0. 

S I S 1 3430 


H01 = 0. 

SIS13440 


HI 1 = 0. 

S I S 1 3450 


HIM = 0. 

S I S 1 3460 

c 


SIS13470 


ARG = VBIAS 

S I S 1 3480 


ARG 1 = ARG * E / HBAR 

S I S 1 3490 


CALL EVAL0IARG1 .CVALl ,0 , NOI V , DELTAV.RN > 

S IS 1 3500 


CURDC = AIMAG(CVAL 1 ) 

SIS13510 

c 


S IS 1 3520 


CURPMP = 0.0 

S IS 1 3530 


CURLOl = 0.0 

S IS 1 3540 


CURL02 = 0.0 

S I S 1 3550 


CURL03 =0.0 

S 1 S 1 3560 


CURL04 = 0.0 

S IS 1 3570 


CURL05 = 0.0 

S I S 1 3580 


CURL06 = 0.0 

S I S 1 3590 


CURL07 = 0.0 

SIS 13600 


CURL08 = 0.0 

S I S 1 3 6 1 0 


CURL09 = 0.0 

SIS1 3620 


CURLO0 = 0.0 

S I S 1 3630 

c 


S I S 1 3 640 


DO 300 I = l.NMAX 

S I S 1 3650 

c 


SIS1 3660 


N = NHAX/2 - I + 1 

SIS13670 


ARG = VBIAS + FLOAT! N ) * CONST 

S I S 1 3680 


ARG 1 = ARG * E / HBAR 

S I S 1 36 90 


CALL EVALO (ARG1 , C VAL 1 , 0 , NO I V , D ELTAV , RN > 

S I S 1 3700 


CALL DERIJ(ARG1 , CVAL 2 . 0 , NO I V . DELTAV . RN ) 

SIS13710 


TMP 1 = A I MAG (CVALl ) 

S I S 1 3720 


TMP2 = REAL (CVALl) 

S I S 1 37 30 


TMP 3 = A I M.AG ( CVAL 2 ) 

S I S 1 37 40 


TMP 4 = REAL ( CVAL2 ) 

S I S 1 3750 
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C 

C 


C 


C 


c 


c 


300 

c 


c 


NPTR = NREQ - N 


CURPMP = 
CURLOl = 
CURL02 = 
CURL03 = 
CURL04 = 
CURL05 = 
CURL06 = 
CURL07 = 
CURL08 = 
CURL09 = 
CURLO0 = 


VAL { NPTR ) **2 * TMP 1 + 
VAL ( NPTR ) * ( VAL ( NPTR- 1 
VAL ( NPTR )*( VAL (NPTR-1 
VAL ( NPTR ) * ( VAL ( NPTR -2 
VAL < NPTR ) * { VAL ( NPTR -2 
VAL (NPTR)*( VAL ( NPTR -3 
VAL ( NPTR )*( VAL (NPTR -3 
VAL { NPTR ) * ( VAL ( NPTR -4 
VAL ( NPTR >*( VAL ( NPTR -4 
VAL { NPTR ) * ( VAL ( NPTR-5 
VAL (NPTR )*( VAL (NPTR-5 


CURPMP 

) +VAL ( NPTR+ 1 > >*TMP1 
) -VAL (NPTR+1 > )*TMP2 
) +VAL ( NPTR+2 ) ) *TMP 1 
) - VAL (NPTR + 2) )*TMP2 
) + VAL (NPTR + 3))*TMP1 
) - VAL (NPTR + 3)) *TMP2 
)+VAL(NPTR+4))*TMPl 
) -VAL ( NPTR + 4 ) >*TMP2 
) +VAL (NPTR + 5))*TMP1 
> -VAL (NPTR+5) )*TMP2 


+ CURLOl 
+ CURL02 
+ CURL03 
+ CURL04 
+ CURL05 
+ CURL06 
+ CURL07 
+ CURL08 
+ CURL09 
+ CURLO0 


GSUM00 
GSUM10 
GSUM01 
GSUM1 1 
GSUM1M 


VAL ( NPTR ) **2 * TMP3 + GSUM00 

VAL ( NPTR ) * (VAL( NPTR + 1 >+VAL (NPTR-1 ) >*TMP3 + GSUM10 
VAL ( NPTR )*( VAL (NPTR+1 ) -VAL ( NPTR- 1 >> *TMP 1 +GSUM01 
( (VAL(NPTR + 1 ) )**2 - ( VAL (NPTR-1 ))**2) *TMP 1 + GSUM1 1 
VAL(NPTR)*(VAL(NPTR+2)-VAL(NPTR-2)>*TMPl + GSUM1M 


BSUM10 = VAL ( NPTR ) * ( VAL ( NPTR+ 1 )-VAL(NPTR-l ) ) * TMP 4 + BSUM10 
BSUM11 = < (VAL(NPTR+1 > >**2 - 2 . * < VAL ( NPTR ) ) **2 + 

1 (VAL(NPTR-1 > >**2) * TMP2 + BSUM11 
BSUM1M = ( VAL (NPTR + 2 ) * VAL ( NPTR ) -2 . *VAL (NPTR+1 > *VAL (NPTR-1 ) + 
1 VAL ( NPTR )* VAL ( NPTR- 2 ) ) *TMP2 + BSUM1M 


APG = ARG * E / 2. / BOLTZ / TEMP 
STO = COTH(ARG) 


H00 = ( VAL ( NPTR ) ) **2 * STO * TMP1 + H00 

H 10 = VAL ( NPTR > * ( VAL (NPTR + 1 ) + VAL (NPTR-1 ) > * STO * TMP 1 
1 + H 10 

H 1 1 = ( (VAL(NPTR+1 ) )**2 + ( VAL ( NPTR- 1 >) **2 ) * STO * TMP 1 + Hll 

HIM = VAL (NPTR+1 ) * VAL(NPTR-l) * STO * TMP1 + HIM 

CONTINUE 


GSUM10 

= 

0.5 * i 

GSUM10 




GSUM01 

= 

GSUM01 

★ 

E 

/ 

HBAR 

/ 

OMEGA 

GSUM1 1 

= 

GSUM1 1 

* 

E 

/ 

HBAR 

/ 

OMEGA 

GSUM1M 

= 

GSUM1M 

* 

E 

/ 

HBAR 

/ 

OMEGA 

BSUM10 

= 

BSUM10 

★ 

0. 

5 




BSUM1 1 

= 

BSUM1 1 

* 

E 

/ 

HBAR 

/ 

OMEGA 

BSUM1M 

= 

ESUM1M 

* 

E 

/ 

HBAR 

/ 

OMEGA 


H £0 = H00 * 2. * E 
H 1 0 = H 1 0 * E 
Hll = Hll * E 
HIM = HIM * 2. * E 


* 0.5 

* 0.5 

/ 2 . 

/ 2 . 


GSOURC 

BSOURC 

GIMAGE 

BIMAGE 

GLOAD 

BLOAD 


= R EAL ( Y 1 ) 

= A I MAG ( Y1 ) 

= PEAL(YIM) 

= 4 IMAG(YIM) 
= RE AL ( Y0) 

= A I MAG ( Y0) 


SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 

SIS 


137S0 
13770 
13780 
13790 
13800 
13810 
13820 
13830 
13840 
13850 
13860 
13870 
13880 
13890 
13900 
13910 
13920 
13930 
13940 
13950 
13960 
13970 
13980 
13990 
14000 
14010 
14020 
14030 
14040 
14050 
1 4060 
14070 
14080 
1 4090 
14100 
14 110 
14120 
14130 
14140 
14150 
14160 
14 170 
14180 
14190 
14200 
14210 
14220 
14230 
1 4240 
14250 
14260 
14270 
14280 
14290 
14300 
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YLO = CMPLX! 1 . 0 , 0 . 0 ) / ZLO 
GPUMP = REAL { YLO ) 

BPUMP = AIMAG ( YLO ) 

PLO = ( !CURL01+GPUMP*VL0}**2+!CURL02-BPUMP*VL0)**2>/8. /GPUMP 
WR I TE ( 6 , 400 ) PLO 

400 F0RMAT!/1X, ’LOCAL OSCILLATOR POWER = ’.E15.8,’ WATTS PER 0N’/> 

C 

CALL GEN3FT(GSUM00,GSUM10,GSUM01 ,GSUM1 1 ,GSUM1M, BSUM10. 

1BSUM11 , BSUM 1 M , H00 , H 10 , H 1 1 . H 1 M , GSOURC , GLOAD , BSOURC . G I MAGE , B I MAGE , 
1BLOAD) 

C 

YSIGU.l) = CMP LX( GSUN1 1 , BSUM 11 ) 

YSIG(1,2> = CMPLX!GSUM10,BSUM10) 

YSIG<1,3) = CMPLX<GSUM1M,BSUM1M) 

YS I G ( 2 , 1 ) = CMPLX! GSUM01 , 0 . ) 

YS I G ( 2 , 2 ) = CMPLX(GSUM00,0. > 

YS I G ( 2 , 3 ) = CMPLX ( GSUK0I ,0. > 

YS I G ( 3 , 1 ) = CMPLX!GSUM1M,-BSUM1M> 

YS I G ( 3 , 2 > = CKPLX(GSUK10,-BSUM10> 

YS I G ( 3 , 3 ) = CMPLXIGSUM1 1 ,-BSUKl 1 ) 

C WR ITE ! 6 , 2832 ) ( ( MK, ML , YSIG ( MK , ML > ,ML = 1 , 3 ) ,MK= 1,3) 

C2832 FORMAT! IX, ’MK= ’ ,I5,’ML=',I5,’YSIG=’,2E15.8) 

C 

NOISE! 1 , 1 > = CMPLX!H11 ,0. > 

NOISE! 1 ,2) = CMPLX! H10.0. ) 

NOISE! 1 , 3 ) = CMPLX! HIM, 0. > 

NO I SE ( 2 , 1 ) = CMPLX ( H 10 ,0. ) 

NO I SE ! 2 , 2 ) = CMPLX ! H00 , 0. ) 

NO I SE ! 2 , 3 ) = CMPLX(H10,0. > 

NOISE(3,l> = CMPLX ! H 1M ,0 . ) 

NO I SE ( 3 , 2 ) = CMPLX(H10,0. ) 

NO I SE < 3 , 3 > = CMPLX(H1 1 ,0. ) 

C 

RETURN 
END 

SUBROUTINE DER 1 0 ! DUM , CVAL , 0 , NO I V , DELTAV , RN > 

COMPLEX CVAL, J(NOIV>, TERP5 
LOGICAL FLAG 

DATA E.HBAR/l .602E-19, 1 .055E-34/ 

THIS SUBROUTINE EVALUATES THE DERIVATIVE 
OF J BY INTERPOLATING AMONGST A VECTOR OF 0 VALUES 

CONST = 1 . / RN 
FLAG = .FALSE. 

ARG = DUM 

IFiARG) 100,400,300 
100 FLAG = .TRUE. 

ARG = -DUM 

300 PTR = ARG / DELTAV / E * HEAR 
I PTR = PTP. + 1. 

IF! IPTR.GE.NOIV-1 ) GO TO 200 
RATIO = At OD( PTR , 1 . ) 

CVAL = TER P5 ! 0 ,N01V, IPTR .RATIO, .TRUE . > /DELTAV + CMPLX ! 0 ., CONST ) 

I F ( FLAG ) CVAL = -1. * CONOG(CVAL) 


SIS14310 
SIS14320 
S I S 1 4 330 
S IS 1 4340 
S I S 1 4350 
SIS14360 
S IS 1 4370 
S IS 1 4380 
S IS 1 4390 
S IS 1 4400 
S IS 1 44 10 
S IS 1 4420 
S IS 1 4430 
S IS 1 4440 
S IS 1 4450 
S IS 1 4460 
S IS 1 4470 
SIS14480 
S I S 1 4 490 
SIS 14500 
SIS14510 
S I S 1 4520 
S IS 1 4530 
S I S 1 4 540 
S IS 1 4550 
S IS 1 4560 
S IS 1 4570 
S IS 1 4580 
S IS 1 4590 
SIS 1 4500 
S IS 1 46 10 
S IS 1 4620 
S IS 1 4630 
S I S 1 4 640 
S I S 1 4 650 
S I S 1 4660 
SIS14670 
S I S 1 4 680 
S I S 1 4690 
S IS 1 4700 
SIS14710 
S I S 1 4 720 
S I S 1 4730 
S I S 1 4740 
S I S 1 4750 
S I S 1 4760 
S I S 1 4 770 
S IS 1 4780 
S I S 1 4790 
S I S 1 4800 
S I S 1 48 10 
S I S 1 4 820 
S I S 1 4 830 
SIS14840 
S I S 1 4850 
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RETURN 




SIS 14860 

200 

CVAL = CMPLXI0. .CONST) 




S IS 1 4870 


IF ( FLAG ) CVAL = -I. * CONOG(CVAL) 




S IS 1 4880 


RETURN 




S I S 1 4890 

400 

CVAL = CMPLXI0. ,0. > 




S I S 1 4 900 


RETURN 




SIS14910 


END 




SIS 14920 

C 





S I S 1 4930 

c 





SIS 14940 

c 





S I S 1 4950 

c 

SUBROUTINE BESO 




S I S 1 4960 

c 





SIS14970 

c 

PURPOSE 




S I S 1 4 980 

c 

COMPUTE THE J BESSEL FUNCTION FOR A GIVEN ARGUMENT 

AND 

ORDES I S 1 4 990 

c 





SIS 15000 

c 





S I S 1 50 1 0 

c 

USAGE 




S I S 1 5020 

c 

CALL BESU(X,N,BU,D, IER) 




S I S 1 5030 

c 





S I S 1 5040 

c 

DESCRIPTION OF PARAMETERS 




S I S 1 5050 

c 

X -THE ARGUMENT OF THE J BESSEL FUNCTION DESIRED 


SIS15060 

c 

N -THE ORDER OF THE 0 BESSEL FUNCTION 

DESIRED 



S I S 1 5070 

c 

BJ -THE RESULTANT 0 BESSEL FUNCTION 




S I S 1 5080 

c 

D -REQUIRED ACCURACY (FRACTIONAL) 




S I S 1 5090 

c 

I ER.-RESULTANT ERROR CODE WHERE 




SIS 15 100 

c 

I ER =0 N0 ERROR 




SIS151 10 

c 

I ER = 1 N IS NEGATIVE 




S I S 1 5 1 20 

c 

I ER = 2 X IS NEGATIVE OR ZERO 




S I S 1 5 1 30 

c 

I ER = 3 REQUIRED ACCURACY NOT OBTAINED 




S I S 1 5 1 40 

c 

I ER = 4 RANGE OF N COMPARED TO X NOT CORRECT (SEE 

REMARKS SIS15150 

c 





S I S 1 5 1 60 

c 

REMARKS 




SIS 1 51 70 

c 

N MUST BE GREATER THAN OR EQUAL TO ZERO 

, BUT IT 

MUST BE 

S I S 1 5 1 80 

c 

LESS THAN 




S I S 1 5 1 90 

c 

20 +10*X-X** 2/3 FOR X LESS THAN OR 

EQUAL TO 

15 


SIS 15200 

c 

90+X/2 FOR X GREATER THAN 

15 



SIS15210 

c 





S I S 1 5220 

c 

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 



S I S 1 5230 

c 

NONE 




SIS15240 

c 





S IS 1 5250 

c 

METHOD 




SIS1 5260 

c 

RECURRENCE RELATION TECHNIQUE 




S IS 1 5270 

c 





SIS15280 

c 





SIS1 5290 


SUBROUTINE BESU(XPPP , N , BJPPP , DPPP , IER > 




SIS 15300 


IMPLICIT REAL*8 (A-H.O-Z) 




SIS15310 


REAL BJPPP 




SIS1 5320 


REAL XPPP, DP P P 




S IS 1 5330 


X = XPPP 




S IS 1 5340 


D = DP P P 




S IS 1 5350 

c 





S I S 1 5360 


BJPPP = 0.0E0 




S IS 1 5370 


BO =0 .0D0 




S IS 1 5380 


I F ( N ) 1 0 , 20 > 20 




S I S 1 5 390 

10 

I ER = 1 




S IS 1 5400 
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RETURN 

20 I F ( X > 30, 300, 3 1 

30 I ER = 2 
RETURN 

300 BOPPP = 0.0E0 

I F ( N . EQ . 0 ) BOPPP = 1.0E0 

IER = 0 

RETURN 

31 IF (X-15.D0>32,32,34 

32 NTEST=20.D0+10.D0*X-X** 2.O0/3.D0 
GOTO 36 

34 NTEST=90. D0+X/2.D0 

36 IF ( N-NTEST ) 40, 38 , 38 

38 I ER = 4 

RETURN 
40 I ER =0 

N 1 =N+ 1 
BPREV=0.0D0 
C 

C COMPUTE STARTING VALUE OF M 

C 

IF ( X-5 . D0 ) 50, 60, 60 
50 MA=X+6.D0 

GOTO 70 

60 MA= 1 . 4D0*X+60. D0/X 
70 MB=N+IDINT(X)/4+2 
MZERO=MAX0(MA,MB ) 

C 

C SET UPPEP LIMIT OF M 

C 

MMAX = F TEST 

100 DO 190 M=MZERO , MMAX , 3 
C 

C SET F ( M ) , F ( M- 1 ) 

C 

FN 1 = 1 . 0D-28 
FM=0.0D0 
ALP HA=0 . 0D0 

I F(M-< M/2 >*2)120, 110, 120 
110 UT = - 1 

GO TO 130 
120 0T= 1 

130 M2=M-2 

DO 160 K= 1 , M2 
MK=M-K 

BK,K = 2 . D0* D FLOAT (MK)*(FM1/X)-FM 

FM=FM1 

FK1 =BMK 

IF (MK-N-1 >150,140,150 
140 EO = BMK 

150 UT=-0T 
S= 1 +OT 

160 ALPHA~ALPHA+BMK*S 
BMK = 2 . 0D0* FM1 /X-FM 
IF (N)1C0, 170, 180 


S IS 1 54 10 
SIS 15420 
S I S 1 5430 
S I S 1 5440 
S I S 1 5450 
S IS 1 5460 
S I S 1 5470 
S I S 1 5480 
S I S 1 5490 
S I S 1 5500 
SIS15510 
SIS15520 
S I S 1 5530 
S IS 1 5540 
S I S 1 5550 
S I S 1 5560 
SIS15570 
S I S 1 5580 
S IS 1 5590 
S I S 1 5600 
SIS15610 
S I S 1 5620 
SIS15630 
SIS15640 
S I S 1 5658 
S I S 1 5660 
SIS15670 
S I S 1 5680 
S I S 1 5690 
S I S 1 5700 
SIS15710 
S I S 1 5720 
S I S 1 5730 
S I S 1 5740 
SIS15750 
S I S 1 5768 
SIS15770 
S I S 1 5780 
S I S 1 5790 
S I S 1 5808 
SIS15810 
SIS15820 
S I S 1 5830 
SIS15840 
S IS 1 5850 
SIS 15860 
S I S 1 5870 
S I S 1 5880 
S I S 1 5890 
SIS 15900 
S I S 1 59 10 
S I S 1 59 28 
S IS 1 5930 
S IS 1 5940 
S IS 1 5950 
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170 B0=BMK 
180 ALPHA=ALPHA+BMK 
BU=BO/ALPHA 

IF (DABS(BJ-BPREV) -DABS ( D*B0 ) >200,200, 190 
190 BPREV=B0 
I ER = 3 

200 BJPPP = BO 
RETURN 
END 

FUNCTION COTH1ARG > 

REAL NUM, NEG 
NEG = -ARG 

NUM = EXP ( ARG ) + EXP ( NEG ) 

DEN = EXP ( ARG > - EXP ( NEG ) 

IF (ABS(DEN) .LT. 1 .E-5) GO TO 10 

COTH = NUM / DEN 

RETURN 

10 COTH = 1.E10 
RETURN 
END 

SUBROUTINE GEN3PT ( GSUM00 , GSUM 10 , GSUM01 , GSUM 1 1 , GSUM 1 M , 
1BSUM10.BSUM1 1 , BSUM 1 M , H00 , H 10, H 1 1 , H 1 M , GSOURC , GLOAD , BSOURC , 
2GIMAGE, BIMAGE, BLOAD) 

COMMON /RBOUT/ VDC , VLO, VODC, VOLO, CURDC, CURPMP , LOSS, CREAL 
C 

REAL L0 , LOSS 
COMPLEX LAM01 , CTM 
COMPLEX Z01 , LAM01M 

DATA TOL/l.E-5/, BOLTZ/1 .380622E-23/ 

THIS SUBROUTINE CALCULATES THE SMALL SIGNAL AND 
NOISE PAR/METERS USING THE VALUES THAT HAVE BEEN 
SET UP USING SUBROUTINE ’THEORY’ 

L0 = 2. * GSUM 10 / GSUM01 

CSI = (GSUM11 - GSUM1K)/(GSUM1 1 + GSUM1M) 

G4MMA = BSUM1M / (GSUKll + GSUM1M) 

ETA = 2. * GSUM01 * GSUM10 / GSUM00 / (GSUM11 + GSUM1M) 

GL = GLOAD / GSUM00 
BETA = BSUM10 / GSUM10 
B 1 = BIMAGE / (GSUM11 + GSUM1K) 

GA = (GSOURC + GIMAGE) / 2.0 / (GSUM11 + GSUM1M) 

GD = (GSOURC - GIMAGE) / 2.0 / (GSUM11 + GSUM1M) 

BA = (BSUM11 + (BSOURC + BIMAGE) / 2.0) / (GSUM11 + GSUM1M) 

BD = (BSOURC - BIMAGE) / 2.0 / (GSUM11 + GSUM1M) 

BL = BLOAC / GSUM00 
C 

TERM1 = ((CSI + GA) * (1.0+ GA) + (BA**2 - GAMMA**2 ) - 
1 GD * *2 - B D * *2 > * (1.0 + GL + BL) 

TERM2 = (ED * (1.0 + CSI + 2.0*GA> - 2.0*GD*BA> * 

1(1.0 + GL + BL ) 

TERM3 = -1.0 * ETA * (CSI + GA + EETA*BA - BETA+GAMMA) 

TERM4 = -1.0 * ETA » (BD - BETA * GD) 

TERM5 = GSUM00 * (GSUM11 + GSUM1M) **2 

TERMG = -1.0 * GSUM01 * ( GSUM 11 +GSUM 1 M ) * (CSI + GA - GD) 


SIS15960 
SIS 15970 
S IS 1 5980 
S I S 1 5990 
SIS 16000 
S IS 1 6010 
S I S 1 6020 
S I S 1 6030 
S I S 1 6040 
SIS 16050 
S I S 1 6060 
S I S 1 6070 
SIS 16080 
S I S 1 6090 
SIS 16 100 
SIS161 10 
S I S 1 6 1 20 
S I S 1 6 1 30 
S I S 1 6 1 40 
S I S 1 6 1 50 
SIS 161 60 
SIS16170 
S I S 1 6 1 80 
S I S 1 6 1 90 
SIS 16200 
SIS16210 
S I S 1 6220 
S I S 1 6230 
S IS 1 6240 
SIS 16250 
S IS 1 6260 
S I S 1 6270 
S I S 1 6280 
S I S 1 62 90 
SIS 16300 
SIS16310 
S I S 1 6320 
S I S 1 6330 
S I S 1 6 340 
S I S 1 6350 
S I S 1 6360 
S I S 1 6 370 
S I S 1 6380 
S I S 1 6390 
SIS 16400 
SIS16410 
S I S 1 6420 
S IS 1 6430 
S I S 1 6 4 40 
S I S 1 6450 
S I S 1 64 60 
SIS16470 
S I S 1 64 80 
S I S 1 6490 
SIS 16500 
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TERM7 = -1.0 * GSUM01 * (GSUM1 1+GSUM1M) * (GAMMA + BD - BA) 

ZR1 = TERM5 * (TERM1 + TERM3 ) 

ZI1 = TERM5 * ( TERM2 + TERM4 ) 

ZR2 = TERM6 
Z 1 2 = TERM7 

Z01 = CMP LX ( ZR2 , Z I 2 ) / CMP LX ( ZR 1 . Z 1 1 ) 

LOSS = 1 . /4 . /GSOURC/G LOAD /CABS ( Z01 ) **2 
LOSS = 10. * ALOG 1 0( LOSS ) 

C 

TERM1 = -GSUM01 * ( GSUM1 1 +GSUM1M > * (CSI+GA-GD) 

TERM2 = -GSUM01 * ( GSUM1 1 +GSUM1M ) * ( GAMMA+BD-BA ) 

TERM3 = (GSUM11 + GSUM1M) ** 2 

TERM4 = (CSI+GA) * (1.0+GA) + BA**2 - GAMMA* *2 - GD**2 - BD**2 
TERM5 = BD * ( 1 ,0+CS I +2 .0*GA ) - 2.0*GD*BA 
TERM6 = TERM3 * TERM4 

TERM7 = TERM3 * TERMS 

LAM01 = CMPLX(TERM1 ,TERM2) / CMPLX (TERMS , TERM7 ) 

TERM8 = -GSUM01 * ( GSUM1 1 +GSUM1M ) * (CSI+GA+GD) 

TERM9 = -GSUM01 * ( GSUM11 +GSUM1M ) * ( BA+B D-GAMMA ) 

LAM01M = CMPLX(TERM8,TERM9) / CMPLXI TERM6 , TERM7 ) 

CTM = ( CAES ( LAM01 M ) **2 + CABS ( LAM01 > ** 2 ) * HI! 

1 * ( LAM01N+CONOG ( LAM01M ) +CONOG ( LAM01 ) +LAM0 1 ) * H10 

2 + H00 + (LAM01M*CONJG(LAM01 > +LAM0 1 *CONOG ( LAM01 M ) )*H1M 
CTM = 1./4./GSOURC/CABSI LAM01 )**2*CTM/BOLTZ 

CREAL = REAL < CTM ) 

Cl MAG = A I MAG (CTM) 

I F ( C I MAG .GT.TOL*CREAL ) GO TO 1000 
C 

WR I TE ( 6 , 80 ) LOSS, CREAL 

80 FORMAT! IX, ’CONVERSION LOSS (USB) = ’.E15.8,’ DB ’/ 

1 1 X , ’ SSB SIS MIXER NOISE TEMPERATURE (USB) = ’,E15.8,’ K ’ ) 
RETURN 
C 

1000 WRITE<6, 1001 ) 

1001 FORMAT* IX, ’PROBLEM DEVELOPED IN TEMPERATURE COMPUTATION’) 

STOP 
END 

COMPLEX FUNCTION TERP 5 ( J , NO I V , I PTR , RAT 1 0 , FLAG ) 

COMPLEX O(NOIV) 

COMPLEX STORE ( 4 ) 

LOGICAL FLAG 
DATA (1/3/ 

THIS COMPLEX VALUED FUNCTION PERFORMS AN INTERPOLATION 
AMONGST THE ARRAY OF 0 VALUES. 

IF FLAG IS TRUE, AN INTERPOLATION OF DERIVATIVES IS DONE 
IF FLAG IS FALSE, AN INTERPOLATION OF THE FUNCTION 
J IS DONE. 

THE INTERPOLATING FUNCTION USED IS A+BX+CEXP ( DX > 

AND THE ARRAY OF J VALUES MUST BE EVENLY SPACED 


NP1 = N + 1 

MID = (N + 1) / 2 


DO 1 K = 
KK = I PTR 



SIS16510 
S I S 1 6520 
S I S 1 6530 
S I S 1 6540 
S IS 1 6550 
S IS 1 6560 
SIS16570 
S I S 1 6580 
S I S 1 6590 
SIS16600 
SIS1 6610 
S I S 1 6620 
S I S 1 6630 
S IS 1 6640 
S IS 1 6650 
SIS16660 
S I S 1 6670 
S I S 1 6680 
S I S 1 6 6 90 
S I S 1 6700 
SIS16710 
S I S 1 6720 
S I S 1 6730 
S I S 1 6740 
S I S 1 6750 
S I S 1 6760 
S I S 1 6770 
S I S 1 6780 
S I S 1 6790 
S I S 1 6800 
SIS16810 
S I S 1 6820 
S I S 1 6830 
S I S 1 6840 
S I S 1 6850 
SIS! 6860 
SIS16870 
S IS 1 6880 
S I S 1 6890 
SIS 16900 
SIS16910 
SIS 1 6920 
S IS 1 6930 
S IS 1 6940 
S I S 1 6950 
S I S 1 6960 
S I S 1 6970 
S I S 1 6980 
S IS ! 6990 
S I S 1 7000 
S I S 1 70 1 0 
SIS 17020 
SIS 17030 
S IS 1 7040 
S I S 1 7050 


1 ) + ( K- 1 ) 
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C 

C 


c 

4 70 
C 

400 


C 


C 

570 

C 

500 


C 
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IFIKK.LE.0) GO TO 2 
STORE ( K ) = 0 < KK ) 

GO TO 1 
KK = 2 - KK 

STORE ( K ) = CONOG ( 0 ( KK > > 

CONTINUE 

YM1 = REAL (STORE! 1 ) ) 

Y0 = REAL ( STOR E < 2 ) ) 

Y 1 = REAL ( STORE ( 3 ) > 

Y2 = REAL ( STOP E ( 4 > ) 

CURV1 = Y2 - 2.0*Y1 + Y0 

CURV2 = Y1 - 2.0*Y0 + YM1 

IF ( CUP VI *CURV2 .LT. 0.0) GO TO 470 

I F ( ABS ( CUP VI ) .LE 0.0 .OR. ABS ( CURV2 ) . LE .0.0) GO TO 470 
D = ALOG(CURVl/CURV2) 

I F ( D . EQ . 0) GO TO 470 
C = CURV1 / ( ( EXP ( D ) - 1.0) ** 2) 

A = Y0 - C 

B = Y1 - A - C*EXF ( D > 

VALR = A + B*RAT 1 0 + C*EXP ( D*RAT 1 0 ) 

IF(FLAG) VALR = B + C*D*EXP < D *RATI 0 ) 

GO TO 400 

VALR = RATIO * (Y1 - Y0> + Y 0 
IF (FLAG) VALR = Y1 - Y0 

YM1 = AI MAG ( STORE ( 1 > ) 

Y0 = A I MAC ( STOR E ( 2 ) ) 

Y 1 = AIMAG ( STORE ( 3 ) ) 

Y2 = A I MAC ( STOR E ( 4 ) ) 

CURV1 = Y2 - 2 . 0* Y 1 + Y 0 

CURV2 = Y1 - 2 . 0*N 0 + YM1 

IF ( CUP.V1*CURV2 . LT.0.0) GO TO 570 

IF (ABS(CUFV1 ) .LE.0.0 .OR. ABS ( CURV2 ) . LE . 0.0) GO TO 570 
D = ALOG(CURVl/CURV2) 

I F ( D . EQ .0) GO TO 570 
C = CURV1 / ( ! EXP ( D ) - 1.0) ** 2) 

A = Y0 - C 

B = Y1 - A - C*EXP ! D ) 

VAL I = A + B*RAT 1 0 + C*EXP ( D *RAT 1 0 ) 

I F ( FLAG ) VALI = B + C*D*EXP ( RATIO*D > 

GO TO 500 

VALI - RATIO * ( Y 1 - Y0) + Y0 
IF (FLAG) VALI = Y1 - Y0 

TERP5 = CMP LX ( VALR , VAL I ) 

RETURN 

END 

SUBROUTINE MRT ( VOL T , CUR , NOP NTS . I BASE 1 , CHI .NOCHI , 

1 RN, FACTOR. FREQ ,Z IMPED, VSOURC , OL DVOL , FLAG. Z0.NOPER) 


S IS 1 7060 
SIS 17070 
S I S 1 7080 
S I S 1 7090 
SIS 17 100 
SIS171 10 
S I S 1 7 1 20 
SIS 17 130 
S I S 1 7 1 40 
S I S 1 7 1 50 
S I S 1 7 1 60 
SIS1 7170 
S I S 1 7 1 80 
S I S 1 7 1 90 
SIS 17200 
SIS17210 
S I S 1 7220 
SIS17230 
S I S 1 7240 
S IS 1 7250 
SIS 17260 
S IS 1 7270 
S IS 1 7280 
SIS 17290 
S I S 1 7 300 
S IS 1 73 10 
S IS 1 7320 
SIS17330 
S I S 1 7340 
S I S 1 7350 
SIS 17360 
SIS17370 
S IS 1 7380 
S I S 1 7390 
S I S 1 7 400 
SIS17410 
SIS17420 
S I S 1 7430 
S I S 1 7440 
S I S 1 7450 
S I S 1 74 60 
SIS17470 
S I S 1 7 480 
S I S 1 74 90 
SIS 17500 
SIS17510 
SIS 17520 
S I S 1 7530 
S I S 1 7540 
S I S 1 7550 
S IS 1 7560 
SIS17570 
S I S 1 7580 
S I S 1 7590 
S I S 1 7 600 



65.32 


FILE: SISCAP FORTRAN A NASA/GISS CONVERSATIONAL MONITOR SYSTEM SP302 


REAL VOLT ( NOPNTS ) , CUR(NOPNTS), CHI(NOCHI), OLDVOL ( NOPNTS > , 
1 ESOURC (128) 

COMPLEX FACTOR ( NOPNTS ) , Z 1 MP ED ! NOP NTS > , VSOURC ( NOPNTS ) , 
1ELEFTU28), ERIGHTU28), RHO(128) 

LOGICAL FLAG 
C 

C THIS SUBROUTINE CONTROLS THE MULTIPLE REFLECTION 

C LARGE SIGNAL ANALYSIS PROCEDURE 

C 

INTER = - I BASE 1 
NODIV2 = NOPNTS / 2 
C 

C INITIALIZE THE LEFT AND RIGHT PROPAGATING WAVES 

C 

DO 10 I = 1 , NOD I V2 
ICORR = NOPNTS -1+2 
ELEFT(I) = CMPLX<0. 0,0.0) 

ERIGHT(I) = 2.0 * < Z0/ ! Z IMP ED ! I )+Z0)*VSOURC! I > ) 

IF ( I .EQ. 1 ) GO TO 10 
ELEFT(ICOFR) = CON JG < EL EFT ! I ) > 

ERIGHT! ICORR) = CCNJG ( ER I GHT( I ) > 

10 CONTINUE 

C 

EL EFT ! NOD I V2 + 1) = CMP LX < 0 .0 , 0 . 0 ) 

EP IGHT ( NOD I V2 + 1) = CMP LX < 0 . 0 , 0 .0 ) 

C 

CALL FFT( INTER, ERIGHT, FACTOR) 

C 

DO 15 I = 1, NOPNTS 
ESOURC ( I ) = ERIGHT! I ) / Z0 
OLDVOL ( I ) = 0.0 
15 CONTINUE 

C 

C CALCULATE THE REFLECTION COEFFICIENT OF THE EMBEDDING 

C NETWORK AT EACH FREQUENCY OF INTEREST 

C 

DO 17 I = 1.NODIV2 

RHO(I) = (ZIMPED(I) - Z0> / (ZIMPED! I ) + Z0) 

ICORR = NOPNTS -1+2 

IF(I.EQ.l) GO TO 17 

RHO( ICORR ) = CONUG ! RHO! I ) ) 

17 CONTINUE 

R HO ( NOD I V2 + 1) = CMPLX(0. 0,0.0) 

C 

WR I TE ! 6 , 604 ) 

604 FORMAT! /IX , ’GENERAL NONLINEAR ANALYSIS USING MULTIPLE’, 

1’ REFLECTION TECHN I QUE ’ / 1 X , ’ BEG I NN I NG THE ITERATION PHASE’) 
C 

DO 1000 11 = 1 , NOP ER 
C 

WRITE (6, 400) II 

400 FORMAT! /IX ,’ ITERATION NUMBER = ’,15) 

C 

C CALCULATE THE TUNNEL JUNCTION VOLTAGE AND CURRENT WAVEFORMS 

C 


SIS17610 
S IS 1 7620 
S I S 1 7630 
SIS1 7640 
S I S 1 7650 
S IS 1 7660 
S IS 1 7670 
S I S 1 76 80 
S I S 1 7690 
S IS 1 7700 
SIS17710 
S IS 1 7720 
S I S 1 77 30 
S I S 1 7740 
S I S 1 77 50 
SIS1 7760 
S I S 1 77 70 
S I S 1 7 780 
S I S 1 7790 
SIS 17800 
S I S 1 7 8 1 0 
SIS 17820 
SIS1 7830 
S I S 1 7 8 40 
S I S 1 7 8 50 
S I S 1 78 60 
SIS17870 
S IS 1 7880 
S I S 1 78 90 
S I S 1 7 900 
SIS17910 
S IS 1 7920 
S I S 1 7930 
SIS 17940 
S I S 1 7 9 50 
S I S 1 7 960 
SIS 1 7970 
SIS17980 
S I S 1 7 990 
SIS 18000 
S IS 1 8010 
S I S 1 8020 
S I S 1 8030 
SIS 18040 
S I S 1 8050 
S I S 1 8060 
S I S 1 8070 
SIS 18080 
S I S 1 8090 
S 1 S 1 8 1 00 
SIS181 10 
SIS1 8120 
SIS18130 
SIS18140 
S I S 1 8 1 50 
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CALL SISC HA (ESOURC, NOPNTS, Z0, VOLT, CUR, CHI, NOCHI , RN , FREQ , I BASE 1 , 

SIS 18 160 


1NOPER, FACTOR) 

S I S 1 8 1 70 

c 


S IS 1 8 1 80 

c 

DETERMINES THE REVISED LEFT PROPAGATING WAVE 

SIS18190 

c 


S I S 1 8200 


DO 20 I = 1 , NOPNTS 

SIS18210 


ELEFT ( I ) = (VOLTII) - CUR(I) * Z0> * 0.5 

SIS 18220 

20 

CONTINUE 

S IS 1 8230 

C 


S IS 1 8240 


CALL FFT! IBASE1 , ELEFT, FACTOR) 

S I S 1 8250 

C 


S I S 1 8260 

C 

DETERMINE THE REVISED RIGHT PROPAGATING WAVE 

S I S 1 8270 

C 


S I S 1 8280 


DO 30 I = 1 , NODIV2 

S I S 1 8290 


ICORR = NOPNTS -1+2 

S I S 1 8300 


ERIGHT(I) = (Z0/(Z IMPED! I ) +Z0 ) *VSOURC ! I ) + RHO( I )*ELEFT( I ) )*2.0 

S I S 1 83 10 


IF! I .EQ. 1 > GO TO 30 

SIS18320 


ERIGHT! ICORR ) = CONOG ! ER IGHT! I > ) 

S I S 1 8330 

30 

CONTINUE 

S IS 1 8340 


ER IGHT! NODIV2 + 1) = CMP LX ( 0 . 0 , 0 . 0 > 

S IS 1 8350 

C 


S IS 1 8360 


CALL FFT! INTER, ERIGHT, FACTOR) 

S IS 1 8370 


DO 40 I = 1, NOPNTS 

S I S 1 8380 


ESOURC! I > = ERIGHT! I ) / Z0 

S IS 1 8390 

40 

CONTINUE 

SIS 18400 

C 


SIS18410 

C 

SEE IF A SOLUTION HAS BEEN OBTAINED 

S I S 1 84 20 

c 


S I S 1 8430 


CALL COMPAR ! VOLT, OL DVOL , NOP NTS .FLAG) 

S I S 1 8440 

c 


S I S 1 8450 


DC 500 00 = 1 .NOPNTS 

S IS 1 8460 

500 

OL DVOL ! 00 ) = VOLT ( 00 > 

SIS1 8470 

C 


S I S 1 84 80 


IF! .NOT. FLAG) GO TO 2000 

S IS 1 8490 

C 


S IS 1 8500 

1000 

CONTINUE 

S IS 1 85 10 

c 


SIS1 8520 


WP ITE!G, 1500) NOPER 

S I S 1 8530 

1500 

FORMAT! /8X ,’ LARGE SIGNAL ANALYSIS DID NOT CONVERGE IN ’ ,15,’ITNS* 

>SIS18540 


RETURN 

S I S 1 8550 

C 


S IS 1 8560 

2000 

WR ITE ( G , 2002 ) II 

SIS18570 

2002 

FORMATI/1X, ’LARGE SIGNAL ANALYSIS COMPLETED IN ’,15,' ITNS’/) 

S IS 1 8580 


RETURN 

S IS 1 8590 


END 

S IS 1 8600 


SUBROUTINE S I SCHA ( ESOURC , NOP NTS , RSOURC , VOLT , CUR , CH I , NOCH I , 

S IS 1 86 10 


1RN.FREQ, I BASE, NOPER, FACTOR) 

S I S 1 8620 

C 


S IS 1 8630 


REAL ESOURC! NOPNTS ) , VOLT ( NOP NTS ) , CUR(NOPNTS) 

S IS I 8640 


REAL CHI !NOCHI ) 

S I S 1 8650 


COMPLEX FACTOR (NOPNTS) 

S IS 1 8660 


COMMON /RBOUT / VDC , VLO , VO DC , VOL 0 , CURDC , CUR PMP , RLOSS , RTEMP 

S IS 1 8670 

C 


S I S 1 8 6 80 


REAL ZTERM ! 5 1 2 ) , PHASE(512), VOLT2I512) 

S I S 1 8 6 90 


COMPLEX UT ERM ( 5 1 2 ) , CCUR(128) 

S IS 1 8700 
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SET UP PARAMETERS FOR THE CALL TO SUBROUTINE 

SIS TO DETERMINE THE JUNCTION VOLTAGE AND CURRENT 

WAVEFORMS 

HDEL = 1.0 / FREQ / FLOAT! NOPNTS ) 

IADJST = 0 

CALL SIS (CHI , VOLT2, NOP NTS, HDEL, ZTERM, UTERM, 

1 RN.NOCH I , NOP ER, PHASE, NOSTOP, IADJST, RSOURC , ESOURC , I STO) 

DO 1000 I = 1 , NOPNTS 

I PNT = (ISTO - 1) * NOPNTS + I - IADJST 
VOLT( I ) = VOLTZ ( I P NT ) 

CUR ( I > = ESOURC ( I > - VOLT! I) / RSOURC 
CCUR ! I ) = CUR! I ) 

CONTINUE 

CALL FFT! IBASE, CCUR, FACTOR) 

CURPMP = CCUR! 1 > 

RETURN 

END 

SUBROUTINE SIS (CHI , VOLT , NOP NTS , HDEL , ZTERM , UTERM , 

1RN, NOCHI , NOP ER , P HASE , NOSTOP , IADJST , RSOURC , ESOURC , ISTO) 

INTEGER R, RR . VCNT 

REAL ZTERK! NOCHI ) , VOLT(NOCHI), CHI(NOCHI), PHASE(NOCHI) 
REAL ESOUPC!NOPNTS) 

COMPLEX UTERM! NOCHI > , SUM, CDUM, CONST, CSUM, CTEMP 
LOGICAL FLAG 

DATA HBAR , E/1.055E-34.1. 602 E- 19/ 

DATA TOL, I TYPE/ 1 . E-5 , 1 / 

DATA NHIGH/100/ 

DATA VTHRSH/0 0E-3/ 


DETERMINE THE SIS TUNNEL JUNCTION CURRENT AND 
VOLTAGE WAVEFORM COMPONENTS 

RNN = l./(l./RN + 1. /RSOURC) 

NOSTOP = NOCHI 
UTERM! 1 > = CMP LX ! 1 . ,0. ) 

VOLT! 1 > = 0. 

ZTERM ( 1 ) = 0. 

PHASE! 1 > = 0. 

CONST = CMPLX !0. , - 1 . ) * E / HBAR * HDEL 
C 

DO 1000 I = 1 , NOP ER 
ISTO = I 
VCNT = 0 
C 

DC 500 RR = 1. NOPNTS 
C 

N = NOPNTS * ( 1-1 ) + RR 
C 


SIS18710 
S IS 1 8720 
SIS 18730 
S I S 1 8740 
SIS 18750 
SIS18760 
S IS 1 8770 
S IS 1 8780 
S I S 1 8790 
S I S 1 8800 
S IS 1 88 10 
S IS 1 8820 
S IS 1 8830 
S I S 1 8840 
SIS 18850 
S I S 1 8 860 
S I S 1 8870 
S I S 1 8880' 
S I S 1 8890 
SIS 18900 
S I S 1 89 10 
S I S 1 8 9 20 
S IS 1 8930 
SIS18940 
S IS 1 8950 
S I S 1 8960 
S I S 1 8970 
S I S 1 8980 
S I S 1 8990 
S I S 1 9000 
S I S 1 9010 
SIS 19020 
S I S 1 9030 
S I S 1 9040 
S I S 1 9050 
S IS 1 9060 
S I S 1 9070 
S I S 1 9080 
S I S 1 909£ 
SIS 19 100 
SIS191 10 
S I S 1 9 1 20 
S I S 1 9 1 30 
S I S 1 9 1 40 
S I S 1 9 1 50 
S I S 1 9 1 60 
S I S 1 9 1 70 
S I S 1 9 1 80 
S I S 1 9 1 90 
S IS 1 9200 
SISI9210 
SIS19220 
S I S 1 9230 
S I S 1 9240 
S IS 1 9250 
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c 

IF NECESSARY. SHIFT ALL VALUES LEFT IN THE STORAGE VECTORS 

SIS19260 

c 



S IS 1 9270 


IFtN.LT. NOCHI) GO TO 658 


SIS19280 


DO 657 NPRO - 2. NOCHI 


SIS19290 


VOLT(NPRO-l) = VOLT ( NPRO ) 


S IS 1 9300 


PHASE! NPRO-1 ) = PHASE(NPRO) 


SIS19310 


UTERM ! NPRO- 1 ) = UTERM(NPRO) 


S IS 1 9320 


ZTERM! NPRO-1 ) = 2TERM ( NPRO ) 


SIS’ 9330 

657 

CONTINUE 


S IS 1 9340 

C 



S I S 1 9350 


IADJST = IADJST + 1 


S I S 1 93 60 

C 



S I S 1 9370 

C 

INITIALLY ESTIMATE THE NEXT VOLTAGE POINT BY 

SIS1 9380 

C 

ASSUMING A CONSTANT VOLTAGE 


S I S 1 9390 

C 



SIS 19400 

658 

VOLT(N-IADOST+l ) = VOLT! N- IADOST ) 


SIS19410 


PHASE (N-IADJST+1 ) = PHASE < N- IADJST 

> + AIMAG(CONST) * 

S I S 1 9420 


1VOLT(N-IADOST+1 ) 


S IS 1 9430 


C/>LL FIXPHA(P HASE (N-IADJST+1 ) ) 


S IS 1 9440 


CTEMP = CMPLX (0. , PHASE (N-IADJST+1 ) 

) 

SIS1 9450 


UTERM(N-1ADJST+1 ) = CEXP(CTEMP) 


S I S 1 94 60 


VFREV = VOLT( N- IADOST+ 1 > 


S I S 1 9470 


VMID1 = VOLT(N-IADJST+l ) 


S I S 1 9 480 


VNID2 = VOLT(N-IADOST+1 ) 


S I S 1 9490 

C 



SIS 19500 


CSUM = CMPLXC0. .0. ) 


SIS19510 

C 



SIS 19526 

C 

ITERATE TO DETERMINE THE CORRECT VOLTAGE 

S IS 1 9530 

C 

AT THE NEXT POINT IN TIME 


S IS 1 9540 

C 



S IS 1 9550 


DO 425 MM = l.NHIGH 


S I S 1 9560 

C 



S I S 1 9570 


SUM = CMPLX! 0 . ,0. > 


S I S 1 9580 

c 



S IS 1 9590 


I F ( MM . NE . 1 .OP. N.EQ.I) GO TO 243 


S IS 1 9600 

c 



S I S 1 9 6 1 0 

c 

PERFORM THE REQUIRED CONVOLUTION 


SIS1962*. 

c 



SIS1 9630 


NMIN1 = N - 1 


S I S 1 9640 


ITPRY = IADOST + 1 


S IS 1 9650 


DO 100 R = ITPRY, NMIN1 


S IS 1 9660 


CSUM = (CHI (N-R+l ) *UTERM! R- IADJST + I > + CH I ( N-R+2 > *UTERM ( R- I ADOST > > S I S 1 9 6 70 


1/2. + CSUM 


S I S 1 9680 

100 

CONTINUE 


SIS1 9690 

c 



S IS 1 9700 

243 

SUM = { C H I ( I >*UTERM(N-IADUST+1 ) + 

CHI (2)*UTERM( N- IADOST) >/2. 

S I S 1 97 10 


1+ CSUM 


S I S 1 9720 

C 



S IS 1 9730 


IF {MM.NE. 1 ) GO TO 101 


S I S 1 9740 


YF = YFUNCT(N.HDEL , CHI .NOCHI ) 


S IS 1 9750 


N2 = ( N+ 1 ) - (N/NOPNTS*NOPNTS) 


S IS 1 9760 


CC = ESOUPC(NZ) 


S I S 1 9770 

101 

CDUM = COf'JG!UTERM (N-IADJST+1 > ) * 

SUM * HDEL + 

S I S 1 9780 


1CONOG(UTERM(N-IADOST+1 ) ) * YF 


S I S 1 9790 


ZTERM( N-IADJST+1 ) = -CC + AIMAG(CDUM) 

S IS 1 9800 
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C 

C 

C 

c 

c 

c 


c 

423 


C 

421 


C 

422 


425 

C 


234 

C 

C 

C 

450 


VTST = < (-1 ./CAP*(ZT£RM(N-IADJST+1 >+ZTERM(N-IADJST> )/2. 
1-1 ./2./CAP/RNN*VOLT!N-IADJST> ) * 

1HDEL + VOLT(N-IADOST) )/( 1 .+HDEL/2./RNN/CAP) 

VTST = -1.0 * RNN * ZTERM(N-IAD0ST+1> 

ITS2 = N - IADJST + 1 

IF (ABS(VTST-VOLT( ITS2) > . LT . TOL *ABS { VOLT ( ITS2) ) ) 

1 GO TO 450 

IF <MM/2*2. EQ.MM) GO TO 421 

VOLT(N-IADJST+l ) = (VTST + VOLT! N- IADJST+ 1 ) ) /2 . 
PHASEIN-IADJST+1 > = P HASE ( N- I ADOST > + AIMAG ( CONST ) * 

1! VOLT (N-IADJST+1 ) +VOLT < N- I ADOST > >/2. 

CALL F IXPHA( PHASE! N-IADJST+1 ) > 

CTEMP = CMPLX ( 0. , P HASE (N-1AD0ST+1 ) ) 

UTERM!N-IADJST+1 > = CEXP(CTEMP) 

GO TO 422 

IF (ABS(VMID2-VPREV) .LT. TOL*ABS ( VMI D2 ) ) GO TO 423 
VG = (VTST - VMID1) / (VMI02 - VPREV) 

I F ( ABS ( 1 . -VG > .LT. TOL) GO TO 423 

VOLT(N-IADJST+l ) = (VMID1 - VG*VPREV) / (1. - VG > 

P HASE (N-IADOST+1 ) = PHASE ( N- IADJST > + AIMAG(CONST) * 
l(VOLT(N-IADJST+l ) + VOLT < N- I ADOST )> /2 . 

CALL F I XP hA ( P HASE ( N-IADJST+1 ) ) 

CTEMP = CMPLX ( 0. , PHASE ! N-IADJST+1 > ) 

UTERM(N-IADJST+1 > = CEXP(CTEMP) 

VMID1 = VTST 

VPREV = VMID2 

VMID2 = VOLT! N-IADOST+1) 

CONTINUE 

WR ITE ( 6 , 234 ) 

FORMAT! IX ******PROBL EM : ITERATION LOOP 1 EXPIRED’) 

STOP 

ITS10 = N - IADJST + 1 

IF (VOLT! ITS10) .LT. VTHRSH ) VCNT = VCNT + 1 


C 

I F ! RR . EQ . NOPNTS ) GO TO 600 
500 CONTINUE 
C 

C A WATCH IS KEPT ON HOW MANY VOLTAGE POINTS IN 

C THE WAVEFORMS ARE NEGATIVE VOLTAGES 

C 

600 I F ( VCNT .FQ. 0) GO TO 620 
WR I TE ( 6 , 6jC 1 ) VCNT, I 

601 FORMAT! IX, ’SIS:VOLTAGE THRESHOLD COUNTER = ’,15,’ITN NO = 
C 

620 IF ( I .LE . ITYPE ) GO TO 1000 

CALL COMP 2 ! VOLT, I , NOPNTS , FLAG , ITYPE , NOCHI . IADJST) 


S IS 1 98 10 
S IS 1 9820 
S I S 1 9830 
SIS19840 
S I S 1 9850 
S I S 1 986£ 
SIS19870 
S I S 1 9880 
S I S 1 9 890 
S I S 1 9900 
SIS19910 
S I S 1 9920 
S I S 1 9930 
S IS 1 9940 
S I S 1 9 950 
S I S 1 9960 
S I S 1 9970 
S I S 1 9980 
S I S 1 9990 
SIS20000 
S IS20010 
SIS20020 
SIS20030 
SIS20040 
S IS20050 
SIS20060 
SIS20070 
SIS20080 
SIS20090 
SIS20100 
SIS201 10 
SIS20120 
SIS20130 
SIS2014P 
SIS20150 
SIS20160 
SIS20170 
SIS20180 
SIS20190 
SIS20200 
SIS20210 
S IS20220 
SIS20230 
SIS20240 
SIS20250 
SIS20260 
SIS20270 
S IS20280 
SIS20290 
SIS20300 
SIS20310 
’,15) S IS20320 
SIS20330 
SIS20340 
SIS20350 
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IF (FLAG) GO TO 2000 
1000 CONTINUE 
C 

WP. I T E ( 6 , 2 3 3 ) 

233 FORMAT! IX, ’******PROBLEM: ITERATION LOOP 2 EXPIRED’) 
NOSTOP = N + 1 
RETURN 
C 

2000 WRITE! 6,65 5) I TYPE, ISTO 

655 FORMAT? IX, ’ S I S : P ER I OD I C I TY OBTAINED WITH I TYPE = ’,15, 

1/1X, ’SISsNO. OF ITERATIONS = ’,15) 

NOSTOP = N + 1 

RETURN 

END 

FUNCTION YFUNCT! N.HDEL , CHI , NOCHI > 

C 

REAL CHKNOCHI ) 

C 

C THIS FORMS PART OF THE CALCULATION OF THE CONVOLUTION 

C FORMULA IN THE SIS EQUATIONS 

C 

YFUNCT = 0.0 
IF(N. GE. NOCHI ) RETURN 
C 

NPLUS1 = N + 1 

NKIN1 = NOCHI - 1 

DO 20 I = NPLUS1 , NMIN1 

YFUNCT = YFUNCT + (CHI ( I )+CHI ? 1 + 1 ) )/2.0*HDEL 
20 CONTINUE 

RETURN 
C 

END 

SUBROUTINE COMP2 ( VOLT , I , NOP NTS , FLAG , I TYPE , NOCH I , I ADUST) 
LOGICAL FLAG 
REAL VOLT ( NOCh I > 

DATA TOL, EPS/1. E-4.1.E-1/ 

C 

C THIS SUBROUTINE COMPARES THE TWO MOST RECENT 

C SUCCESSIVE VOLTAGE WAVEFORMS IN THE SIS WAVEFORM 

C ITERATION PROCESS 

C 

EMAX = 0.0 
I P NT = 0 
FLAG = .TRUE. 

IF! I .LE. I TYPE) GO TO 1000 
C 

CMAX = 0. 

DO 100 M = 1, NOP NTS 

M!' = (I - 1) * NOP NTS + M - IADJST 

IF (ABS! VOLT(Mt:) ) .GT. CMAX ) CMAX = ABS ( VOLT ( MM ) ) 

100 CONTINUE 
C 

I F ( CMAX .LT. l.E-12) RETURN 
C 

= 1, NOP NTS 


S IS20360 
S IS20370 
SIS20380 
SIS20390 
SIS20400 
SIS20410 
S IS20420 
SIS20430 
SIS2044C 
SIS20450 
SIS20460 
S IS20470 
SIS20480 
S I S 204 90 
S IS 20500 
SIS20510 
S IS<-0520 
S ISi.0530 
SIS20540 
SIS2055P 
S IS20560 
SIS20570 
SIS20580 
SIS20590 
SIS20600 
SIS20610 
S I S 20620 
S IS20630 
S IS20640 
S IS20650 
S IS20660 
S IS20670 
S IS20680 
SIS20690 
S IS20700 
SIS20710 
SIS20720 
S IS20730 
SIS20740 
S IS20750 
SIS20760 
SIS20770 
S1S20780 
SIS20790 
SIS20800 
SIS20810 
SIS 20820 
SIS20830 
S IS20840 
SIS20850 
S IS20860 
S IS20870 
SIS20880 
SIS2089C 
SIS20900 


DO 500 U 
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IPOINT = ( I - 1 ) * NOPNTS + J - IADUST 
ISEEK = IPOINT - NOPNTS * I TYPE 
I F < ISEEK . LE .0 > GO TO 500 

IF(ABS(VOLT( IPOINT) ) .LT. EPS*CMAX > GO TO 500 
ERROR = (VOLT( I PO I NT ) -VOLT l ISEEK) >/VOLT( IPOINT) 
ERROR = ABS { ERROR ) 

IF (ERROR .LT. EMAX) GO TO 400 
EMAX = ERROR 
IP NT = 0 

400 IF (ERROR .LT. TOL > GO TO 500 
FLAG = .FALSE. 

500 CONTINUE 
RETURN 
C 

1000 FLAG = .FALSE. 

RETURN 
END 

SUBROUTINE FIXPHA(X) 

DATA PI/3. 1415925/ 

THIS PERFORMS A MODULO FUNCTION ON THE PHASE 
CALCULATIONS IN THE SIS EQUATIONS 

IF (ABS(X> .LT.2.0*PI ) RETURN 
NOREV = X / 2.0 / PI 
X = X - 2.0 * PI * FLOAT ( NOREV ) 

RETURN 
END 


SIS20910 
S IS20920 
SIS20930 
SIS20940 
SIS20950 
S IS209G0 
S I S 20970 
SIS20980 
S I S 20990 
SIS21000 
SIS21010 
SIS21020 
SIS21030 
SIS21040 
SIS21050 
SIS21060 
S I S2 1070 
SIS21080 
S I S2 1090 
SIS21 100 
SIS21 1 10 
SIS21 120 
S1S21 130 
SIS21140 
SIS21 150 
SIS21 160 
SIS21 170 
SIS21 180 
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